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PREFACE 


This  report  and  the  model  it  describes  are  part  of  the  ongoing  program  in 
automated  fare  collection  and  transit  revenue  research  supported  by  the 
Transportation  Systems  Center  (TSC)  and  the  Urban  Mass  Transportation 
Administration  (UMTA) . The  work  was  sponsored  by  the  Office  of  Systems 
Engineering,  within  UMTA’s  Office  of  Technical  Assistance. 

The  need  for  developing  a reliability-based  performance  and  cost  analysis 
technique  was  identified  in  conjunction  with  the  American  Public  Transit 
Association  (APTA)  Fare  Collection  Reliability  Liaison  Board,  which  consists 
of  representatives  from  the  rail  transit  systems  within  APTA  and  provides 
guidance  to  the  TSC  Transit  Revenue  program.  The  objective  of  this  study  is 
to  develop  a technique  which  can  be  used  by  transit  systems  in  specifying  and 
procuring  fare  collection  equipment. 

The  author  wishes  to  acknowledge  the  generosity  of  the  Miami  Dade  County 
Transit  Authority  and  the  Metropolitan  Atlanta  Rapid  Transit  Authority  in 
making  available  specification,  performance,  and  cost  data  for  the  sample 
analyses  in  this  report.  While  the  analyses  themselves  do  not  necessarily 
represent  either  MDCTA  or  MARTA,  they  are  designed  to  be,  thanks  to  the 
supplied  data,  representative  of  situations  which  arise  in  actual  practice. 
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EXECUTIVE  SUMMARY 


The  collection  of  transit  system  fares  has  been  receiving  increased  attention  as 
fares  rise  and  Federal  operating  subsidies  decrease.  Transit  authorities  are 
becoming  more  concerned  about  ways  to  maximize  revenue  and  minimize  costs  while 
providing  equitable  fares  and  reliable,  convenient  service  for  passengers.  Fare 
collection  methods  have  a significant  impact  on  total  transit  costs,  amount  of 
revenue  generated,  and  passenger  service.  One  study  showed  that  fare  collection 
costs  range  from  1%  to  31/S  of  passenger  revenue  at  rail  transit  systems. 
Moreover,  revenues  generated  from  fares  can  vary  from  40^6  to  90/S  of  total 
transit  costs.  Fare  collection  systems  must  therefore  be  selected  only  after 
careful  examination  of  their  cost,  revenue,  and  service  effects. 

Automated  fare  collection  offers  the  potential  for  reducing  costs  by  minimizing 
the  need  for  personnel  to  perform  cumbersome,  repetitive  functions.  However, 
the  newer  and  more  complex  a piece  of  equipment,  the  more  likely  it  is  to  have 
frequent  failures,  which  can  lead  to  significant  passenger  delay,  lower 
throughput  capacity,  and  general  frustration.  Efforts  are  underway  to  increase 
the  reliability  of  automated  fare  collection  equipment.  However,  it  becomes 
imperative  to  know  just  how  much  of  a reliability  increase  is  required,  as  the 
amount  of  reliability  increase  contemplated  makes  a significant  difference  in 
the  cost  of  such  an  effort,  as  well  as  its  likelihood  of  success. 

Furthermore,  because  of  the  significant  cost  of  fare  collection  equipment,  the 
number  of  equipment  units  to  acquire  for  a given  station  also  becomes  quite 
important.  Too  many  units  can  increase  total  system  cost  considerably,  while 
too  few  can  lead  to  significant  passenger  congestion  and  delay  problems.  In 
addition,  it  is  important  to  reasonably  assess  total  system  costs,  and  thus 
control  them  by  comparing  the  costs  and  passenger  performance  for  various 
possible  system  specifications. 

This  report  describes  an  analysis  technique  designed  to  help  transit  systems  to 
make  more  effective  investment  decisions  in  selecting  fare  collection  methods, 
systems,  and  equipment  to  best  fit  their  needs,  minimize  costs,  and  provide 
equitable  and  convenient  service  to  passengers.  Software  implementing  this 
analysis  technique  is  available  in  the  form  of  a user-friendly  computer  model. 
The  model  is  designed  to  be  used  by  existing  transit  systems,  as  they  enhance  or 
expand  their  fare  collection  systems,  as  well  as  by  new  transit  systems,  as  they 
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consider  their  fare  collection  needs.  In  either  case,  the  model  can  uncover 
important  cost-saving  or  service-enhancing  implications  in  station  design. 

The  model  allows  transit  systems  to  assess  the  performance  and  cost  of  various 
fare  collection  systems,  including  those  with  entry  processing  only,  those  with 
both  entry  and  exit  processing,  and  barrier-free  systems.  A transit  authority, 
consultant,  or  analyst  can  use  this  software  to  examine  the  fare  collection 
system  in  various  useful  ways: 

o (Fare  collection  system  evaluation)  — An  existing  or  proposed  fare 
collection  system  can  be  evaluated  to  see  if  current  or  proposed 
specifications  on  sizing,  reliability,  and  other  features  will  lead  to 
acceptable  performance  and  cost. 

o (Sensitivity  analysis)  — Proposed  changes  in  specifications  can  be 
analyzed  to  determine  their  impact  on  cost  and  performance,  in 
particular  whether  the  changes  will  bring  (or  keep)  the  cost  and 
performance  within  acceptable  levels. 

o (Specification  determination)  — Specifications  necessary  to  meet  given 
cost  and  performance  constraints  can  be  determined. 

o (Tradeoff  analysis)  — Tradeoffs  can  be  made  between  different 

specifications,  tightening  one  specification  while  relaxing  another; 
between  various  types  of  costs,  decreasing  one  cost  while  increasing 
another;  or  between  cost  and  performance. 

Such  analyses  can  uncover  very  useful  findings.  For  example,  sample  analyses 
shown  in  this  report,  conducted  on  a derived  fare  collection  system  based  on 
data  from  actual  systems,  found  the  following: 

o The  equipment  reliability  specifications  could  be  sharply  reduced,  by 
nearly  an  order  of  magnitude,  without  significantly  affecting 
performance. 

o A proposed  minor  reliability  improvement  at  some  additional  acquisition 
cost  turned  out  to  decrease  corrective  maintenance  costs  to  the  extent 
that  the  total  costs  decreased,  while  performance  meanwhile  improved. 
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o Unless  the  equipment  reliability  is  very  poor,  a spares  margin  of  only 
one  additional  unit  (above  the  minimum  number  necessary  to  process  the 
average  peak-hour  arrival  rate  of  passengers  in  the  absence  of  failures) 
was  enough  to  yield  acceptable  performance. 

The  fare  collection  analysis  software  can  thus  provide  useful  insights  to  a 
variety  of  purposes,  such  as: 

Determination  of  the  number  of  machine  units  to  deploy  at  a station. 
Determination  of  reliability  and  maintainability  specifications  for  fare 
collection  equipment. 

Assessing  the  impact  of  changes  in  passenger  demand. 

Evaluating  changes  in  maintenance  policies. 

Evaluating  changes  in  fare  collection  procedures. 

The  technical  approach  of  the  software  is  to  model  the  operation  of  the  fare 
collection  system  as  a multiple-server  queue,  with  passengers  as  customers, 
machine  units  as  servers,  and  on  a first-come-first-served  service  discipline. 

A key  facet  of  the  approach,  and  the  reason  that  off-the-shelf  queuing  models 
cannot  be  used,  is  that  the  number  of  servers  (machine  units)  changes  as  the 
machine  units  fail  and  are  repaired. 

The  performance  measures  consist  of  passenger  congestion  (queue  length)  and 
passenger  delay  time.  These  are  expressed  as  average  values,  as  well  as 
probability  frequency  distributions.  The  cost  measures  consist  of  equipment 
acquisition  costs,  spares  costs,  equipment  operating  costs,  scheduled- 
maintenance  costs,  and  corrective-maintenance  costs,  all  of  which  are  computed 
on  an  annualized  basis. 

Three  kinds  of  input  data  are  required  for  the  fare-collection  analysis; 
hardware  data,  passenger  flow  data,  and  cost  data  (if  cost  analyses  are 
desired).  The  hardware  and  passenger-flow  data  required  are: 

Passenger  arrival  rate 
Group  size  (optional) 

Passenger  processing  rates 
Failure  rates 
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Repair  times 

Number  of  machine  units 

Division  of  passenger  flow  (for  multiple  service  areas  only) 

For  cost  analyses,  the  required  input  data  are: 

Acquisition  cost  per  unit 

Useful  life  of  the  unit 

Discount  rate 

Spares  ratio 

Operating  cost  per  unit 

Annual  hours  of  scheduled  maintenance 

Hourly  pay  rate  for  repair  personnel 

Annual  passenger  volume  at  the  fare  collection  area 

To  demonstrate  the  use  of  the  software,  sample  fare  collection  system 
assessments  are  described.  These  analyses  are  based  on  passenger  demand, 
equipment  performance,  reliability  and  maintainability,  and  cost  data  from 
actual  transit  systems.  Passenger-performance  and  cost  results  are  obtained, 
and  fare  collection  system  evaluations,  sensitivity  analyses,  specification 
determinations,  and  tradeoff  analyses  based  on  these  results  are  described. 
Conclusions  of  these  analyses  are  given.  Among  the  outputs  of  the  analysis  are 
a sensitivity  graph  showing  the  effect  of  equipment  reliability  on  passenger 
delay,  and  an  equivalent-cost  tradeoff  graph,  which  shows,  for  various  increases 
in  equipment  reliability,  by  how  much  the  resulting  equipment  acquistion  cost 
can  increase  without  increasing  the  overall  annual  system  cost. 

Further  information  on  the  analysis  software,  including  source  code,  is 
available  from  the  author  or  the  Transportation  Systems  Center  upon  request. 
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1 . INTRODUCTION 


The  collection  of  transit  system  fares  has  become  more  sophisticated  in 
recent  years,  as  transit  authorities  turn  to  more  flexible  fare  structures. 
Instead  of  a single  flat  fare  to  use  the  system,  the  fare  now  often  depends 
on  the  passenger's  origin  and  destination  points.  This  is  done  for  various 
reasons,  such  as  to  make  the  fares  more  equitable  in  relation  to  the  actual 
distance  travelled  on  the  system,  to  allow  for  fare  advantages  to  certain 
areas  or  to  certain  types  of  passengers  (such  as  senior  citizens  or 
students),  or  to  allow  for  accommodation  among  the  various  governmental 
bodies  and  other  groups  which  support  the  system  (1,2, 3, 4)* 

In  order  to  collect  the  fares  on  a system  with  such  a complex  fare  structure, 
a transit  property  can  no  longer  use  a simple  coin-  or  token-activated 
turns tyle.  More  sophisticated  equipment  is  necessary  to  ascertain  boarding 
and  exiting  information  for  each  passenger  and  to  encode  it  on  tickets  for 
the  necessary  processing,  to  vend  tickets  and  collect  money,  etc.  This  more 
sophisticated  function  could  conceivably  be  done  by  human  beings.  However, 
in  most  systems  to  use  personnel  to  the  extent  that  would  be  necessary  would 
be  too  costly.  Therefore,  transit  systems  have  turned  to  more  sophisticated 
fare  collection  machinery,  using  data  processing  and  electronics,  to  operate 
the  fare  collection  system  (1,2,4). 

However,  a major  problem  arises  in  the  use  of  this  equipment,  that  of 
reliability.  Quite  often,  the  more  complex  a piece  of  equipment,  and  the 
newer  its  technology,  the  more  likely  it  is  to  have  relatively  frequent 
failures  (due  to  some  extent  to  electronic  complexity,  but  mostly  due  to 
mechanical  aspects  such  as  ticket  transport,  coin  and  bill  detection  and 
transport,  etc.).  Transit  systems  using  this  equipment  have  indeed  often 
experienced  high  failure  rates,  leading  to  significant  passenger  delay, 
passenger  pass-throughs  (i.e.,  the  emergency  gates  are  opened  and  passengers 
are  allowed  into  the  system  with  the  fare  either  collected  manually  or  not  at 
all),  lower  throughput  capacity,  and  general  frustration  (5,6).  Efforts  are 
underway  to  increase  the  reliability  of  fare  collection  equipment,  through 
better  off-the-shelf  components,  more  reliable  components,  redundancy 
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techniques,  or  improved  maintenance  procedures  (5, 6, 7, 8).  The  question  that 
arises,  however,  is  by  just  how  much  should  the  reliability  be  improved? 

While  improvements  in  reliability  are  necessary,  they  are  also  expensive, 
both  in  terms  of  money  and  time.  Under  some  circumstances,  such  as  the  ones 
shown  in  the  following  examples,  equipment  improvements  may  not  be  cost- 
effective: 

o The  wrong  service  area  is  being  improved:  Quite  clearly,  one  should 

not  go  to  any  expense  to  improve  the  reliability  of  a subsystem  whose 
failures  are  not  strongly  impeding  the  system  (e.g.,  if  most  of  the 
delay  occurs  at  the  ticket  vendors,  don't  improve  the  gates).  However, 
the  process  leading  to  a delay  can  be  subtle.  For  example,  suppose  most 
of  the  delay  in  a fare  collection  system  presently  occurs  at  the  ticket 
vendors,  and  occurs  there  because  the  ticket-vendor  reliability  is  poor. 
One  might  expect  that  system  delay  could  be  reduced  by  improving  the 
ticket-vendor  reliability.  Suppose,  however,  that  all  passengers  must 
pass  through  the  gates,  and  that  the  gates,  while  not  showing 
significant  delays,  are  presently  processing  just  about  as  many 
passengers  as  they  can  handle.  If  the  reliability  of  the  vendors  is 
improved,  what  may  happen  is  that  the  passenger  flow  to  the  gates  may 
increase  beyond  their  capacity  to  process  it,  causing  large  delays  at 
the  gate  area.  In  other  words,  the  effort  and  expense  to  improve  the 
vendor  reliability  would  be  ineffective  in  this  case,  since  it  would 
just  shift  the  delay  from  the  ticket  vendors  to  the  gates. 

o The  reliability  is  improved  too  much:  As  the  reliability  of  a 

subsystem  improves,  its  frequency  of  failure  of  course  decreases.  After 
some  point,  it  no  longer  falls  often  enough  for  the  failures  to 
significantly  affect  system  operation.  Any  expense  spent  to  further 
improve  reliability  beyond  this  point,  even  if  successful,  is  a wasted 
expense. 

o Measures  other  than  reliability  improvement  may  be  more  efficient: 
Improving  the  equipment  reliability  is  not  the  only  way  to  reduce  the 
impact  of  failures  on  system  performance.  Faster  recovery  times  in  case 
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of  failures  (i.e.,  maintainability  or  recoverability)  or  having  more 
units  available  for  service  (i.e.,  redundancy)  may  improve  system 
performance  even  with  the  same  reliability.  While  these  measures  also 
have  their  costs,  they  may  be  less  expensive  than  reliability 
improvements,  particularly  if  less  of  an  improvement  is  needed  through  a 
maintainability  or  redundancy  approach. 

o System  failure  may  not  be  the  main  problem:  Though  there  may  be 

congestion  and  delay  in  a system,  and  though  there  may  be  failures,  the 
failures  may  not  be  the  primary  cause  of  the  delay.  For  example,  if  a 
station  is  served  by  a major  feeder  bus  line  which  periodically  sends 
large  numbers  of  passengers  simultaneously  into  the  fare  collection 
system,  delays  will  occur  which  will  not  be  primarily  due  to  failures. 

If  so,  reducing  the  failures  that  exist  will  not  really  help  the 
problem. 

In  order  to  properly  answer  the  question,  "By  how  much  should  reliability  be 
improved?",  one  needs  some  way  to  find  out  the  passenger  delay  in  a fare 
collection  system,  given  information  on  its  reliability,  maintainability,  the 
number  of  machine  units  (redundancy),  nominal  processing  rate,  and  passenger 
demand.  With  such  a method  one  can  investigate  the  effect  on  delay  of  not 
only  reliability  improvements,  but  improvements  in  maintainability,  operating 
policy,  and  the  number  of  machine  units  as  well,  and  derive  the  proper  mix 
and  extent  of  improvements  necessary. 

This  report  describes  a method  we  have  developed  to  analyze  this 
interrelationship  among  reliability,  maintainability,  the  number  of  machine 
units,  and  passenger  delay.  The  basic  approach  is  a model  to  simulate  the 
flow  of  passengers  through  the  fare  collection  system.  This  is  done  by 
treating  the  system  as  a network  of  queues,  with  a queue  at  each  service  area 
and  the  passengers  moving  from  one  service  area  to  the  next  (a  "service  area" 
is  a specific  set  of  machine  units,  such  as  coin  and  bill  changers,  ticket 
vendors,  gates,  etc.).  Superimposed  on  this  network  is  the  failure/recovery 
process,  by  which  units  fail  at  a rate  according  to  their  reliability  and  are 
repaired  according  to  the  failure  response  and  repair  times. 
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In  addition,  we  have  developed  an  analytical  (mathematical)  model  to  evaluate 
passenger  congestion  and  delay  at  a single  service  area,  given  failures  and 
repairs.  The  analytical  model  obtains  the  congestion  and  delay  values  not  by 
simulation,  but  rather  by  solving  the  underlying  queue  length  probability 
equations  directly.  The  model  investigates  a single  service  area,  with 
certain  assumptions  on  the  probability  distributions  being  needed  in  order  to 
solve  the  resulting  equations,  which  somewhat  limits  the  scope  of  the  model. 
However,  within  that  scope,  the  model  requires  much  less  computer  time  for  a 
large-scale  analysis  than  does  a simulation,  and  also  provides  the  actual 
underlying  queue  length  distribution,  which  a simulation  provides  only  after 
running  for  a large  amount  of  simulated  time.  With  the  analytical  model, 
therefore,  many  analyses  can  be  carried  out  with  the  computer  resources 
necessary  for  one  simulation  run,  and  so  the  analytical  model  should  be  used 
in  the  situations  for  which  it  is  applicable.  One  such  manner  is  to  use  the 
analytical  model  to  conduct  the  initial  investigations,  in  which  many  runs 
may  bemade,  which  many  different  parameter  values.  From  the  results  of  these 
investigations,  if  desired,  a smaller  number  of  scenarios  of  interst  can  be 
identified  and  analyzed  further  by  means  of  the  simulation  model. 

Section  2 describes  the  types  of  results  which  the  models  provide,  both 
direct  results  such  as  passenger  congestion  and  delay,  and  indirect  results 
such  as  sensitivity  and  tradeoff  analyses.  It  also  describes  how  transit 
systems  can  make  use  of  these  model  results  for  planning,  procurement, 
maintenance,  or  analyses  of  operating  policies. 

Section  3 develops  the  detailed  technical  approach  of  the  models.  The 
concept  of  the  fare  collection  system  as  a network  of  queues  is  outlined. 

The  simulation  model  is  then  described,  both  for  a single  service  area  (e.g., 
the  ticket  vendors)  and  for  entire  fare  collection  systems.  This  is  followed 
by  a description  of  the  analytical  model. 

Section  4 describes  the  data  required  for  fare  collection  system  analysis, 
procedures  to  collect  them,  and  possible  problems  and  cautions.  To 
demonstrate  how  the  modeling  process  operates  and  the  sensitivity  and 
tradeoff  analyses  which  result,  several  examples  are  provided  in  Section  5, 
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including  one  based  on  an  actual  transit  system  currently  in  development. 

The  final  section,  Section  6,  outlines  further  work  to  verify  and  improve  the 
models  and  to  produce  a user-oriented  package  for  transit  systems. 
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2.  THE  ANALYSIS  PROCESS 


2.1  Types  of  Dependability  Analyses 

The  direct  output  of  the  fare  collection  dependability  model  is  information 
on  the  congestion  (queue  length)  and  passenger  delay  in  the  fare  collection 
system,  given  the  system  configuration  and  passenger  demand.  This  is 
produced  for  each  of  the  service  areas  (i.e.,  ticket  vendors,  gates,  etc.)  in 
the  system,  as  well  as  delay  for  the  overall  system.  The  information  is 
given  in  terms  of  probability  frequency  distributions  for  congestion  and  for 
delay  at  each  service  area,  from  which  means  and  variances  are  obtained. 

This  output  can  give  rise  to  four  different  kinds  of  analyses: 

Evaluation 

Sensitivity  Analysis 
Specification  Determination 
Tradeoff  Analysis 

In  evaluation,  a given  fare  collection  system  is  examined.  The  required 
information  about  the  system  is  collected  and  entered  into  the  model  as  input 
data.  The  results  are  an  estimate  as  to  how  well  the  system  performs  under 
the  given  passenger  demand,  reliability,  maintainability,  and  number  of 
machine  units. 

Once  a system  has  been  evaluated,  one  may  naturally  wish  to  know  how  the 
results  would  change  if  one  or  more  of  the  input  parameters  were  different 
from  their  indicated  values  (this  is  particularly  true  if  there  is  some  doubt 
as  to  the  values  of  some  of  the  input  parameters).  One  can  then  make  several 
runs  of  the  model  with  differing  values  for  a given  input  parameter,  and  see 
what  changes  occur  to  system  congestion  and  delay.  This  is  called 
sensitivity  analysis,  as  it  measures  the  sensitivity  of  congestion  and  delay 
to  changes  in  input  parameters  such  as  reliability,  maintainability,  or 
number  of  machine  units. 
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Specification  determination  is  a process  which  determines  the  values  of  given 
input  parameters  (such  as  reliability  or  maintainability)  necessary  to 
achieve  a desired  level  of  performance.  These  values  then  become 
specifications  for  those  parameters.  Specification  determination  is  the 
reverse  of  sensitivity  analysis,  in  that  while  sensitivity  analysis  takes 
given  values  of  the  input  parameters  and  determines  the  impact  of  these 
values  on  system  performance,  specification  determination  takes  a given  level 
of  system  performance  and  determines  the  input  parameter  values  which  are 
necessary  to  achieve  it. 

Another  aspect  one  may  wish  to  know  about  the  system  is  tradeoffs,  i.e.,  how 
could  one  trade  off  between  two  input  parameters  while  achieving  the  same 
overall  result.  For  example,  if  the  equipment  reliability  declined,  by  how 
much  would  maintainability  have  to  improve  to  obtain  the  same  overall 
performance?  This  type  of  analysis,  called  tradeoff  analysis,  differs  from 
sensitivity  analysis  in  that  tradeoff  analysis  examines  the  interaction 
between  two  input  parameters,  whereas  sensitivity  analysis  examines  the 
interaction  between  an  input  and  an  output  parameter. 

The  relationship  among  these  four  different  types  of  analyses  is  shown  in  the 
conceptual  graph  given  in  Figure  1 . 

2.2  Possible  Uses  of  Dependability  Analyses  by  Transit  Systems 

The  results  produced  by  the  fare  collection  model  can  be  used  by  transit 
systems  for  a number  of  different  purposes.  Examples  of  these  are  the 
following: 

Determination  of  required  number  of  machine  units 
Reliability  and  maintainability  specifications 
Impact  of  changes  in  passenger  demand 
Effect  of  maintenance  policy  changes 
Effect  of  changes  in  fare  collection  method 
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MAINTAINABILITY 


PASSENGER  DELAY  CURVES 
(ARROW  SHOWS  DIRECTION  OF 
DECREASING  AVERAGE  PASSENGER 
DELAY) 


RELIABILITY 


A.  Evaluation  — How  does  the  system  perform  for  a given  reliability 
and  maintainability? 

B-] . Sensitivity  Analysis  — What  happens  to  the  average  delay  when 
maintainability  improves? 

82*  Sensitivity  Analysis  — What  happens  to  the  average  delay  when 
reliability  improves? 

C.  Specification  Determination  — How  high  must  the  reliability  be 
to  obtain  the  desired  level  of  average  delay? 

D-|.  Tradeoff  Analysis  — If  reliability  declines,  by  how  much  must 
maintainability  improve  in  order  to  retain  current  performance? 

D2.  Tradeoff  Analysis  — If  reliability  improves,  by  how  much  can 
maintainability  requirements  be  relaxed  while  still  retaining 
current  performance? 


FIGURE  1.  TYPES  OF  FARE  COLLECTION  DEPENDABILITY  ANALYSES 
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o Determination  of  Required  Number  of  Machine  Units 


An  important  consideration  in  designing  a fare  collection  system  is 
determining  how  many  machine  units  are  required  in  each  service  area,  i.e., 
how  many  coin  changers,  ticket  vendors,  gates,  etc.,  are  needed.  This  is  an 
important  consideration,  since  these  machines  are  expensive  (note  — a cost 
module  is  described  in  Chapter  6 which  can  be  used  to  estimate  these 
expenses).  While  an  insufficient  number  will  adversely  affect  station 
performance,  an  excessive  number  will  cause  a large  unnecessary  expense, 
particularly  if  this  happens  for  a large  number  of  stations. 

By  conducting  a sensitivity  analysis  between  passenger  delay  and  the  number 
of  machine  units  in  place  at  a particular  station,  one  can  choose  the  minimum 
number  of  units  such  that  the  anticipated  passenger  delay  remains  within 
acceptable  limits. 

o Reliability  and  Maintainability  Specifications 

The  level  of  specification  for  the  reliability  and  maintainability  of  the 
fare  collection  system  units  is  very  important.  If  the  specifications  are 
too  low,  the  resulting  reliability  and  maintainability  will  be  too  low, 
producing  unacceptable  passenger  delays  and,  hence,  expensive  retrofits  or 
even  more  expensive  reprocurements.  However,  if  the  specifications  are  too 
high,  the  units  will  become  unnecessarily  expensive  to  obtain,  or  even 
infeasible  to  produce  at  all,  leading  to  unreasonable  expectations  followed 
by  failure  to  meet  specifications,  cost  overruns,  and  possible  lawsuits. 
Either  overspecification  or  underspecification  thus  gets  the  procurement 
process  off  to  a poor  start. 

With  the  model,  one  can  evaluate  a fare  collection  system  under  a set  of 
proposed  reliability  and  maintainability  specifications.  If  the  resulting 
congestion  and  delay  is  unacceptable,  one  can  then,  through  sensitivity  and 
tradeoff  analyses,  increase  one  or  both  specifications  until  an  acceptable 
performance  is  reached.  If,  on  the  other  hand,  the  congestion  and  delay  is 
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acceptable,  one  can  then,  again  by  means  of  sensitivity  and  tradeoff 
analyses,  decrease  the  specifications,  as  long  as  the  performance  remains 
acceptable,  in  order  to  decrease  system  costs. 

o Impact  of  Changes  in  Passenger  Demand 

Another  use  of  the  model  results,  particularly  for  new  systems,  is  to 
determine  the  accuracy  required  in  predicting  passenger  demand.  If  a 
significant  increase  in  passenger  demand  does  not  cause  a large  increase  in 
passenger  delay,  then  the  accuracy  of  the  demand  estimate  is  not  that 
critical.  However,  should  a small  increase  in  demand  cause  a large  rise  in 
passenger  delay,  then  the  demand  estimate  must  be  much  more  accurate. 

Knowing  this  sensitivity  of  system  delay  to  passenger  demand  thus  allows  the 
demand 

prediction  to  be  carried  out  at  an  appropriate  level  of  accuracy. 

An  additional  use  of  the  sensitivity  of  the  fare  collection  system  to 
increase  (or  decrease)  demand  is  to  find  out  in  advance  what  changes  need  to 
be  made  in  the  system  in  the  future  due  to  changes  in  demand.  These  changes 
could  arise  fromincreased  use  of  the  system  or  from  significant  changes 
either  in  the  transit  system  itself  or  its  feeder  systems  (buses  or  parking 
lots) . 

o Effect  of  Maintenance  Policy  Changes 

Maintenance  of  fare  collection  equipment  can  be  carried  out  in  a number  of 
ways.  One  can  have  an  attendant  present  at  every  station  (or  only  at 
heavily-used  ones)  to  correct  minor  failures  and  to  call  in  immediate  aid  for 
more  significant  ones.  One  can  make  repairs  immediately  as  failures  occur, 
or  let  failures  accumulate  until  enough  occur  at  a given  location  to  justify 
sending  a repairman.  One  can  forego  preventive  maintenance,  or  have  such 
maintenance,  with  either  shorter  or  longer  periods  between  successive 
preventive  maintenance  actions.  Since  maintenance  can  significantly  affect 
system  performance,  it  is  important  to  be  able  to  verify  the  impact  of 
changing  the  manner  in  which  it  is  conducted.  By  running  the  model  under 
various  maintenance  schemes,  one  can  compare  them  in  terms  of  resulting 
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system  performance,  thus  allowing,  for  example,  a comparison  of  the  minimum 
reliability  or  number  of  machine  units  needed  for  adequate  performance  under 
different  candidate  maintenance  policies. 

o Effect  of  Change  in  Fare  Collection  Method 

The  model  may  be  used  to  evaluate  changes  in  fare  collection  media,  i.e., 
ticket-only,  token-only,  cash-only,  cash-or-ticket , cash-or-ticket-or-pass, 
etc.  The  choice  of  media  dictate  the  type  of  arrangement  of  machinery  used 
in  the  system.  For  example,  one  can  compare  a system  in  which  the  gates  will 
only  accept  tickets,  thus  simplifying  the  gate  design  in  return  for  requiring 
heavier  use  of  the  ticket  vendors  and  coin  changers,  with  a system  in  which 
the  gates  will  accept  cash,  thus  increasing  the  cost  and  possibly  decreasing 
the  reliability  of  the  gates,  but  reducing  the  use  of  the  ticket  vendors  and 
coin  changers.  In  this  way,  different  fare  collection  methods  can  be 
evaluated  before  expensive  machinery  is  ordered  or  installed. 

The  model  can  evaluate  systems  with  either  entry  control  only  (i.e.,  no  exit 
processing  of  passengers  necessary),  both  entry  and  exit  control  (i.e., 
passengers  must  use  a ticket  or  other  medium  to  enter  the  system  at  their 
origin  stations,  as  well  as  to  leave  the  transit  system  at  their  destination 
stations),  or  barrier-free  control  (i.e.,  passengers  may  enter  and  exit  the 
system  without  impedance  by  gates,  but  must  obtain  and  validate  a ticket,  and 
must  produce  such  a validated  ticket  on  demand). 


11 


3.  TECHNICAL  APPROACH 


The  basic  approach  to  the  model  is  to  investigate  the  operation  of  the  fare 
collection  system  as  a multiple-server  queue,  with  passengers  as  customers, 
machine  units  as  servers,  and  a first-come-first-served  service  discipline. 

In  addition  to  the  normal  queue  features,  there  is  the  additional  aspect  that 
the  number  of  servers  (machine  units)  itself  decreases  and  increases  as  the 
units  fail  and  are  repaired,  respectively. 

Two  types  of  simulation  models  have  been  developed;  a basic  model  which 
covers  a single  service  area  (such  as  a gate  or  ticket  vendor),  and  an 
extended  model  which  covers  several  service  areas  in  tandem.  In  addition,  an 
analytical  model  has  been  developed  which  covers  the  single  service  area 
without  the  use  of  simulation. 

Note:  The  random  variables  in  the  models  as  presented  are  assumed  to  be 

exponentially  distributed.  This  is  done  for  tractability  and  to  provide  a 
reasonable  generic  probability  distribution  if  detailed  data  for  such  a 
distribution  is  not  available.  The  simulation  models  can  be  routinely 
changed  to  different  distributions,  however,  if  such  data  is  in  fact 
available. 

3.1  The  Single-Service-Area  Simulation  Model 

Figure  2 shows  a situation  in  which  there  is  a single  service  area.  Several 
machine  units  make  up  the  service  area,  each  of  which  can  serve  one  passenger 
provided  they  are  not  otherwise  busy  or  out  of  service  because  of  failure. 

The  model  is  an  event-oriented  simulation  in  which  the  next  event  to  be 
processed  is  the  earliest-occurring  of  the  four  basic  events:  passenger 

arrivals,  passenger  departures,  equipment  breakdown,  and  equipment  repair. 
Events  are  processed  until  the  time  of  the  prospective  next  event  is  no 
longer  within  the  time  period  being  simulated.  Each  of  the  four  events  is 
described  in  detail  in  Sections  3.1.1  ~ 3.1.4. 
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3.1.1  Passenger  Arrivals 


The  arrival  of  passengers  in  the  system  is  assumed  to  be  a sequence  of 
independent  events,  in  which  the  time  between  successive  arrivals  is  an 
exponentially-distributed  random  variable.  (While  the  fact  of  an 
exponential  distribution  is  not  essential  to  the  model,  the  independence  of 
successive  arrivals  is.)  Passengers  are  assumed  to  arrive  in  groups.  For 
passengers  exiting  the  station  after  arriving  on  a train  or  for  those 
entering  the  station  after  arriving  on  a feeder  bus,  the  group  size  can  be 
large.  However,  for  passengers  entering  the  station  after  arriving  on  foot, 
or  by  taxi,  or  by  private  car,  the  group  size  can  be  small,  possibly 
consisting  of  only  one  passenger. 

Upon  arrival,  a passenger  immediately  begins  service  if  a machine  unit  is 
available.  If  so,  a departure  time  is  calculated  and  put  into  the  departure 
stack,  which  is  an  array  of  departure  times  of  passengers  in  service,  sorted 
in  ascending  order,  that  is  used  to  determine  the  time  of  the  next  departure. 
If  no  machine  is  available,  however,  the  passenger  enters  a queue,  which  is 
maintained  on  a first-in-first-out  basis. 

Passengers  who  find  an  available  unit  upon  arrival  and  so  immediately  enter 
service  have  delay  times  of  zero.  Passengers  entering  service  from  the  queue 
have  their  (non-zero)  delay  times  calculated  at  the  time  they  enter  service. 
The  delay  time  is  the  interval  between  the  arrival  time  and  the  time  of  entry 
into  service. 

3.1.2  Passenger  Departures 

The  service  (or  processing)  time  of  passengers  is  also  assumed  to  be  an 
independent  random  variable  with  an  exponential  distribution.  The  departure 
time,  calculated  upon  entry  into  service,  is  equal  to  the  time  of  arrival 
plus  the  service  time.  Upon  departure,  the  record  of  the  departure  time  is 
removed  from  the  departure  stack. 
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3.1.3  Equipment  Breakdown 


The  reliability  of  the  equipment  being  modeled  is  given  in  terms  of  the  MCBF, 
or  mean  cycles  between  failures.  The  reliability  is  given  in  terms  of  cycles 
rather  than  time  because  the  equipment  (especially  the  mechanical  parts)  is 
exposed  to  failure  by  actual  use  rather  than  simply  by  the  passage  of  time. 
Upon  the  departure  of  a passenger  from  service,  a random  draw  is  made  with 
probability  1/MCBF  that  the  machine  unit  just  used  breaks  down.  If  it  does 
not  break  down,  it  becomes  available  for  service  and,  if  a queue  exists,  the 
first  passenger  in  the  queue  enters  service.  If  it  does  break  down,  however, 
a repair  time  for  it  is  calculated  (see  "Equipment  Repair"  below)  and  it  is 
placed  in  the  repair  stack, which  is  an  array  of  return-to-service  times  of 
failed  units,  arranged  in  ascending  order.  The  number  of  machine  units 
available  is  decreased  by  1 if  a breakdown  occurs,  and  no  passengers  are 
brought  in  from  the  queue. 

3.1.4  Equipment  Repair 

The  time  necessary  to  repair  a failed  unit  is  also  assumed  to  be  an 
independent  random  variable  with  an  exponential  distribution.  Note  that 
because  a failed  unit  impedes  system  operation  from  the  time  it  fails  until 
the  time  it  returns  to  service,  "repair  time"  includes  not  only  the  time  to 
actually  do  the  repair,  but  also  the  administrative  downtime,  i.e.,  the  time 
necessary  to  detect  a failure,  have  a repair  crew  get  to  the  unit,  and 
diagnose  the  failure.  The  return-to-service  time  is  calculated  upon 
breakdown  as  the  sum  of  the  time  of  breakdown  and  the  total  repair  time,  and 
entered  into  the  repair  stack.  Upon  repair,  the  unit  is  removed  from  the 
repair  stack,  and  the  number  of  units  available  is  increased  by  1 . If  a 
queue  exists  when  a unit  returns  to  service,  the  first  passenger  in  the  queue 
enters  service. 

3.2  The  Multiple-Service-Area  Simulation  Model 

An  actual  fare  collection  system  at  a station  will  have  not  one,  but  several 
service  areas  which  passengers  may  use  during  the  fare  collection  process. 

For  example,  a station  may  have  coin  changers,  ticket  vendors,  and  gates.  An 
arriving  passenger  may:  1)  have  the  proper  ticket  and  proceed  directly  to 


15 


the  gate;  2)  may  need  to  purchase  a ticket  and  therefore  proceed  first  to 
the  ticket  vendors  and  then  to  the  gate  area;  or  3)  may  need  the  proper 
change  to  buy  a ticket  and  so  proceed  through  all  three  service  areas.  The 
multiple-service-area  model  (Figure  3)  investigates  this  type  of  situation. 
This  model  contains  three  service  areas  (a  change  to  a different  number  of 
areas  is  easy  to  implement):  service  area  #1  (the  coin  changers,  in  the 

above  example),  service  area  #2  (the  ticket  vendors),  and  service  area  #3 
(the  gates).  The  probability  is  p-|  that  an  arriving  passenger  begins  service 
at  area  #1,  P2  that  he  begins  service  at  area  #2,  and  P3  that  he  begins 
service  at  area  #3.  (The  probabilities  p-|,  P2,  and  P3  are  specified 
parameters  whose  sum  is  1.)  After  completing  service  at  his  original  service 
area,  the  passenger  continues  to  the  next  service  area  downstream  (i.e.,  from 
area  #1  to  area  #2,  and  from  area  #2  to  area  #3),  until  he  completes  service 
at  area  #3  (the  gate),  upon  which  he  departs  the  system.  Congestion  (queue 
length)  and  delay  times  are  tabulated  at  each  service  area,  and  the  delay 
time  for  a passenger  is  the  sum  of  the  delay  times  at  each  service  area  the 
passenger  uses. 

A departure  from  service  area  #1  creates  a simultaneous  arrival  at  service 
area  #2.  Similarly,  a departure  from  service  area  #2  creates  a simultaneous 
arrival  at  service  area  #3.  A departure  from  service  area  #3  (the  gates)  is 
a pure  departure,  as  the  passenger  then  leaves  the  system. 

3.3  The  Analytical  Model 

A simulation  model  is  a reasonably  straightforward  way  to  investigate  a fare 
collection  system.  However,  simulation  does  have  drawbacks.  The  underlying 
probability  distributions  for  congestion  and  delay  are  not  specifically 
determined  by  a simulation.  Rather,  random  draws  on  these  distributions  are 
taken  for  each  passenger  simulated  and,  after  some  number  of  passengers  have 
been  processed,  the  frequency  distributions  of  congestion  and  delay  of  these 
passengers  are  taken  as  an  approximation  to  the  actual  distributions. 

However,  these  are  only  approximations,  and  are  subject  to  a number  of 
statistical  sensitivities.  Furthermore,  because  each  passenger  must  be 
generated  and  processed  individually,  a simulation  of  this  type  will  require 
a substantial  amount  of  computer  time  if  a large  number  of  passengers  need  to 
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FIGURE  3.  PASSENGER  FLOW  THROUGH  A FARE  COLLECTION  SYSTEM  WITH  THREE  SERVICE  AREAS 
(COIN  CHANGERS,  TICKET  VENDORS,  AND  GATES) 


be  processed.  This  would  be  necessary  not  only  for  stations  with  a high 
passenger  demand,  but  also  for  stations  which  need  to  be  simulated  for  a long 
period  of  time  due  to  the  statistical  sensitivities  mentioned  above. 

One  example  of  the'  latter  is  a station  where  the  flow  of  incoming  passengers 
is  at  or  near  processing  capacity.  In  such  a station,  large  queues  will 
develop  when,  due  to  random  fluctuations,  a surge  of  passenger  arrivals 
occur.  Furthermore,  these  large  queues  will  take  a long  time  to  dissipate, 
since  the  station  has  very  little,  if  any,  extra  processing  capacity.  These 
surges  thus  cause  considerable  congestion  and  delay;  in  fact  they  cause  most 
of  the  total  congestion  and  delay  in  this  case.  Because  these  random 
passenger  arrival  surges  do  not  occur  very  often,  and  because  the  model  needs 
a reasonably  large  sample  of  surges  to  provide  accurate  statistical  results, 
the  simulation  must  be  run  for  a long  period  of  time. 

Another  example  is  a station  with  machine  units  which  have  low  failure  rates, 
but  long  repair  times.  Failures  therefore  do  not  happen  often,  but 
significantly  affect  the  system  when  they  do.  Because,  similar  to  the  above 
case,  the  model  needs  to  have  a large  enough  sample  of  these  failures  in 
order  to  provide  statistical  results  of  sufficient  accuracy,  we  need  to  run 
the  simulation  for  a long  period  of  time  in  this  case  as  well. 

We  therefore  have  developed,  in  addition  to  the  simulation  model,  an 
analytical  model  to  analyze  a service  area  for  a fare  collection  system.  The 
analytical  model  directly  formulates  and  solves  the  mathematical  equations 
for  the  underlying  queue  length  probabilities,  rather  than  obtaining  an 
approximate  solution  by  simulation.  As  mentioned  in  Section  1,  while  the 
model  only  considers  a single  service  area  (with  no  grouping  of  passengers), 
rather  than  the  entire  fare  collection  system,  and  requires  certain 
assumptions  on  the  arrival,  processing,  and  repair  rate  distributions  (i.e., 
that  they  are  Independent  and  exponentially  distributed) , it  does  obtain  the 
congestion  and  delay  distributions  directly.  Hence,  the  computer  time 
required  by  the  model  does  not  depend  on  the  passenger  arrival  rates,  degree 
of  saturation  (passenger  arrival  rate  compared  to  service  capacity),  or 
failure  and  repair  rates.  The  analytical  model  can  therefore  be  used  as  the 
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initial  investigating  tool  for  large  problems,  allowing  many  runs  to  be  made 
with  many  different  parameter  values,  without  using  an  inordinate  amount  of 
computer  time.  From  these  runs,  a far  smaller  number  of  scenarios  of 
interest  can  be  identified  and  investigated  with  the  simulation  model. 

The  approach  to  such  a model  is  a modification  of  that  derived  by  Neuts  and 
Lucantoni  (9)  for  the  multiple-server  exponential  queue  with  a randomly- 
varying  number  of  servers.  As  with  the  simulation  model,  the  service  area  is 
represented  as  a queue  with  as  many  servers  as  there  are  machine  units. 
Passengers  are  assumed  to  arrive  according  to  a Poisson  process  (i.e.,  the 
interarrival  time  between  any  two  successive  passengers  is  independent  of  any 
other  passenger  interarrival  time  and  is  exponentially  distributed),  and 
their  processing  times  are  also  assumed  to  be  independent  and  exponentially 
distributed.  Similarly,  times  to  failure  and  repair  times  for  the  machine 
units  are  independent  and  exponentially  distributed. 

The  modification  to  the  Neuts-Lucantoni  model  is  that  the  rate  at  which 
failures  occur  in  a service  area  is  proportional  not  to  the  number  of 
functioning  machine  units,  as  specified  in  (9),  but  rather  the  minimum  of  the 
number  of  functioning  machine  units  and  the  queue  length.  In  other  words, 
machine  units  cannot  fail  while  they  are  idle.  The  modification  results  in  a 
straightforward  change  in  the  arguments  in  (9),  which  are  then  used  to 
directly  calculate  the  probabilities  (xij),  where 

Xij  = Probability  that  the  queue  length  (number  of  customers  either  being 
served  or  waiting  for  service)  is  i and  that  j machines  are 
functioning. 
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The  queue  length  probabilities  (Ij^)  are: 


Ij^  = Probability  that  there  are  i customers  either  being  served  or 
waiting  for  service 


and  the  mean  queue  length  (congestion)  is 

00 

L = ^ ili  (2) 

i=0 

The  mean  time  in  the  system  (i.e.,  the  time  spent  both  in  service  and  waiting 
for  service)  is  given  by  Little's  Formula  (10,  P.  60)  as 

W = L/A  , (3) 

and  the  mean  delay  time  as 

Wd  = (L/A)  - (1/y)  = W - (1/y)  , (4) 

where  A is  the  arrival  rate  and  y the  processing  rate  of  customers. 

By  the  above  process  the  analytical  model  computes  the  probability 
distribution  of  the  queue  length  (congestion),  and  the  mean  congestion  and 
delay. 
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4.  USING  THE  MODEL 


4.1  Data  Requirements 

The  data  requirements  of  the  model  are  of  two  kinds:  hardware  and  passenger 

flow.  The  hardware  data  include  reliability  and  maintainability  data,  as 
well  as  the  passenger  processing  rate  per  machine  unit  and  the  number  of 
units  provided,  for  each  of  the  service  areas  in  the  station.  The  passenger 
flow  data  include  the  passenger  arrival  rate,  group  size,  and  division  of 
passenger  flow  to  the  various  service  areas.  While  average  values  of  these 
data  are  acceptable  (one  can  assume,  as  the  present  models  do,  an  exponential 
distribution  with  the  given  average  value  as  the  mean),  for  simulation  models 
actual  frequency  distribution  information  is  preferable  if  available,  so  that 
the  variability  of  the  data  from  the  average  can  be  taken  into  account. 

Specifically,  the  data  requirements  are  as  follows: 

Passenger  arrival  rate  (one  parameter  for  entire  system)  — The  hourly  rate 
at  which  passengers  arrive  at  the  fare  collection  system  during  the  peak 
period.  This  may  be  established  by  a priori  demand  estimation  techniques 
or  by  an  actual  passenger  count.  Assuming  a stable  (l.e.,  steady-state) 
mean  flow  during  the  peak  period,  a good  estimate  of  the  actual  frequency 
distribution  of  flow  can  be  obtained  by  taking  successive  passenger 
counts  over  many  brief  periods  of  time  (such  as  five  minutes). 

Group  size  (one  parameter  for  the  entire  system)  — The  size  of  a group  of 
arriving  passengers  (this  can  be  given  either  as  a single  deterministic 
value  or  as  a random  variable  defined  by  a probability  distribution). 
Passengers  exiting  from  the  transit  system  and  those  arriving  in  large 
vehicles  such  as  commuter  vans  or  buses  will  generally  come  or  go  in 
large  groups,  while  passengers  arriving  on  foot  or  in  small  vehicles 
such  as  cars  will  usually  arrive  in  small  groups  (frequently  consisting 
of  just  one  passenger).  This  information  is  not  easy  to  obtain,  as  it 
requires  direct  observation  of  passengers,  and  careful  observation  at 
that,  since  groups  may  form,  merge,  divide,  or  dissipate  rapidly. 
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Passenger  processing  rate  (one  parameter  for  each  service  area)  — The 
hourly  rate  at  which  a machine  unit  in  the  service  area  can  process 
passengers,  and  hence  the  unit's  capacity  to  handle  passenger  flow.  This 
is  a difficult  value  to  obtain  under  actual  operating  conditions.  While 
design  capacit.ies  are  normally  given  for  a machine  unit,  this  merely 
represents  the  theoretical  maximum  rate  at  which  passengers  can  be 
processed.  The  actual  rate  in  the  field,  even  assuming  the  unit  is  never 
idle,  will  be  significantly  less,  because  of  passenger  time  lost 
accessing,  fumbling  for,  or  dropping  coins  or  tickets,  reading  and 
comprehending  (or  misunderstanding)  instructions,  making  decisions  while 
using  the  unit,  making  errors  and  correcting  them,  and  simply  using  the 
unit  at  a slower  rate  than  its  capacity.  The  actual  processing  rate  is 
thus  difficult  to  specify  and  estimate  beforehand,  and  so  can  only  be 
obtained  by  an  actual  count  which  simultaneously  obtains  for  a given  time 
period  the  number  of  passengers  processed  and  the  amount  of  time  the  unit 
was  actually  busy.  The  latter  data  are  not  routinely  collected,  so  that 
a special  effort  must  usually  be  made  to  obtain  the  passenger  processing 
rates.  Furthermore,  in  a new  or  significantly  upgraded  fare  collection 
system  such  data  by  definition  does  not  exist,  and  so  can  only  be 
estimated  beforehand,  as  difficult  and  imprecise  as  that  process  is  to 
accomplish. 

Failure  rates,  or  reliability  (one  parameter  for  each  service  area) — The 
rate  at  which  failures  occur  to  a machine  unit  which  make  it  unable  to 
process  passengers.  Because  of  the  mechanical  nature  of  the  equipment, 
and  the  resulting  fact  that  the  basic  measure  of  stress  to  the  unit,  and 
hence  exposure  to  failure,  is  the  use  of  the  unit  by  an  individual 
passenger,  i.e.,  one  cycle,  the  measure  of  failure  rate  is  best  given  in 
terms  of  Mean  Cycles  Between  Failures  (MCBF) . If  necessary,  a Mean  Time 
Between  Failures  (MTBF)  measure  can  be  used  (where  the  time  used  is  hours 
in  actual  use,  rather  than  elapsed  clock  time),  from  which  MCBF  can  be 
obtained  by  dividing  MTBF  by  the  mean  passenger  processing  time. 

For  an  existing  system,  the  MCBF  can  be  obtained  by  direct  observation,  or 
by  combining  passenger  use  records  and  failure  logs.  For  new  or 
significantly  upgraded  systems,  since  such  data  does  not  exist, 
estimates,  predictions,  specifications,  or  manufacturer's  descriptions  of 
reliability  must  be  used. 
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Repair  times,  or  maintainability  (one  parameter  for  each  service  area)  — 
The  elapsed  time,  in  hours,  between  the  failure  of  a machine  unit  and  its 
return  to  service.  Notice  that  there  are  other  concepts  of  repair  times, 
such  as  the  time  it  takes  to  carry  out  the  actual  repair  on  the  unit,  the 
time  between  the  discovery  of  the  failed  unit  and  its  return  to  service, 
etc.  However,  the  repair  time  defined  above  is  the  appropriate  one  for 
our  purposes,  as  it  denotes  the  length  of  time  the  system  will  function 
at  a reduced  capacity  due  to  the  failure. 

For  an  existing  system,  repair  time  can  be  obtained  by  direct 
observation,  or  by  consulting  repair  logs.  If  the  latter,  care  must  be 
taken  to  include  an  estimate  of  the  time  necessary  to  detect  a failure 
which  has  occurred,  as  the  logs  will  at  best  only  give  the  time  from 
detection  to  return-to-service.  (In  addition,  if  the  logs  only  give  the 
times  to  actually  repair  the  units,  estimates  must  also  be  made  of  the 
time  required  for  a repairman  to  travel  to  the  failed  unit,  as  well  as 
any  additional  travel  time  which  may  be  necessary).  For  new  or 
significantly  upgraded  systems,  the  repair  times  must  be  estimated. 

Number  of  machine  units  (one  parameter  for  each  service  area)  — The  number 
of  machine  units  nominally  available  for  passenger  use,  in  the  absence  of 
failures.  In  a small  station,  space  considerations  may  restrict  the 
number  of  machine  units  possible  in  a given  service  area,  even  if  a 
greater  number  would  be  warranted  by  dependability  considerations. 

Division  of  passenger  flow  to  service  areas  (one  parameter  for  each  service 
area)  — The  proportion  of  arriving  passengers  who  begin  their  use  of  the 
fare  collection  system  in  that  particular  service  area  (see  Figure  3). 

In  an  aggregate  sense,  passenger  flow  division  can  be  obtained  by 
determining  flow  rate  counts  for  each  of  the  service  areas  and  comparing 
them  to  the  overall  passenger  arrival  rate  (or  the  total  passenger  flow 
through  the  gates,  which  should  be  equivalent  unless  the  system  is 
overloaded) . 
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4.2  Suggested  Precautions 


Any  model,  including  this  one,  is  not  an  exact  replication  of  reality,  but 
rather  an  approximation.  As  such,  it  exhibits  differences  from  reality  which 
can  mislead  an  unwary  user.  The  following  are  some  precautions  which  a user 
should  take  into  account. 

A Monte  Carlo  simulation  model  depends  on  a random  number  generator  to 
produce  the  random  variation  in  input  data  that  is  so  important  in  putting 
the  simulated  system  "through  its  paces."  Virtually  all  computer  systems 
have  such  a generator,  which  produces  numbers  which  vary  randomly  between  0 
and  1.  However,  some  of  these  generators  can  be  defective.  In  the  early 
stages  of  this  effort,  inexplicable  discrepancies  occurred  between  expected 
and  actual  results.  It  turned  out  that  the  random  number  generator  being 
used  by  the  computer  tended  to  generate  high  numbers  (values  near  1) 
significantly  more  often  than  low  numbers  (values  near  0).  A different 
random  number  generator  was  substituted,  and  the  discrepancies  disappeared. 
This  is  no  isolated  occurrence,  since  such  problems  have  been  described  in 
the  literature.  One  must  therefore  make  sure  that  the  random  number 
generator  being  used  is  an  accurate  one. 

The  collection  of  data  may  also  introduce  some  inaccuracies.  For  example,  to 
obtain  the  data  for  passenger  processing  rates,  one  must  simultaneously  count 
the  number  of  passengers  using  the  machine  units  and  measure  the  total  time 
that  these  unit's  are  busy  (the  processing  rate  is  the  quotient  of  these  two 
values).  If  arrivals  occur  frequently  but  intermittently,  the  busy  periods 
of  the  machine  units  will  be  numerous  but  brief.  It  is  often  difficult  to 
tell  exactly  when  a busy  period  begins  or  ends,  and  errors  of  several  seconds 
may  be  Introduced,  either  because  of  this  or  merely  through  observer 
inattention  (this  particular  data  collection  effort  requires  much 
concentration  by  the  observer).  Because  of  this  difficultly,  the  data 
collection  process  may  introduce  significant  inaccuracies  into  the  data  when 
the  busy  periods  of  the  machine  units  are  numerous  but  brief. 
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The  actual  probability  distribution  of  the  data  (i.e.,  whether  it  is  normal, 
exponential,  or  some  other  shape)  may  be  quite  important  to  system 
performance,  as  different  distributions  have  different  probabilities  of 
generating  extreme  values  (i.e.,  values  which  are  extremely  different  from 
the  mean,  and  especially  values  which  are  extremely  larger  than  the  mean). 

It  is  often  these  extreme  values  (such  as  a random  overwhelming  surge  of 
passenger  arrivals,  or  a random  rash  of  failures)  which  cause  system  delay. 
Nevertheless,  quite  often  only  the  mean  value  is  available  for  the  desired 
data,  and  so  some  particular  type  of  distribution  (such  as  the  exponential 
one)  must  be  assumed.  If  this  is  so,  a possibly  significant  source  of  delay 
may  be  missed. 

Also,  quite  often  only  estimates  are  available  for  much  of  the  required  data. 
While  this  is  especially  so  for  new  or  significantly  upgraded  systems,  which 
practically  by  definition  have  not  had  opportunity  to  accumulate  actual 
performance  records,  it  also  holds  for  some  data,  such  as  passenger 
processing  rates  for  machine  units,  in  existing  systems  as  well.  Estimates, 
however,  may  be  inaccurate  to  a greater  or  lesser  degree,  which  would  then 
affect  the  results  of  the  analysis. 

The  final  precaution  applies  to  the  simulation  models.  Clearly,  breakdowns 
and  repairs  of  machine  units  occur  far  less  frequently  than  do  passenger 
arrivals  and  departures.  If  the  simulation  is  run  for  only  a short  simulated 
time  period,  only  a few  breakdowns  and  repairs  will  have  occurred.  This  will 
make  the  results  quite  volatile  and  unreliable,  since  they  will  be  highly 
sensitive  to  the  addition  or  subtraction  of  only  a few  breakdowns,  or  even  a 
change  in  the  time  of  occurrence  (i.e.,  early  or  late  in  the  simulated  time 
period)  of  the  breakdowns  which  do  occur.  These  kinds  of  events  can  take 
place  quite  readily,  simply  due  to  the  random  fluctuations  provided 
(intentionally)  by  the  random  number  generator.  Furthermore,  if  the 
simulation  is  run  for  a short  simulated  time  period,  the  starting  and  ending 
effects  can  be  significant  (starting  effects  are  what  happens  between  the 
start  of  the  simulation  and  the  time  it  "settles  down"  to  a steady-state 
situation,  while  ending  effects  take  into  account  the  delay  of  passengers 


25 


still  in  the  system  when  the  simulation  ends).  Therefore,  the  simulation 
models  must  be  run  for  a long  enough  time  (or  enough  separate  runs  must  be 
made)  for  a sufficient  amount  of  breakdowns  and  repairs  to  have  occurred,  and 
for  starting  and  ending  effects  to  be  no  longer  significant.  This  may 
require  a large  amount  of  computer  time. 

With  these  precautions  in  mind,  the  models  can  be  used  by  transit  systems  as 
a useful  tool  in  determining  fare  collection  system  requirements. 
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5 . SAMPLE  RUNS 


To  demonstrate  the  model  and  its  application  to  fare  collection  system 
analysis,  a number  of  sample  runs  were  performed.  For  the  single-service- 
area  simulation  model,  hypothetical  test  data  were  generated  to  describe  a 
ticket-vending  facility  under  various  circumstances.  In  addition,  to 
demonstrate  the  use  of  the  multiple-service-area  model  for  actual  fare 
collection  system  analysis,  sample  runs  were  also  carried  out  based  on 
preliminary  data  obtained  from  the  Miami  Dade  County  Transit  Authority 
describing  their  stations  and  fare  collection  system  (the  sample  runs, 
described  below,  thus  use  preliminary  configurations  and  are  not  intended  to 
reflect  the  final  configuration  of  the  actual  Miami  system).  A base  run  was 
made  using  a representative  station,  followed  by  various  sensitivity  runs  to 
assess  the  impact  of  possible  changes  in  reliability,  maintainability,  and 
the  number  of  machine  units  at  the  station.  Additionally,  a single  service 
area  of  the  representative  station  was  used  to  provide  a test  example  for  the 
analytical  model.  In  a similar  manner  as  with  the  simulation  model  example, 
a base  run  was  made,  followed  by  various  sensitivity  and  tradeoff  analyses. 

5.1  Ticket-Vendor  Simulation  (Single-Service-Area  Model) 

In  this  analysis  it  is  assumed  that  a ticket  vendor  has  low  reliability,  and 
that  significant  delays  are  occurring  as  a result.  The  aim  is  to  assess  the 
base  case,  then  examine  a number  of  other  cases  to  see  how  delay  is  affected 
by  various  alternative  approaches.  These  cases  are  defined  as  follows  (see 
Table  1 for  the  data): 

1.  The  base  case,  in  which  there  is  a problem  with  low 
reliability  (note  the  MCBF  of  120). 

2.  An  extra  ticket  vendor  is  added. 

3.  The  reliability  of  the  ticket  vendors  is  improved. 

4.  The  maintainability  of  the  ticket  vendors  is  improved. 

5.  The  reliability  is  improved  as  a tradeoff  against  worse 
maintainability . 

6.  The  maintainability  is  improved  as  a tradeoff  against  worse 
reliability. 
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7.  Passengers  arrived  in  much  larger  group  clusters  than  in  the  base 
case . 

8.  Case  7,  but  an  extra  ticket  vendor  is  added. 

9.  There  are  twice  as  many  ticket  vendors,  but  also  twice  as  many 
passengers. 

10.  An  unusually  large  passenger  flow  occurs  initially,  then  recedes. 
To  handle  this,  an  extra  ticket  vendor  is  added. 


TABLE  1.  DATA  FOR  SINGLE-SERVICE-AREA  ANALYSES 


Group 

Passenger 

Ser- 

Group 

Arrival 

Arrival 

vice 

Maintain- 

Size 

Rate 

Rate 

Rate 

Reli- 

ability 

(per 

(per 

(per 

Number 

(per 

ability 

(MTTR; 

Case 

hr . ) 

hr . ) 

hr . ) 

of  Units 

hr.) 

(MCBF) 

seconds) 

1 . 

2 

250 

500 

3 

300 

120 

720 

2. 

2 

250 

500 

4 

300 

120 

720 

3. 

2 

250 

500 

3 

300 

480 

720 

4. 

2 

250 

500 

3 

300 

120 

180 

5. 

2 

250 

500 

3 

300 

480 

2880 

6. 

2 

250 

500 

3 

300 

30 

180 

7. 

20 

25 

500 

3 

300 

120 

720 

8. 

20 

25 

500 

4 

300 

120 

720 

9. 

2 

500 

1000 

6 

300 

120 

720 

10. 

2 

800/250* 

1600/500* 

4 

300 

120 

720 

*Arrival  rates  are  at 

first  figure 

for  first 

360  seconds  of  simulated  time. 

then 

at  second  figure 

for  the  remaining  time 

• 

Underlining 

indicates 

changes  from 

the  base 

case  (Case 

1) 

Cases  2,  3,  and  4 are  sensitivity  analyses;  Cases  5 and  6 tradeoff  analyses; 
and  Cases  7,  8,  9,  and  10  are  variants  of  the  base  case.  The  results  of  the 
model  runs  for  the  various  cases  are  shown  in  Table  2. 
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TABLE  2.  RESULTS  OF  SINGLE-SERVICE-AREA  ANALYSIS 


Case 

1 

2 3 

4 

5 

6 7 

8 9 10 

Mean 

Queue  Length 

15.2 

4.7  3.8 

4.0 

28.0 

6.1  112.3 

23.0  5.6  9.4 

Mean 

Passenger  Delay 
(Seconds) 

65.5 

17.2  11.9 

13.6 

185.8 

34.4  348.5 

77.9  9.0  31.3 

In  the  base  case 

(Case 

1 ) , the  mean 

delay  is 

about 

one  minute  per 

passenger,  which 

is  a large  amount  of  delay,  and  indicates  that  there  is  indeed  a problem  at  this 
service  area.  The  three  cases  in  which  an  improvement  is  made  in  the  physical 
system  (Cases  2,  3,  and  4),  all  improve  the  situation  equally  well,  so  that  in  a 
choice  among  them,  the  least  expensive  alternative  should  be  used.  The  two  tradeoff 
cases  and  the  base  case  (Cases  1,  5,  and  6)  demonstrate  that,  in  this  particular 
system,  if  any  tradeoff  is  made,  it  should  be  in  favor  of  maintainability.  Case  7 
demonstrates  the  adverse  sensitivity  of  the  system  to  heavily  bunched  arrivals, 
which  is  significantly  improved  by  the  addition  of  the  extra  unit  in  Case  8.  Case  9 
indicates  that  the  improvement  in  adding  extra  units  is  more  than  linear,  i.e., 
twice  as  many  units  can  handle  more  than  twice  as  many  passengers.  Finally,  Case  10 
indicates  that  if  the  system  has  an  extra  (fourth)  unit,  it  is  moderately  capable  of 
handling  surges  in  passenger  demand. 

5.2  Full-System  Performance  Simulation  (Multiple-Area  Model) 

The  Miami  Dade  County  Transit  Authority  is  in  the  process  of  designing,  procuring, 
and  constructing  a rapid  transit  system  for  the  Greater  Miami  area.  Specifications 
have  been  developed  for  the  various  fare  collection  facilities  and  contracts  awarded 
to  supply  these  facilities.  Estimates  have  also  been  obtained  on  passenger  flow  and 
machine  capacity.  We  shall  use  these  estimates,  along  with  the  reliability  and 
maintainability  specifications,  as  input  data  for  a sample  model  run  based  on  data 
derived  from  an  actual  system  (2,11). 

The  fare  on  the  Miami  system  is  collected  at  entry  gates  which  accept  either  coins, 
passes,  or  (at  some  stations)  transfers.  Some  gates  are  set  for  full-fare 
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passengers,  while  others  are  set  for  reduced-fare  passengers.  Passes  are  sold 
off-site,  while  transfer  dispensers  and  parking-lot-receipt  machines  operate  in 
the  station  area. 

The  station  analyzed  in  the  model  run  is  derived  from  the  Dadeland  North  station 
during  peak-hour  operation  (this  station  was  selected  because  it  is  a relatively 
important  station  which  has  enough  passenger  demand  to  result  in  significant 
congestion  and  delay  if  enough  machine  units  fail).  For  simplicity,  the  sample 
station  uses  gates  which  accept  (no  transfers),  processes  full-fare  only,  and 
includes  only  the  gates  and  change  machines.  No  exit  processing  of  passengers 
is  required,  so  only  the  entering  passenger  flow  will  be  examined. 

In  order  to  provide  a multiple-service-area  example,  the  station  design  will  be 
based  on  the  preliminary  design  configuration  for  the  Miami  system,  which 
provides  for  change  machines  to  change  bills  and  coins  (the  final  design 
configuration,  which  does  not  include  change  machines,  can  be  analyzed  most 
efficiently  by  using  the  single-service-area  simulation  model  or  the  analytical 
model).  The  station  layout  is  shown  in  Figure  4.  The  estimated  peak-hour 
passenger  flow  at  the  sample  station  is  5400  passengers  per  hour.  Of  these,  30? 
use  the  change  machines,  while  70?  go  directly  to  the  gates.  Passengers  are 
assumed  to  arrive  one  at  a time,  not  in  groups. 

There  are  five  gates  at  the  station.  Each  gate  has  a physical  capacity  to 
process  1800  passengers  per  hour.  As  mentioned  previously,  the  actual  field 
processing  capacity  is  lower  than  this  and  is  difficult  to  obtain.  As  the  Miami 
system  is  not  yet  in  revenue  service,  there  is  no  actual  field  data  available, 
nor  are  there  detailed  analyses  of  the  field  passenger-processing  capacity  of 
these  types  of  gates.  (There  are  data  collection  efforts  underway  at  ^i^ferent 
transit  systems  to  obtain  such  information.  As  discussed  earlier,  it  is  not 
straightforward  to  measure  the  field  processing  capacity).  A rough  rule-of- 
thumb  for  field  processing  capacity  of  75?  of  the  physical  processing  capacity 
is  assumed  for  this  analysis.  Therefore,  the  gate  processing  rate  used  in  the 
model  runs  is  .75  x 1800,  or  1350  passengers  per  hour. 
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PROCESSING  RATE  360/HR. 
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FIGURE  H.  SCHEMATIC  LAYOUT,  PASSENGER  FLOWS,  AND  PROCESSING  RATES  TRACK 

FOR  EXAMPLE  TRANSIT  STATION 


There  are  also  six  coin  changers.  Each  can  process  bills  at  12  transactions  per 
100  seconds,  and  coins  at  20  transactions  per  100  seconds.  The  transactions 
involve  bills  85^  of  the  time  and  coins  15^  of  the  time.  The  mean  transaction 
rate  is  therefore  12*(0.85)  + 20*(0.15)  = 13.2  per  100  seconds,  or  7.92  per 
minute,  or  approximately  475  per  hour.  Using  the  rule-of-thumb  field  processing 
rate  ratio  of  75%,  we  obtain  an  estimated  field  transaction  rate  of  360 
passengers  per  hour  for  the  changers. 

The  reliability  specifications  are  an  MCBF  of  60000  for  the  gates  and  80000  for 
the  changers.  For  maintainability,  the  specifications  are  that  the  time  to 
carry  out  the  actual  repair  on  a unit  shall  not  exceed  0.6  hours  (36  minutes) 
for  the  gates  and  0.3  hours  (18  minutes)  for  the  changers.  We  shall  assume  that 
the  time  to  detect  a failure  and  get  a repairman  to  the  failed  unit  averages  0.2 
hours  (12  minutes).  Therefore,  the  specified  mean  total  downtime  due  to  a 
failure  (MTTR)  shall  be  0.8  hours  (48  minutes)  for  the  gates  and  0.5  hours  (30 
minutes)  for  the  changers. 

The  analysis  focuses  on  the  gates,  investigating  the  effects  of  changes  in  the 
number  of  gates  and  their  reliability  and  maintainability.  Ten  cases  are 
examined  (the  data  for  the  gates  are  tabulated  in  Table  3): 

1.  The  base  case,  with  specifications  as  described  above. 

2.  The  gate  reliability  is  10000  MCBF. 

3.  The  gate  reliability  is  3000  MCBF. 

4.  The  gate  reliability  is  1000  MCBF. 


Case 

4, 

but 

with 

6 

gate 

units 

• 

Case 

4, 

but 

with 

a 

mean 

total 

downtime 

of 

0.6 

hours. 

Case 

4, 

but 

with 

a 

mean 

total 

downtime 

of 

0.3 

hours. 

Case 

4, 

but 

with 

a 

mean 

total 

downtime 

of 

0.2 

hours. 

Case 

4, 

but 

with 

a 

mean 

total 

downtime 

of 

0. 1 

hours . 

Case 

3, 

but 

with 

only  4 

gate 

units . 
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TABLE  3.  DATA  FOR  MULTIPLE-SERVICE  AREA  ANALYSES 


Gate 

Gate 

Arrival 

Number 

Processing 

Rate 

of 

Maintainability 

Rate 

(per 

Gate 

Reliability  (mean 

total  downtime) 

(per) 

Case 

hour) 

Units 

(MCBF)  (in  hours) 

hours) 

1 

5400 

5 

60000 

0.8 

1350 

2 

5400 

5 

10000 

0.8 

1350 

3 

5400 

5 

3000 

0.8 

1350 

4 

5400 

5 

1000 

0.8 

1350 

5 

5400 

6 

1000 

0.8 

1350 

6 

5400 

5 

1000 

0.6 

1350 

7 

5400 

5 

1000 

0.3 

1350 

8 

5400 

5 

1000 

0.2 

1350 

9 

5400 

5 

1000 

0.1 

1350 

10 

5400 

4 

3000 

0.8 

1350 

Underlining  indicates  changes  from  the  base 

case  (Case  1) 

The  results  for 

the  gates 

are  as  follows  (Table  4); 

TABLE  4.  RESULTS  OF  MULTIPLE 

-SERVICE  AREA  ANALYSES 

Case 

1 

2 3 4 

5 6 7 

8 9 

10 

Mean 

Queue  Length  3 . 2 

3.2  9.2  40.2 

3.3  28.0  20.4 

7.3  5.3 

(Congestion) 

Mean 

Passenger  Delay  0.5 

0.5  6.6  34.1 

2.6  24.2  16.7 

4.6  2.6 

(in 

seconds , 

exclud- 

ing 

processing  time) 

•Infinity  (queue  length  exceeds  500) 

(Note  — Except  for  Case  10,  the  mean  congestion  and  delay  at  the  coin  changers 
(not  the  gates)  remain  more-or-less  constant  over  all  the  cases  at  mean  queue 
lengths  and  mean  passenger  delays  of  approximately  6.0  and  3.0,  respectively.) 
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There  are  a number  of  conclusions  which  can  be  drawn  from  these  results: 

(Evaluation  of  given  situation  — Case  1)  No  serious  delay  problems  are 
expected  from  the  fare  collection  system  as  specified. 

(Sensitivity  analysis  of  gate  reliability  — Cases  1-4)  (Figure  5)  The 
specification  for  gate  reliability  can  be  significantly  reduced  from  its 
original  level  of  60000  MCBF  without  seriously  affecting  delay.  In  fact, 
the  reliability  can  decrease  by  almost  an  order  of  magnitude  without 
serious  impact.  Delays  start  becoming  significant  when  the  MCBF  reaches 
3000,  and  become  a problem  when  the  MCBF  reaches  1000. 

(Sensitivity  analysis  of  increased  number  of  gates  under  conditions  of  low 
reliability  — Cases  4 and  5)  Adding  one  additional  gate,  making  six  units 
in  all,  when  the  gate  reliability  is  low  (1000  MCBF),  is  equivalent,  in 
terms  of  the  resulting  delay,  to  entirely  solving  the  problem  of  low 
reliability,  i.e.,  improving  the  reliability  to  effectively  60000  MCBF. 
(Note  that  with  cost  data  showing  how  cost  varies  with  reliability,  an 
effective  tradeoff  analysis  can  be  made  to  decide  whether  adding  a gate  or 
improving  gate  reliability  would  be  more  cost-effective.) 

(Sensitivity  analysis  of  maintainability  under  conditions  of  low 
reliability  — Cases  4 and  6-9)  (Figure  6)  A delay  problem  due  to  low 
reliability  can  be  solved  for  this  system  by  improving  maintenance 
response,  but  the  improvement  must  be  consider-able  (even  an  improvement 
from  0.8  hrs.  to  0.1  hrs.  does  not  completely  restore  the  performance  of 
the  base  case).  Notice  that  the  nominal  detection  and  response  time  for 
the  system  is  0.2  hours,  so  that  Cases  8 and  9 can  only  be  brought  about  by 
an  improvement  in  this  time,  not  just  by  an  improvement  in  the  repair  time 
itself. 

(Sensitivity  analysis  of  decreased  number  of  gates  under  conditions  of 
marginal  reliability  — Cases  3 and  10)  The  system  cannot  operate  with 
fewer  than  five  gates.  If  failures  occur  under  a four-gate  operation,  the 
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MEAN  PASSENGER  DELAY  (SECONDS) 


SENSITIVITY  OF  DELAY  TO  GATE  RELIABILITY 

FIGURE  5.  GRAPH  OF  GATE  SENSITIVITY  TO  RELIABILITY 
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MEAN  PASSENGER  DELAY (SECONDS) 


MEAN  TOTAL  DOWNTIME  (HOURS) 


FIGURE  6.  GRAPH  OF  GATE  SENSITIVITY  TO  MAINTAINABILITY 
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system  will  sustain  catastrophic  congestion  and  delay.  (In  fact,  even 
without  the  impact  of  failures,  the  four-gate  system  is  operating  at 
capacity,  so  that  large  congestion  and  delay  can  arise  simply  due  to  random 
surges  of  passenger  arrivals . ) 

(Tradeoff  analysis  of  reliability  vs.  maintainability  under  conditions  of 
low  reliability  — Cases  3,  and  8)  An  increase  in  the  reliability  (of 
Case  4)  from  1000  to  3000  is  approximately  equivalent  in  delay  Impact  to  an 
improvement  in  the  maintainability  from  0.8  hours  to  0.2  hours. 

Remark  — The  precautions  described  previously  apply  of  course  to  the 
interpretation  of  this  example,  especially  the  precaution  on  the  volatility  of 
the  results  when  only  a few  breakdowns  have  occurred.  For  example,  with  an  MCBF 
of  3000,  an  average  of  1.8  breakdowns  would  be  expected  in  a peak  hour  in  which 
5400  passengers  arrive  (5400/3600=1.8).  Of  course,  however,  in  an  actual 
simulation  run  the  number  of  breakdowns  would  be  some  (Poisson-distributed) 
integer-valued  random  variable  of  0,  1,  2,  3,  or  more,  depending  not  only  on  the 
underlying  probability  distribution  for  the  breakdown  process,  but  also  on  the 
particular  sequence  of  random  numbers  on  which  the  specific  run  is  based.  The 
most  likely  number  of  breakdowns  to  occur,  2,  is  indeed  what  happens  in  the 
actual  sample  run.  However,  in  another  run  using  the  same  input  data,  the 
number  of  breakdowns  could  be  1 or  3,  since  each  run  is  based  on  a new  sequence 
of  random  numbers.  This  would  significantly  change  the  results. 

The  way  to  make  the  results  less  volatile  is  to  run  a number  of  independent 
peak-hour  runs  and  combine  the  results  (or  make  a single  run  for  a long  time 
period).  For  example,  if  five  independent  runs  were  made  of  the  peak-hour 
situation  with  an  MCBF  of  3000  (or  if  one  run  were  made  for  a five-hour  period 
under  peak  conditions),  the  expected  number  of  breakdowns  would  be  9.0.  To  have 
8 or  10  breakdowns  instead  of  9 would  clearly  affect  the  delay  result  much  less 
than  to  have  1 or  3 breakdowns  instead  of  2.  The  results  would  therefore  be 
less  subject  than  before  to  random  fluctuations  in  the  simulation  process.  To 
make  these  five  runs,  however,  would  use  five  times  as  much  computer  time.  For 
the  simulation  with  an  MCBF  of  3000,  this  would  mean  over  22  minutes  of  CPU  time 
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instead  of  under  5 minutes.  Since  the  cost  of  such  amounts  of  computer  time 
would  have  been  prohibitive,  only  one  1-hour  time  run  was  made  for  each  case 
examined.  In  a real  situation,  however,  either  longer  or  multiple  runs  must  be 
carried  out. 

5.3  The  Analytical  Model  (Single-Service-Area  Analysis) 


In  this  section  the  analytical  model  described  in  Section  3.3  is  used  to  examine 
the  coin  changers  in  the  Miami-system-based  example  discussed  above.  The 
analytical  model  considers  a single  service  area  (i.e.,  the  coin  changers) 
rather  than  the  complete  system.  However,  this  is  a reasonable  approach  in  this 
example,  since  the  coin  changers  receive  customers  directly  from  outside  the 
system,  not  from  any  other  service  area  within  the  fare  collection  system 
itself. 

Since  5400  passengers  per  hour  arrive  at  the  station  and  30/J  use  the  coin 
changers,  the  mean  arrival  rate  at  the  changers  is  1620  passengers  per  hour.  As 
described  in  the  previous  section,  there  are  six  changer  units  which  process 
passengers  at  a rate  of  360  passengers  per  hour  each,  with  a reliability  of 
80000  MCBF  and  a maintainability  of  0.5  hours  mean  total  downtime. 


Eight  cases  are  considered  (the  data  are  tabulated  in  Table  5): 


1.  The  base  case,  with  specifications  as  described  above, 

2.  The  reliability  is  8000  MCBF. 

3.  The  reliability  is  4000  MCBF. 

4.  The  reliability  is  3000  MCBF. 

5.  Case  4,  but  with  a mean  total  downtime  of  0.4  hours. 

6.  Case  4,  but  with  a mean  total  downtime  of  0.3  hours. 

7.  Case  4,  but  with  a mean  total  downtime  of  0.2  hours. 

8.  Case  4,  but  with  a seventh  coin  changer  unit  added. 
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TABLE  5.  DATA  FOR  ANALYTICAL  MODEL  ANALYSES 


Changer 

Changer 

Arrival 

Number 

Processing 

Rate 

of 

Maintainability 

Rate 

(passengers 

Changer 

Reliability 

(mean  total  downtime; 

(passengers 

Case 

per  hour) 

Units 

(MCBF) 

in  hours) 

per  hour) 

1 

1620 

6 

80000 

0.5 

360 

2 

1620 

6 

8000 

0.5 

360 

3 

1620 

6 

4000 

0.5 

360 

4 

1620 

6 

3000 

0.5 

360 

5 

1620 

6 

3000 

0.4 

360 

6 

1620 

6 

3000 

0.3 

360 

7 

1620 

6 

3000 

0.2 

360 

8 

1620 

7 

3000 

0.5 

360 

The  resulting  mean 
delay  is  obtained 

congestion  and  delay 
by  using  formula  (4), 

for  these  examples  are  shown 
Section  3*3): 

in  Table  6 (the 

TABLE  6. 

RESULTS  OF 

ANALYTICAL  MODEL 

ANALYSES 

Case  1 

2 

3 

4 

5 

6 

7 

8 

Mean  Queue  Length  5.9 
(Congestion) 

8.1 

18.0 

40.9 

12.8 

8.1 

6.8 

6.3 

Mean  Passenger  3.1 

Delay  (in 
seconds,  not 
including 
processing  time) 

8.0 

30.0 

80.9 

18.4 

8.0 

5.1 

4.0 

Some  of  the  conclusions  which  can  be  drawn  from  these  results  are: 


(Evaluation  of  given  situation  — Case  1)  The  coin  changer  service 
area  as  specified  should  not  have  serious  delay  problems. 
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(Sensitivity  analysis  of  coin  changer  reliability  — Cases  1-4)  (Figure 
7)  The  coin-changer  reliability  specification  can  be  significantly 
reduced  from  its  original  level  of  80000  without  significantly  degrading 
the  system.  A decrease  of  an  order  of  magnitude,  to  8000,  only 
marginally  increases  the  delay.  Delays  start  becoming  significant  when 
the  MCBF  reaches  4000,  and  become  a problem  when  the  MCBF  reaches  3000. 

(Sensitivity  analysis  of  coin  changer  maintainability  under  conditions 
of  low  reliability  — Cases  4-7)  (Figure  8)  An  improvement  in  mean 
maintainabilty  from  0.5  hours  to  0.4  hours  will  significantly  reduce  the 
mean  delay  resulting  from  low  reliability  (3000  MCBF).  Further 
improvements  beyond  this  occur  more-or-less  proportionally,  i.e., 
decreasing  the  maintenance  recovery  time  by  one-half  will  decrease  the 
mean  delay  by  one-half. 

(Sensitivity  analysis  of  increased  number  of  gates  under  conditions  of 
low  reliability  — Cases  4 and  8)  Adding  an  additional  unit  to  the 
system  when  the  reliability  is  low  (3000  MCBF)  solves  the  reliability 
problem,  i.e.,  the  mean  delay  becomes  similar  to  that  of  the  original 
system  (which  has  the  specified  reliability  of  80000  MCBF) . 

(Tradeoff  analysis  of  reliability  vs.  maintainability  under  conditions 
of  low  reliability  — Cases  1,2, 4, 6, and  7)  Improving  the  maintenance 
recovery  time  (of  Case  4)  from  0.5  hours  to  0.3  hours  is  equivalent  to 
improving  the  reliability  from  3000  MCBF  to  8000  MCBF.  Improving  the 
recovery  time  from  0.5  hours  to  0.2  hours  is  nearly  equivalent  to 
improving  the  reliability  from  3000  MCBF  to  the  full  original 
specification  of  80000  MCBF. 
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MEAN  PASSENGER  DELAY  (SECONDS) 


FIGURE  7.  GRAPH  OF  CHANGER  SENSITIVITY  TO  RELIABILITY 
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MEAN  PASSENGER  DELAY (SECONDS) 
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FIGURE  8.  GRAPH  OF  CHANGER  SENSITIVITY  TO  MAINTAINABILITY 


6.  THE  COST  MODULE 


We  have  dealt  so  far  with  the  analysis  of  the  passenger  congestion  and  delay  in 
a fare  collection  system,  given  the  configuration  of  the  station,  passenger 
demand,  and  the  processing  rates,  reliability,  and  maintainability  of  the 
equipment.  However,  one  would  often  like  to  have  some  idea  of  the  cost  of  the 
various  alternatives,  since  this  would  supplement  the  information  on  passenger 
dependability  performance  and  thus  provide  an  enhanced  analysis  of  the 
comparative  strengths  and  weaknesses  of  each  alternative.  To  that  end,  we  have 
developed  a cost  module  which  computes  the  annual  costs  relevant  to  fare 
collection  dependability  (these  include  equipment  acquisition,  spares  provision, 
equipment  operation,  and  scheduled  and  corrective  maintenance).  The  module 
accepts  input  data  on  equipment  acquisition  costs,  operation  costs,  discount 
rates,  spares  requirements,  annual  required  scheduled  maintenance,  pay  rates  for 
repair  personnel,  and  annual  passenger  volume,  and  uses  these  data,  along  with 
data  previously  supplied  to  the  model,  to  calculate  specific  and  total 
annualized  costs  at  the  station  under  considera-tion.  The  module  thus  makes 
possible  such  analyses  as  cost/performance  evaluations,  sensitivity  assessments 
of  costs  to  changes  in  specifications,  tradeoffs  between  cost  and  performance, 
and  tradeoffs  between  different  types  of  costs. 

6.1  Module  Description 

Among  the  large  variety  of  costs  incurred  in  operating  a transit  system,  five 
types  play  a role  in  fare  collection  dependability  analysis: 

Capital  costs  — This  represents  the  price  to  acquire  the  fare  collection 
equipment.  It  includes  the  direct  cost  of  the  equipment  itself,  as  well 
as  the  costs  of  financing,  specification  development,  procurement, 
delivery,  installation,  etc. 

Spares  costs  — In  order  to  assure  having  enough  functioning  equipment 
available  for  adequate  operation,  one  must  have  more  equipment  on  hand 
than  required  for  service,  with  the  excess  kept  as  spares  to  replace 
units  which  fail  in  the  field.  Spares  cost  is  the  cost  of  this  extra 
equipment.  The  requirement  for  spares  is  stated  as  a given  percentage  of 
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the  number  of  field  units,  and  so  spares  cost  will  be  a similar 
percentage  of  the  capital  cost. 

Operating  cost  — This  is  the  cost  to  operate  the  fare  collection 

equipment.  I-t  consists  of  such  costs  as  energy,  media  (such  as  tickets 
or  tokens),  routine  service  other  than  maintenance  (such  as  ticket 
filling,  ticket  removal,  coin-change  filling,  or  revenue  removal),  etc. 

Scheduled  maintenance  — The  equipment  must  be  maintained  periodically  to 
keep  it  in  good  operating  order.  Often,  such  maintenance  is  required  to 
keep  valid  the  warranties  on  the  equipment,  and  in  any  event  it  will  tend 
to  enhance  equipment  reliability,  thus  improving  performance  and 
decreasing  repair  costs. 

Corrective  maintenance  — This  is  the  cost  incurred  to  repair  equipment 
which  fails  while  in  service.  This  cost,  mainly  for  personnel,  transportation, 
and  parts,  clearly  depends  on  the  reliability  and  maintainability  of  the 
equipment,  as  well  as  the  number  of  units  and  passenger  demand  at  the  station. 

Note  that  the  above  costs  apply  to  a given  type  of  fare  collection  equipment. 

If  the  analysis  covers  multiple  service  areas  (coin  changers,  ticket  vendors, 
and  gates,  for  example),  there  will  be  several  sets  of  these  costs,  one  for  each 
type  of  service  area. 

The  cost  module  calculates  annual  costs  of  each  of  the  above  five  types  at  the 
station  under  consideration.  To  obtain  these  costs,  the  following  data  are 
required  (in  addition  to  the  information  already  provided  to  the  dependability 
model  Itself): 

1.  Acquisition  cost  per  unit: 

The  total  price  to  acquire  one  machine  unit  of  the  fare  collection 
equipment. 

2.  Useful  life  of  the  unit: 

The  number  of  years  the  unit  will  provide  adequate  service.  After  the 
useful  life  period  has  elapsed,  one  will  expect  to  have  to  replace  the 
unit,  thereby  incurring  a new  acquisition  cost. 
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3.  Discount  rate  (exclusive  of  inflation): 

The  "time  value  of  money";  the  ratio  of  the  value  of  one  dollar  now  to 
the  value  of  one  dollar  a year  from  now.  The  discount  rate  is  used  to 
convert  acquisition  cost,  a single  expense  incurred  at  the  start  of  the 
useful  life  of  the  equipment,  to  annualized  capital  cost,  the  equivalent 
annual  expense.  This  conversion  is  necessary  since  the  other  costs  in 
the  module  are  provided  and  calculated  on  an  annual  basis.  The 
annualized  capital  cost  represents  the  annual  repayment  which  would  be 
required  for  a loan  equal  to  the  acquisition  cost,  with  an  interest  rate 
equal  to  the  discount  rate  and  a repayment  period  equal  to  the  useful 
life  of  the  equipment. 

4.  Spares  ratio: 

The  number  of  additional  units  to  be  bought  as  spares,  expressed  as  a 
percentage  of  the  number  of  units  required  for  active  service. 

5.  Annual  operating  cost  per  unit: 

The  cost  for  energy,  media,  routine  service,  etc.,  for  a single  machine 
unit  for  one  year. 

6.  Annual  hours  of  scheduled  maintenance: 

The  number  of  labor-hours  needed  to  perform  the  required  scheduled 
maintenance  on  a single  unit  for  one  year. 

7.  Hourly  pay  rate  for  repair  personnel: 

The  costs  for  repair  personnel,  including  benefits,  supplies,  and 
overhead. 

8.  Annual  passenger  volume  at  the  station: 

The  total  volume  of  passengers  passing  through  the  fare  collection  area 
during  the  year.  This  determines  the  usage  of  the  equipment,  and  hence 
the  number  of  failures  and  consequent  corrective  repair  actions. 

The  five  cost  types  are  calculated  from  the  above  input  data  as  follows  (a 
dagger  "t"  indicates  data  from  the  dependability  model  itself,  rather  than  the 
above  data) ; 
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I.  Annualized  Capital  Cost 


ACAP  = ACQ  » (r  / (1-(1+r)-t))  * (n) 

where 

ACQ  = acquisition  cost 

r = discount  rate 

t = useful  life 

n = number  of  units  at  the  station"^ 


II.  Annualized  Spares  Cost: 


where 


SPRS  = s » (ACAP) 


s = spares  ratio 


III.  Operating  Cost: 


where 


OPER  = n * (UOPR) 

UOPR  = annual  operating  cost  per  unit 


IV.  Cost  of  Scheduled  Maintenance: 


where 


SCHD  = h * n * w 

h = annual  hours  of  scheduled  maintenance  per 

unit 

w = pay  rate  for  repair  personnel 


V.  Cost  of  Corrective  Maintenance: 


CORR  = (VOL)  » (p)  * (1/MCBF)  » (MTTR)  * (w) 
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where 


VOL 

P 


MCBF 

MTTR 


= annual  passenger  volume  at  station 
= passenger  split;  i.e.,  percentage  of  total 
passengers  who  use  the  service  area 

a 

containing  the  equipment  under  consideration' 
= Mean  Cycles  Between  Failures 

■4* 

= Mean  Time  To  Repair 


6.2  Sample  Cost  Analysis 

To  demonstrate  the  cost  module,  we  have  carried  out  sample  runs  based  on  the 
full-system  example  in  Section  5.2.  The  cost  and  other  additional  data 
necessary  for  the  cost  module  are  based  on  cost  information  for  the  faregates  at 
the  Metropolitan  Atlanta  Rapid  Transit  Authority.  (The  numbers  used  in  these 
examples,  however,  represent  hypothetical  situations  and  therefore  do  not 
reflect  actual  costs  at  Atlanta.  In  fact,  it  is  because  these  are  hypothetical 
situations  that  we  can  reasonably  combine  the  Atlanta-based  cost  data  with  the 
Miami-based  performance  and  demand  data.) 


The  passenger  arrival  and  processing  rates,  the  number  of  gate  units,  and  the 
gate  reliability  and  maintainability  are  the  same  as  those  provided  for  the 
full-system  example.  To  recapitulate,  these  values  are: 


Passenger  arrival  rate 
Passenger  processing  rate 
Number  of  machine  units 
Reliability 
Maintainability 

Percentage  of  passengers  using  gates 


= 5400  passengers/hour 

= 1350  passengers/hour 

= 5 

= 1000  MCBF 

= 0.8  hours  MTTR 

10056 


The  cost  and  other  additional  data  are: 


Acquisition  cost 
Useful  life 
Discount  rate 
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$33000/unit 
10  years 
1056 


Spares  ratio 
Annual  operating  cost 


= 5.5% 

= $41 00/unit 


Annual  hours  of  scheduled  maintenance  = 36  hours/unit 


Pay  rate  for  repair  personnel 
Annual  station  passenger  volume 


= $15. 66/hour 

= 2,448,000  passengers/year 


Eight  cases  are  considered  in  these  sample  analyses: 

1.  The  base  case,  with  data  as  shown  above  (this  case  is  equivalent  to  Case 
4 in  section  5.2) . 

2.  An  extra  unit  is  added. 

3.  A minor  improvement  increases  the  reliability  to  166?  MCBF,  while 
increasing  the  acquisition  cost  by  $1000. 

4.  The  scheduled  maintenance  is  doubled  to  12  hours  annually,  and  the 
equipment  reliability  improves  so  as  to  keep  the  total  annual  costs 
unchanged  from  the  base  case. 

(Cases  5-8):  The  following  reliability  improvements  occur,  and  the 

acquisition  cost  increases  so  as  to  leave  total  annual  costs  unchanged  from 

the  base  case. 

5.  The  reliability  improves  to  1667  MCBF  (as  in  Case  3 above). 

6.  The  reliability  improves  to  3000  MCBF. 

7.  The  reliability  improves  to  10000  MCBF. 

8.  The  reliability  improves  to  3000  MCBF,  and  the  useful  life  improves  to 


The  results  of  the  analysis  are  as  follows  (underlined  values  represent  the 
answers  sought  in  each  particular  case): 


15  years. 
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Case 

1 

2 

3 

4 

5 

6 

7 

8 

Reliability  (MCBF) 

1000 

1000 

1667 

1100 

1667 

3000 

10000 

3000 

Number  of  gate  units 

5 

6 

5 

5 

5 

5 

5 

5 

Sched.  maintenance 

hours 

36 

36 

36 

72 

36 

36 

36 

36 

Useful  life  (years) 

10 

10 

10 

10 

10 

10 

10 

15 

Acquisition  cost 

$33000 

$33000 

$34000 

$33000 

$47300 

$56900 

$65200 

$70400 

Annual  Costs 

Capital 

$26853 

$32224 

$27667 

$26853 

$38489 

$46301 

$53005 

$46279 

Spares 

1477 

1772 

1521 

1477 

2116 

2546 

2918 

2545 

Operating 

20500 

24600 

20500 

20500 

20500 

20500 

20500 

20500 

Scheduled  Main- 

tenance 

2819 

3383 

2819 

5638 

2819 

2819 

2819 

2819 

Corrective 

maintenance 

30688 

30668 

18398 

27880 

18398 

10223 

3067 

10223 

Total  annual  costs 

$82317 

$92647 

$70905 

$82348 

$82322 

$82389 

$82359 

$82366 

Mean  passenger  delay 

34.1 

2.6 

14.2 

25.6 

14.2 

6.6 

0.5 

6.6 

(seconds) 


From  the  above  analyses,  one  can  draw  a number  of  conclusions: 

(Case  1)  The  base  case  shows  a very  high  mean  passenger  delay  (this 
delay,  of  course,  is  the  same  as  Case  4 of  section  5.2).  Some  action 
needs  to  be  taken  to  reduce  delay.  Furthermore,  because  of  the  low 
reliability,  the  corrective  maintenance  costs  are  quite  high;  in  fact 
they  are  of  the  same  magnitude  as  the  capital  and  the  operating  costs. 


(Case  2)  By  adding  a sixth  gate  unit,  the  mean  passenger  delay  drops  off 
sharply,  becoming  virtually  insignificant.  However,  the  total  annual 
costs  rise  by  some  $10000,  or  12.5/6,  mostly  due  to  increased  capital  and 
operating  costs.  Despite  the  greater  number  of  machines,  corrective 
maintenance  costs  remain  the  same,  since  each  machine  handles  fewer 
passengers  (the  corrective  maintenance  cost  is  unchanged  because  it 
depends  on  total  annual  failures  at  the  station,  which  in  turn  depends  on 
machine  reliability  and  total  passenger  volume,  both  of  which  are 
unchanged  from  Case  1 ) . 
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(Case  3)  The  minor  equipment  improvement  is  definitely  worthwhile  at  a 
cost  of  $1000.  In  fact,  not  only  does  the  mean  passenger  delay  decrease, 
but  total  annual  cost  decreases  as  well,  since  the  savings  in  corrective 
maintenance  cost  due  to  fewer  failures  outweigh  the  increase  in  capital 
and  spares  cost.  Case  3 thus  represents  not  a tradeoff  from  Case  1,  but 
in  fact  a clear  improvement. 

(Case  4)  From  a strictly  cost  standpoint,  a doubling  of  the  scheduled 
maintenance  effort  is  worthwhile  if  it  yields  at  least  a 10^  increase  in 
reliability.  However,  such  a small  reliability  increase  still  leaves  a 
large  mean  passenger  delay. 

(Case  5)  The  minor  improvement  of  Case  3 represents  a clear  advancement 
even  at  costs  more  expensive  than  $1000,  up  to  a cost  of  $14,300. 
Therefore,  the  conclusion  of  Case  3 will  hold  up  even  against  severe  cost 
overruns. 

(Cases  5,6,  and  7)  As  the  equipment  reliability  improves,  the  corrective 
maintenance  costs  decrease.  One  can  therefore  sustain  increased 
acquisition  costs  in  order  to  obtain  this  improved  reliability  and  still 
have  the  situation  of  better  passenger  dependability  for  the  same  or 
lower  cost.  This  situation  of  absolute  improvement  holds  as  long  as  the 
acquisition  cost  remains  below  a given  cutoff  figure,  which  is  plotted  in 
Figure  9 against  the  associated  equipment  reliability  (for  acquisition 
costs  above  this  figure,  the  situation  becomes  a cost/performance 
tradeoff  rather  than  an  absolute  improvement).  For  an  improvement  in 
reliability  from  1000  MCBF  to  1667  MCBF  (a  range  of  improvement  possibly 
achievable  by  a small  enhancement  or  retrofit),  the  acquisition  costs  can 
increase  by  up  to  4355  while  still  yielding  an  absolute  improvement.  For 
an  improvement  in  reliability  from  1000  MCBF  to  3000  MCBF  purchase  of  new 
or  replacement  equipment),  the  acquisition  costs  can  increase  by  up  to 
72/5.  For  an  improvement  in  reliability  from  1000  MCBF  to  10000  MCBF  (a 
range  of  improvement  which  likely  would  need  an  advance  in  the  state  of 
the  art  compared  to  the  equipment  used  in  the  base  case),  the  acquisition 
costs  can  increase  by  up  to  97/5. 
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flCOUlSmON  COST  ($) 


(COMBINATIONS  OF  ACQUISITION  COST  AND  EQUIPMENT  RELIABILITY 
WHICH  GIVE  RISE  TO  EQUAL  TOTAL  ANNUAL  COSTS) 


EOUIVRLE.NT-COST  TRRDEOFF  CURVE 

FIGURE  9.  EQUIVALENT-COST  TRADEOFF  CURVE 
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Note  that  the  acquisition  costs  increasing  by  more  than  the  given  limits 
for  the  specified  reliability  does  not  imply  that  the  improved  equipment 
would  not  be  a better  buy  than  the  base-case  equipment.  Rather,  it  would 
simply  imply  that  total  annual  costs  would  increase.  Since  the  mean 
passenger  delay  is  significantly  better  for  the  improved  equipment  (the 
decreases  in  delay  are  67%,  85%,  and  98%,  respectively,  for  the  three 
cases),  the  additional  annual  cost  could  in  fact  be  worthwhile.  The 
analysis  in  this  case  would  show  the  tradeoff  between  annual  cost  and 
mean  passenger  delay,  and  thus  facilitate  the  decision  between  the  two 
choices  of  equipment. 

(Cases  6 and  8)  In  order  to  obtain  an  increase  in  the  useful  equipment 
life  from  10  years  to  15  years  (assuming  an  equipment  reliability  of  3000 
MCBF),  it  is  feasible  to  spend  up  to  a 24%  increase  in  the  acquisition 
cost.  Since  the  equipment  reliability  and  hence  the  mean  passenger  delay 
is  unchanged,  this  does  not  represent  a tradeoff  between  cost  and 
dependability,  but  rather  a clear  improvement  or  worsening  of  the  current 
situation,  according  to  whether  the  acquisition  cost  increases  by  less  or 
more  than  24%. 


52 


7.  SUMMARY  AND  FURTHER  DEVELOPMENT 


In  this  report,  we  have  discussed  the  concept  and  importance  of  fare  collection 
system  dependability  analysis,  as  well  as  the  types  of  such  analyses  and  their 
uses  to  transit  systems.  Analysis  software  consisting  of  simulation  and 
analytic  queuing  models  has  been  developed  to  assess  the  performance  and  cost  of 
fare  collection  systems.  Sample  analyses  based  on  data  derived  from  actual 
transit  systems  have  been  presented  to  demonstrate  the  use  of  this  software  for 
fare  collection  system  analysis.  (If  these  sample  analyses  had  been  those  of  an 
actual  transit  system,  such  a system  would  have  realized  that  they  could  reduce 
their  reliability  specifications  by  nearly  an  order  of  magnitude  with  little 
degradation  of  performance,  and  also  implement  a reliability-improving  retrofit 
while  recouping  the  expenditure  involved  through  lower  corrective  maintenance 
costs. ) 

This  report  and  the  accompanying  software  fulfills  the  development  of  the  fare 
collection  performance  model.  However,  the  modeling  effort  is  really  not 
complete  until  transit  systems  can  actually  use  the  model  to  carry  out  their 
fare  collection  system  analyses.  The  next  step,  therefore,  is  to  facilitate  the 
use  of  the  model  by  transit  operators  and  planners.  Towards  this  end,  we  shall 
enhance  the  model’s  efficiency  and  user-friendliness,  and  demonstrate  the  model 
by  applying  it  to  current  fare  collection  planning  efforts.  This  includes  the 
following  activities: 

User-friendliness 

Improvements  in  model  efficiency 
Interactive  menu  and  query  system  for  users 
Graphic  output 
User's  manual 

Obtaining  of  comments  and  feedback  from  potential  transit  users 
Model  Demonstration 

Support  of  fare  collection  system  planning  and  design  efforts 
Implementation  of  models  at  existing  transit  systems 
Holding  workshops  and  conference  presentations  on  the  model  and 
its  use 
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APPENDIX 


PROGRAM  LISTINGS  FOR  FARE  COLLECTION 
DEPENDABILITY  ANALYSIS  MODELS 


A-1/A-2 


A.1  LISTINGS  FOR  SIMULATION  MODEL 


Main  Program  MAIN. 

Subroutines 

DPTIME 

QUEUE 

BRAKDN 

RANDU 

BEGIN 

INPUT 

ENDING 

CSTMOD 


A-3/A-4 


CQMmgN7S1/D£P(3,50C),NCU5T(3),PASS(3,500) 

1 ,TuLAY(3,500),FUNIT(3,5CC) 

CaMMur.7S2/AP.?TM(3,500),NGUE:UE(3),rA3Q(3,500) 

1 , TULAyQ(3,50:-),FUNlTG(3,500) 

CntMMQN/33/REP(3,30  ),ND0WN(3) 

Cav,r.^nN/S4/I5,lMULT(1000  ),IQU£(10  00),IUNTT(1000)/lLOC2(10  00) 
C3MMCN/S'S/?JM,  ARPATE,NARR,SRVaAR,MCBF,RPRBAR,TIMEND,DLYMA>: 

1 /PPTIM,  A LT,Pr:OB,I.S£  ED, LAST, I COST/I  PRINT 

2 ,CUNIT,LIFE,DISCRT,SPARE3,f]PER,HRSM?4T,i\RATE,PASSYR 
ETMERSIUN  L AB EL ( 5 ), LSI AR { 51 ), QLST AR ( 1 01 ) 

I'lViENSION  C'JNIT(3),LIFE(3),SPARES(3),DPER(3),HRSMNT(3) 
DIMENSION  rtP,ATE(3),COST(3,6),QLPRQB(500) 
l.nEGcR  QLDI3TC3,0:50C) 

I El  EGER  PASSYR 

i PI  EGER  PAS3,?ASQ,A!JNIT,FU.NIT,FUNITQ 

integer  UiUE,QL(3),BLANk,STAR,QSTAR,QLSTAR,COUNT(3),NM(3), 

1 P K I N T,  M U LT  ( 3 ) , Di\ E, F R OR WH, TJ N I T,  MC BF  ( 3 ) , N P2  ( 3 ) , NP  ASS  ( 3 ) 

INTEGER  PERIOD, RPT 

REAL  LOSSPRC3),  TOTQL(3),  P ROB ( 3 ) , RPRBAR ( 3 ), TOTDLY( 3 ) 

PEAL  ULOAP(3),  SkVbAR(3),  SRVPTEO),  DLYBAR(3) 
data  label/ 'APR  TV', 'DEFRT', 'RPAIR', 'BRKD N ', 'CONTN '/ 

DATA  p P [ J M V:  H / 0 / , □ N E / 1 / 

DATA  NMUiNlT/3/,HLf  NK/'  '/,  ST  AR/ '^'/,(JSTAR/ '#'/ 

C READ  IT  THE  IwPUT  THROUGH  SUBROUTINES  BEGIN  AND  INPUT 


20  C ALL  P-EGI  'K  FROm,.jh,  N-»^,  ARR  » T’^',  NARR,SRVBAR,MCBF,  RPRBAR, 

1 TIMEND, DLYM AA,PRT1m, ALT,PR03,ISEED,PRINT,NMUNIT,CUNIT 

2 ,LIFL,dISCFT, SP AKES,QPER,HPSMVT,WRATE,PASSYR,ICQST,IPRINI) 


3  ‘M'ST  iir  TnE  1N->UT  VARIABLES  ( "M,  ARR  ATE,  N ARR,  SR  VBAR,  MC6F, 
: PPkHAI  ,1  IMEND,DLYN!  AX,pt>TlM,  ALT,PROR,ISEED),  ARE  LISTED 

C T.  S U t P 0 o’  T T tl  E ' fi  n G I N ' , 

: iO-ER  0/7'-  VARIABLES  AND  CONSTANTS: 

C TVER  ^>A.1E  DESCRIPTION 


C 

c 

c 

c 


1 

1 


c 

r 


c 


n 

r 


c 


K 

R 

R 

K 

1 

I 

r 

1 

1 

T 

1 


NEXT  ARRIVAL 

/.RrBAP  mean  arrival  TIM?  IN  SECONDS 
ARPTM  ARRIVAL  TIME 
BLANK  BLANK  VALUE  FOR  GRAPHS 
P^.K  2500  RANDOM  NUMBERS  FOR  BREAKDOWNS 

CIUNT  IN  CONJUNCTION  WITH  'MULT',  AIDS  IN  PRINTING  THE 
QUEUE  LENGTH  GRAPH  EVERY  10  SECONDS 
D^LTV  DELAY  TIME 
Ib.r'  DEPARTURE  ARRAY 

U'='P1  NEXT  DEPARTURE-IS  nsED  ONLY  FOR  PRINTING  PURPOSES 

AND  DOES  NOT  FIGURE  INTO  ANY  CALCULATIONS 

DLVy.AP  mean  delay 

DTIME  DEPARTURE  TIME  FOR  THAT  PASSENGER 
TXI'Ea  DETERMINES  IF  AN  ARRIVAL,  DEPARTURE  OR  REPAIR 
SHOULD  BE  PROCESSED  NEXT  (VALUES  ARE  1,  2 OR  3 
WITH  0 AT  THE  START) 

LA?  EL  ALPliAHF.TIC  ARRAY  OF  TITLES  FOP  THE  CURRENT 
PROCEDURE 

LEVENT  LAST  EVENT 

LuC  LOCATION  Ob  '*'  PITKIN  GRAPH  TO  BE  PRINTED  ON  UNIT  4 

LOCI  LOCATION  OF  'B'  WITHIN  GRAPH  TO  BE  PRINTED  ON  UNIT  7 

LOC2  LOCATION  OF  'i'  WTTMIN  GRAPH  TO  BE  PRINTED  ON  UNIT  5 
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C ] 

C :> 

C I 

C I 

C 

C I 

C I 

C 3 

C 1 

C 1 

C 1 

C 1 

C 1 

C I 

C 

C I 

c 

C I 

C K 

C K 

c 

c 

c i-; 

C 

c 

C ] 

C R 


r 

r* 


r* 

c 


■> 


l.uSi 

LC3SFP 

LSTA.-: 

"■'ULT 


NCUST 
N t.l  !*'  N 
N 'v'  U I T 

f.P  ASS 

N i u (.1  f: 

?4Sif,TCH 

PERIOD 

QL 

U L h A 
Q L.  T A P 
GST  5-. 

RRP 
r SP] 


PP'I  I.’iE 
Si^fVTp 
SR V PTE 

S T P 
T ‘M  IP 
TuTPLY 
T'ITrL 
b . I T 


LOSS  Of  SEPVICF 

LOSS  Uf  SERVICE  PROBABILITY 

GRAPH  FOR  UNIT  4 

IN  CONJUNCT  ION  V.TTH  'COUNT',  AIDS  IN  PRINTING  THE 
UUEUE  LENGTH  GRAPH  EVERY  10  SECONDS 

NUi^dER  of  CUSTCmERS-ONLY  used  IN  SUBROUTINE  'UPTIME' 
OIUMHER  Or  BREAKDOWNS-ONLY  USED  IN  SUBROUTINE  'BRAKDN' 
NUMBER  Of  UNITS 

PASSENGER  NUMBER  - USED  IN  MEAN  QUEUE  LENGTH 
PASSENGER  NUMBER- 1 NVOLVED  IN  MEAN  DELAY 
NUMBER  IN  QUEUE-ONLY  USED  IN  SUBROUTINE  'QUEUE' 

A SWITCH  USED  IN  SETTING  NEW  ARRIVAL  RATE 

PERIOD  POP  THE  QUEUE  LENGTH  GRAPH  (USUALLY  10  SECONDS) 

QUEUE  LENGTH  IK  SERVICE  OR  WAITING 

MEAN  QUEUE  LENGTH 

GRAPH  FOR  UNIT  b AND  UNIT  7 

t.  VALUE  FOf.  GRAPHS 

USED  FOR  placement  Of  # IN  GRAPHS 

REPAIR  ARRAY 

NEXT  kEPAIk-THIS  VARIABLE  IS  ONLY  USED  FOR  PRINTING 
PURPOSES  AND  DLES  NOT  FIGURE  INTO  ANY  CALCULATIONS 
ARRAY  OF  RANDOM  NUMBERS  FOR  BREAKDOWNS 
REPAIR  TIME 
SERVICE  TIME 

SERVICE  RATE-EACH  'MACHINE  SERVES  X NO.  OF  PASSENGERS 

PER  HOUR 

* VALUE  FOR  GRAPHS 

TIME  MINIMUM,  TI«E  ARRIVAL 

TOTAL  DELAY 

TOTAL  QUEUE  LENGTH 

THE  NUMBER  OF  UNITS  (OR  TYPES  OF  MACHINES); 

CURRENTLY  DEFINED  AT  3,  BUT  TO  RETURN  TO  A SINGLE 
UNIT,  JUS"^  GIVE  AN  ARRIVkL  PROBABILITY  OF  1.0  FOR 
unit  3,  AND  THE  otheR  2 UNITS  WILL  NEVER  BE  USED 


AST  — 0 
;:sRTC”=o 
Ns  r=  :< 

1 5-N 
I 'iDrx=o 
J = u 
K = C 

DO  45  UNIT=1,N’HUNIT 
:'JC1J5T(UNIT)  = 0 
i.OUEUE(  UMT)  = 0 
NDOWN(UNIT)=0 
NPASS(UNIT)=0 
QL(’JNTT)  = 0 
NPV(UNIT)=0 
T.;TDLY(UNIT)  = 0. 

T3TQL(UNIT)=0. 
f'.ULT(UNlT)  = 0 
C’JUNT(UNIT)  = 0 
URP(UNIT,1)=1,E9 
R':P(UNIT,1)  = 1.L9 

LOSSPRCUNIT  ) = MCBF  (UNI  T)/‘=’L0AT(MCBF  (UNIT  )-l) 
DO  45  T=0,d00 
.^LB1ST(UNIT,  I ) = 0 


A-6 


criKTiNUF: 


45 
C 

C GCT  RAM[)OM  MJmREC:  GENERATOR  STARTED 
C 

IF( lStED*Ea.99999)GQ  TC  51 
CALL  Rf.NDUC  1SELD,IX1,R) 

CALL  RANDU(ISEED+40,IX2,S) 

CALL  RAuDo( ISEED+dO,lX3,S) 

CALL  RANDU( IS5ED+12C,  IX4/S) 

CALL  RAN'DU(  IS£ED+lfeO,IX5/S  ) 

GO  TO  52 

51  TYRE  ^12 

512  F:ir-1AT{'  ENTER  LAST  KAADQM  NUMBER  SEEDS  FROM  PRIOR  RUN'/ 

1 7X'(FIVE  INTEGER  VALUES;  ON  ONE  LINE  SEPARATED  BY  SPACES)') 
ACCEPT  511,1X1,1X2,1X3,1X4,1X5 
511  fOi.MAT(51) 

CALL  RANDUCIXl, lY,k) 

1X1=1 Y 
C 

C DL'T’L  N-aPE  FIRST  ARRIVAL 
C 

52  tKRLAP=3b00./ ARRATE 
ARP--ARRBAR-^ALQG(l.-R) 

C CHMPUTr:  SERVICE  RATE 
C 

['0  44  UNIT=1,NMUN'IT 

Sr:;RTE(UNIT)=l.E9 

ir  ( SRVBAR(UNIT)  .EQ.O.  )GfJ  TO  44 

SRVRTE(UNIT)=3600./SRVBAP(UNIT) 

44  CONTU'UE 

C TMTIALIZE  GRAPn  ARRAYS 
DC  5 1=1,51 

5 LSTAr  (1  ) = BLAM’< 

DO  l^  1=1,101 
IB  CLoTA-(  I)=rir.  ANK 

1TE(4,302) 

3C2  FOhMATClHl,'  PASSMGR','  AT  ARR I V AL ', 2X 'SER VI CE ', 3X 'DELA Y ' 

1 /,^X,  'U  NT  T',  4X, 'TIMS ',5X,' ""IMS',  4X 'TIME'/) 

'.‘.R1TE(7,400) 

400  rORMATdMl,'  P A SS ',2X, 'UNIT  ', 'ARRIVAL ',2X' QUEUE '/1 5X'TI  ME 

1 4X, 'LENGTH'/) 

480  FGRMATdHl,'  TIME ',2  X' QUEUE', 4X, 'AT ',/7  X 'LENGTH  UNIT'//) 

wRlTF,(6,2C2) 

202  r Or'MATdHl,2X'TIME','  EV ENT ',  6X 'NEXT',  7X'NEXT',  5X 'NEXT' 

1 , 4X  'uU  E'J  E ',3X  'M«CHIN  LS  AT  '/1 9X  'ARRIVAL  ',  2X  'DEPARTURE  ' 

2 ,2X'REPATR',2X'LENGTH',2X,'IN  SERVICE  UNIT'/) 

1TL(8,490) 

490  EC.-MATC'  DEP  ARTURS ',  1 3X 'I  NIT  I AL ',  3 X'TOTAL '/ 

1 3 X'TIMf:',5X'P  ASS ENGER  ',4X 'UNIT', 5X 'DELAY'/) 

DO  19  LD=10,30,10 

Ln4=LD-*-4 

LD7=LD+7 

LD8=LD+6 

^^'■m’C{LD4,302) 

.':!aTS{LD7,400) 

WRITE(LDB,49C) 

19  CONTINUE 

C 
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PkIM  OUT  GkAPH  HEAbll^GS  FOR  INTERACTIVE  GRAPHS 
IE(PRTNT.I:;c.4)  TYPE  13C2 

1302  FOkMATCIHI,'  PASS  AT  ARRIVAL  SERVICE  DELAY  ’ ,J , 

1 ' UwlT  TIP.E  TIME  TIME  '/) 

IE(PkI NT. EQ.7)  TYPE  1400 

140C  EOKMATdnl,'  PASS  UN  ARRIVAL  0 UEU  E',/ , 1 1 X, 'T I ME  LENGTH*/) 

1 F( Pk I NT.EQ. 5)  TYPE  14b0 

1480  F Ji-:MAT(  Ihl,  ' TIME  QUEUE  A T*, /, 7X, 'LE  NG TH  UNIT*/) 

1 F(PRI NT.  Eul.  o )TYPE  20  2 
lF(pPINT.FQ.d)TYPE  490 

START  SIMULATION 

/vPtN  INDEX  = 1,  DO  ARRIVAL 

INDEX  = 2/  DO  DEPARTURE 
'.'.'HLr  IuDEX  = 3,  jO  REPAIRS 

LEVFNT=1NDEX 
I EA=2 

[ = A'"Ii,l(nFP(l,l),DLP(2,l),DEP(3,l)) 

F=AMIN1 (RE?(1,1),RLPC2,1),PEP(3,1)) 

IP( ARR.LT.D)INDEX=1 
T->1  J=i.MINl(  ARP,D) 

Ir CR.LE.TMIN)INDEX=3 

t^;t\=amin1(r,tmin) 

Di<A'\  A UNIT  F'JF  THE  FIRST  PASSENGER  TO  ARRIVE  AT 

NOTt.  Tn/.T  THIS  SECTION  IS  ONLY  USED  ONCE  AT  THE  VERY  BEGINNING 

Ir (L^VENT.GT.O)  GO  TO  42 
C:LL  RANnU(IX2/lY,S) 

I X2  = J Y 
NH  nIT=1 

IF(S.GT.PPOb(l  ))  NUNIT  = 2 
IF(  S.CT.  PRaS(l)-^PRuR(2))  MUNIT=3 
ij  F I T=  N u N I T 
Gi  TO  43 

C THIS  UNIT  IS  THE  ONE  ASSOCIATED  WITH  TMIN 
42  UNi'T'  = i 

If  (INDEX.  ) UNTT=NUNIT 

ir(  J NDEX.EQ.2.  AND.DEP(  2,1  ) .LT.DEPCl,  1 ) .'AND  .DEP  ( 2,  1 ) - LT.  DEP  ( 3, 1 ) ) 

1 UN1T=2 

1F( TNDEX.LQ.2. AND.DEP(3,1).LT.DEP(1,1). AND.DEP(3,1).LT.DEP(2,1)) 

1 U N I T=  3 

1F( Ti<^EX.EQ.3.AND.RFP(2,l).LT.REP(l,l).AND.REP(2,l).LT.REP(3,l)) 

1 u:;it=2 

IF( I NDEX. Eu. 3. A ID. HEP (3,1 ).LT. REP (1,1). AND. REP (3,1). LT. REP (2,1)) 

1 UFIT-3 

43  DEFI =DLPCUNIT,1 ) 

1F(DEP(UNIT,1 ).GT.9.£8)DEP1=0. 

EEP1=REP(UNIT,1) 

lF(RiP(UNIT,l).GT.9.E8)R»^Pl=0. 

1F( IIDEX.EQ.I )TOTQL(UHIT)=TOTQL(UNIT)+.5*NARR*(NARft-l) 

1 -tNARP*UL(  UMT) 

^'^ITE(o,201)TMlN,LABEL(lNDEX),  ARR,DEP1,REP1,QL(UNIT),NM(UNIT) 

1 , ■]  t:  i T 

IF( PRINT. EU.o)  TYPE  20 1,TMT n,L AD  EL ( I NDEX) , ARR,DEP 1, REPl, 
i QLCUNIT),  NM(UNIT),UNIT 

FCkmaT(1XF6.1,2XA5,3XF7.1, 4XF7.1,2XF7.1,4XI4,5XI4,5X,I4) 


201 


o o o noon 


C IS  IT  TI^-*l  FOk  Th£  SIMULATION  TO  END? 

C 

1FC^MIN.GT.T1MSND)G0  TO  < 

IS  IT  TIME  FOR  PRIMING  A 'COLUMN'  OF  THE  QUEUE  LENGTH  GRAPH? 
THIS  IS  SAMPLING  GRAPH,  THAT  IS,  IT  PRINTS  EVERY  10  SECONDS 

PER  100=10 

IFCTMIN  .LT.  MULl(UNiT))GG  TO  500 
IFtLiVENT  .EQ.  0)  UUE=QL(UNIT) 

IFCLEVENT  .EQ.  1)  QUE=QL(UNIT)-1 

If (QU£.EQ.-1)QUE=0 
1F(L£VENT  .EQ.  2)  QUL=tL(UNIT)  + 1 

IFCLlVENT  .eg.  3)  QU£=GL(UKIT) 

IF((TMIN-HULT(UNIT))  .LT.  PERIOD)  GO  TO  510 
LOC2=M1NO ((QUE  + i),101) 

GLSTA-'(L0C2  ) = GSTAR 

TilE  MAXIMUM  NUMBER  OF  CHARACTERS  FOR  THIS  GRAPH  IS  55 

IF (PFINT.  EQ.5.  A:';D.LnC2.GT.55)  QLSTAR(  55)  = QSTAR 
1I=M1NOCLOC2,55) 

If  (PRINT.  EQ. 5)  TYPE  1520,  »-'ULT(UNIT),  QUE,UNIT, 
i (^iLSTARdJ),  JJ  = 1,II) 

1520  F0RMAT(1X,I5,1X,T4,2X,I2,2X,55A1) 

1F(?RIIIT.EQ.5. AND.LQC2.GT.55)  QLSTAR( 55 )=BL ANK 

IFC'-'MUNIT.EQ.  l)GO  TO  91 

15=15+1 

IMULT( I5)=MULT(UNIT) 
iQU£CI5)=QUE 
iUkIT(15)=UMT 
IL[jC2(  I5)=LDC2 

520  t li.-‘-AT(1aI4,2(3X, 14), 3X1015.1) 

n GL3TAR(L0C2 )=BLANK 

COUMT(  U.NIT)  = C0UNT(UNIT)  + 1 
i-'.ULT(UNlT  )=CQUNT(UN'IT  ) * PFRIQD 
IF((TMTN  - HUL'i(L'NlT))  .RE.  PERIOD)  CO  TO  515 
51C  LOC2=M1MO (CQUE  + 1),101) 

1F(NMUMT,EQ.1)GU  to  92 
15=15+1 

1MULT(I5)=MULT(UN1T) 

1QUE(I5 )=QUE 
inNIT(lD)=UNIT 
1LUC2(I5)=L0C2 
52  QLSTAR(LGC2)=QSTAR 

IF(PkINT.EQ.5. AND.L0C2.GT. 55)  QLSTAR( 55)=QSTAR 
II=MIN0(LUC2,55) 

IF(PRINT.£Q,5)  TYPE  1520,  MULT ( UNI T ), QU E, UNIT, 

1 (QLSTAR( IJ), IJ=1,II) 

IE (Pkl NT. EQ.5. AND. LOC2.GT. 55)  QLSTAR(55  )=BLANK 

C!LSTAR(L0C2)=3LANK 

CnUNT(UNIT)=CUUNT(UNlT)  + 1 

MULT(UNIT)  = C0UNT(UMT)  •*  PERIOD 

C DETERMINE  IF  THE  CHANGE  IN  ARRIVAL  TIME  IS  TO  BE  MADE 

n 

V, 

300  IF(TMl£,.LE.PRTIM)GO  TO  8 

IF(NSATCH.SQ.1)G0  TO  8 
N5WTCH=1 
ARkATE=ALT 
A R K P A R = 3 6 0 0 . / A R P A T £ 
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ITE(  4,31C  )PkTIM,NARR,ARRATE 
IFCPRINT.Evj..-^)  TYPE  310,  PRTIM,  NARR,  ARR  ATE 
310  FTkMAT(//'  the  simulation  has  now  COVERED', F5.0, 'SECONDS. 

1 ARRIVALS  ARE  IN  GROUPS  OF', 14, 'AT  A RATE  OF ', F5 . 0, 'PER  HOUR') 

GO  TO  (1,2, 3), INDEX 

GENERATE  NEXT^  ARRIVAL 

call  RAMDUC 1X1, IY,S) 

IX1=IY 

ARR=ARR-ARRBAR*ALQG(1.-S) 

GEMlK'ATE  UNIT  FUR  NEW  ARRIVALS  AFTER  CURRENT  ARRIVALS  HAVE 

HEEN  PROCESSED 

NOTr.  THAT  IF  NARP  (NUMREP  OF  ARRIVALS)  IS  GREATER  THAN  1, 

SACi;  ARRIVAL  WILL  STILL  TAKE  PLACE  AT  THE  SAME  SERVICE  UNIT 

CALL  RANDU( 1X2, IY,S) 

1X2=1Y 
nUUT=1 

IF(S.GT.PROb(D)  NUNIT  = 2 
IF(S.GT.PRQR(1)+PPQB(2))  NUNIT=3 
C 

C PnOCESS  arrivals  in  GROUPS 

c 

DO  11  1=1, N ARP 

C This  IS  ONE  OF  TWO  PLACES  WHERE  'NP2'  IS  INCREMENTED; 

C IT  IS  a;  exact  count  Or  THE  NUMBER  OF  ARRIVALS  AT  EACH  UNIT. 

r 

V, 

C WHEN  THIS  N’UMHER  IS  PRINTED  OUT  AT  THE  END,  REMEMBER 
C T;.h'I  all  Or  THEM  MAY  NOT  bE  PRINTED  OUT  IN  THE  DETAILED 
C SUM.'iAKIES  OF,  FOP  EXAMPLE,  REPORT  3.  THIS  IS  BECAUSE 
C ^lEV  COT  PUT  IN  THE  ARRIVAL  QUEUE  (ARRTM)  AND  WHEN  'TMIN' 

C -!J  T :RF:  cutoff  time,  'lilEY  WERE  LEFT  HANGING  IN  THE  QUEUE. 

Np  T = Np  T + i 
Tnp:LAY  = 0. 

1F( NP^/500*500. EQ. NPT)TYPE  997, NPT 
997  FiVMAT(lXI4, ' PASSENGERS  HAVE  ARRIVED') 

HP2 ( UNIT )=NP2( UNIT)  + 1 
LOC1=MINO((QL(UNIT)+1),101) 

QLSTARCLOCl )=QSTAR 

aP1TE(7,401  )NP2(UNIT),UNTT,TMIN,QL(UNIT),(QLSTAR(II),II  = 1,LQC1) 
NQL=QL(UMT) 

QLDIST(UMT,NQL)=QLDIST(UNIT,NQL)  + 1 

I F( PRINT. EQ.7. AND.LOCl.GT. 49)  QLST AR( 49 )=QSTAR 

1I=M1NO(LOC1,49) 

JF(?R1NT.EG.7)  TYRE  1 401, NR2 ( U NIT  ), UNIT,TMIN, QL (UNIT ), 

1 (QLSTAR(IJ),IJ=1,II) 

14C1  FOKf-iAT(lX,  15,  I2,1X,F7.1,1X,I4,1X,49A1) 

If (PRINT.EQ.7. AND.LOCl .GT.49)  QLSTAR(49 )=BL ANK 
LD=7+10*UNIT 

1F( NMUNIT.GT. 1) WRITE (LD, 40 1)NP2(U NIT), UNI T,TMIN,QL( UNIT) 

1 ,(OLSTAR(II), 11=1, LOCI) 

401  F GRMAT( 1XI4,1X, 14,3XF  7.1,3X14,3X10141) 

QLSTAR(L0C1 )=BLANK 

1F(QL(UMT)  .GE.NM(UMIT))GQ  TO  12 
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C GSNhRATE  SERVICE  TIME 

C NG  DELAYS  FOR  THESE  ARRIVALS 

C 

CALL  RANDU(IX3,IY,U) 

I X3=IY 

SER  V TM=-  S RV  B A R ( U M T ) * ALOG  ( 1 . - U ) 

NO  uCLAYS  EOc  THESE  ARRIVALS 

DriME=TMlN+3LRVTM 
DELAY=0. 

NPASS(UNIT)=i',PASS(U!aT)-*-l 
LSTAk(l)=STAR 

TnTDLY(UNIT)  = TnTDLY(UNIT)-*-DELAY 

WP1TE(4,301)NPAS5(UNIT),HNTT,TM1N,SERVTM,DELAY,LSTAR(1) 
1F(PRINT.LQ.4)  type  1301,  NPASS(UNIT),  UNIT,TMI N, SERVTM, 
1 bELAY,(LSTARCIJ),IJ=l,l) 

1301  FOr.VAT(lX,I5,l2,3F6.1,lX,39Al) 

Ln=4+10*UNIT 

IF  ( NmUNI  T,  GT.  1 ) ViRITE  (LD,30 1)NP  ASS  (UNIT), UNIT, TMIN, SERVTM 
1 , DELAY, LSTAR(l) 

301  F(Jr.MAT(lXI5,2.X,I4,2X,3F8.1,2X51Al) 

LSTAR (1 )=3LANK 

call  DPTIME(l,DTlMh,,NPT,TDELAY,UNIT,UNIT) 

GO  TO  11 

n 

C POT  AN  ARRIVAL  INTO  THE  QUEUE  BECAUSE  THE  QUEUE  EXCEEDS 
C THE  NUVdEk  of  AVAILABLE  MACHINES  FOR  THAT  UNIT. 

12  CALL  QnEUE(l,TMIN,NPT,TDELAY,UNIT,UNIT) 

11  QL(UNIT)=QL(UN1T)+1 

GO  TO  5 

C p.-.OCFSS  .NEXT  DEPARTURE 
A WL(nNIT)=QL(UNIT)-l 

C To  IS  IS  THE  ONLY  jEI'AFTUPE  TIME  THAT  PASSENGERS  LEAVE  THE 
Z Dh.t'AR'iuRE  TIME  QUEUE 

t:.elay=tllay(unit,1) 

Jr!JNlT=FUNIT(UNIT,l) 

NPS=?ASS(UN1T,1) 

CALL  DPTIME(2,0.,0,0.,0,IJNIT) 


C CHECK  WHETHER  A BREAKDOWN  HAS  OCCURRED 

CALL  RAND’J(IX4,  IY,BRK) 

TX^=IY 

LaSS=LOSSPR(UNIT)»BRK 
IF(La3S.EQ.C)GO  TO  21 
TYPE  899, UNIT, TMIH 

899  FUnMAT(5X'BREAKD0WN  AT  UNIT', 16, AT  TIME  ',F8.1) 

: GE derate  REPAIR  TIME 

CALL  RANL'U{IX5,1Y,RPK) 

1X5=1 Y 

R P 1 1 M E = T M I ;N  - R P R E A R ( U N I T ) * A [ . 0 G ( 1 . - R P R ) 

C 

C A MACHINE  ENTERS  THE  BREAKDOWN  QUEUE 
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CtlL  BRAKD.Md,  kPTlME,UNIT) 

NM(I.lMT)  = Ny(UMT)-l 
1 F(  ^-IM(UNIT)  .LT.C)  TYPE  995 
995  POK'PAK'  NM  IS  NEGATIVE.  STOP.') 

IE(NM(UNIT) .LT.O)  STOP  10 
l)KPl=DSP(UNIT,l  ) 

IF (UEP( UNIT, 1 ).GT. 9. L8)DE?1  = 0. 

W'>-  TTE(6,'201)TMIu,LAnEL(4),  ARR,DEP1,REP(UNIT,1  ),QL(UNIT), 

1 NM(UNIT),nNIT 

If  (PRINT.  E'^.  6)  TYPE  20 1,  TMTN,L  ABEL  (4  ),  ARR,DEP1,  REP  (UNIT,  1 ) , 

1 QL(UMT),NM(UMT),UNIT 
GO  *^0  49 

21  IF(CL(UMT)  .LT.NM(UN1T))G0  TO  49 

C^LL  SANDUC 1X3,IY,V) 

1 X3=1Y 

S F R V TM=-  S R V B A R ( U MT  ) * AL on  { 1 . - V ) 

nTTMS=TMlN+S£RVTM 

I:  EL  t Y=T  y I N-  A R P TM  ( UNI  T , 1 ) 

f.PG  = PA3Q(UNIT,l) 

TwDLY=TDLAY(KUNIT,l  ) + DELAY 

iNP  ASS  (UNIT  )=NP  ASS  (UNIT  ) + l 

LGC=PIN1 ((DEL AY/DLYM AX* 50. >+1.99,51.) 

LSTAR(LUC)=STAR 

TOTDLY(UNTT)  = TOTDLV(UMT)+DELAY 

WP  TTE( 4,301 )NPASS( UNIT), UNIT, ARRTM( UN IT,1),SERVTM 
1 ,Dj.LaY,  (LSTAR(II),II  = l,LOC) 

1F(PPINT.EU.4.  A\'D.LOC.GT.39)  LSTAR  ( 39)  = STAR 
1 T=MINO(LOC,39) 

1F(^KINT. -Q.4)  TYPE  13C1,  NP ASS ( UNI T), U NIT, ARRTM( UNIT, 1) , 

1 S^PVTM,DELAY, (LSTAF(1J),IJ=1,II) 

Tr  ( pr  XNT.  t-C(.4.  AND.LOC.GT.3  9)  LSTAR  ( 39  )=  BL  ANK 
Lb=  1 + 1C*UNIT 

I F ( '’MU N I T . GT  . 1 ) Wk  I TE  ( LD,  30 1 ) NPASS  ( UNIT)  , UNI T 
1 , AkRTM( UNIT,1),SERVTM, DEL AY, (LSTAR (II), 11=1, LOG) 

LSI  i.R(LaC)  = ?LANF 

call  :'Pti‘^£(i,:jtime,npg,'^qdly,funitu(unit,1),unit) 

CALL  .^UEnE(2,0.,0,0.,0,UNIT) 

r' 

C RlIUG'v'E  ap  arrival  from  the  queue 

49  irCUMT.EU.  NMUNIT)  GO  TO  5 0 
UN1T=UMT+1 
D.EP1=DEP(  UNIT,1) 

1F(DEP(UNIT,1).GT.9.E8)DEP1=0. 

R‘-'Pl=kEP(UNIT,l) 

1F(PEP(UN1T,1).GT.9.E8)P='P1=C. 

I r:iTE(b,2Cl)  TkIN,LABEL(5),ARR,DEP1,kEP1,QL(UNIT), 

1 NM(UMT),  UNIT 

IF (pPTNT.Ey.b)  TYPE  201,  TMI N,L ABEL (5), ARR, DEP1,REP1,QL( UNIT), 

1 ^y(UNIT),  UNIT 
C 

C Tn  F.RE  HAS  BEEN  A DEPARTURE  AND  NOW  IT  IS  TIME  FOR  THIS 
C CUSTOMER.  TO  MOVE  INTO  THE  NEXT  UNIT.  HERE  DEPARTURE  TIME 
C I.S  THE  NEXT  ARRIVAL  TIME.  THIS  IS  STILL  VARIABLE  THIN, 

C M.)T  VARIABLE  ARR. 

NP2(UNIT)=NP2(UNIT)+1 

LGC1=MIN0((QL(UNIT)+1),101) 

QLSTAR.(LUC1  ) = QSTAR 

#*RITE(7,401)  NP2(UNIT),UNIT,TMIN,QL(UNIT),(QLSTAR(II),II  = 1,L0C1) 
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NQL=QL(UMT) 

QLDT3T(UNTT,NQL)=QLD1ST(TJNTT,NQL)  + 1 
iF(PRINI.E;f:<.7.AI\ID.L0Cl.GT.49)  QLSTAR(49)  =QSTAR 
II=M1MO(LOC1,49) 

IF(PkINT.EQ.7)TY?E  1401,  NP2  (UNIT), UNIT, TMIN,QL(UNIT), 

1 (QLSTAR(IJ),IJ=1,II) 

LL'  = 7+10*UNIT 

IFCMMUNIT.GT.l ) WRITE (LD, 401)  NP2 (U NIT), U NIT,TM I N, QL(UNI T) 

3 , (WLSTARdl  ), 11  = 1, LOCI) 

Ir(PRINT.EQ,7. ANr.LOCl.GT.49)  QLSTARC 49)=BLANK 
QLSTARCLOCl )=BLANK 
1F(QL(UNIT) .GE.NM(UNIT))  GO  TO  112 
CALL  RANDU(IX3,IY,Y) 

IX3=IY 

i)ERVTM=-SRVBAR(UMT)*ALOG(l.-Y) 

DT1ME=TMIN+SERVTM 

uELAY=0. 

Nr'A3S(UiNiIT)=NPASS(UNIT)+l 
LSTAkCl )=STAR 

TliTOLY(UNIT)  = TOTDLY(UNIT)  + DELAY 

WR1TS(  4,301)  N?ASS(UNIT),  UNIT, TMI N,SEP.VTM,DEL AY,LSTAR(  1 ) 
IF(PP,INT.EQ.4)  TYPE  1301,  NP  ASS  (UNIT  ),UNIT,TMI  N,SERVTM,  DELAY, 
1 (LSTAR(IJ),IJ=1,1) 

LD=4+10*UNIT 

1F(NMUNIT.GT.1)WRITE(LD,301)  NPASS(UNIT),  UNIT,TMIN,SERVTM 
1 ,OELAY,LSTAR(l ) 

LSTA^(1)=BLANK 

CALL  DPTIMSd, JTIM£,NPS,TDELAY,JFUNIT,UNIT) 
wL(UNIT)  = QL(UMT)  + l 
GO  T!J  9 
C 

C PUT  Sn^EONE  IN  THE  UUEUE 

112  CALL  UULUE(l,TMIN,NPS,TDELAY,JfUNIT,UNlT) 

C.L('JNIT)  = QL(UNIT)+1 

GH  TO  9 

5 0 LC)C  = N.aNl  ((TDELAY/DLYMAX’*5  0.  ) + 1.99,51.  ) 

L?TAR(LOC )=STAR 

ITE(B,308)TiMl  N,  NFS,  JFUNI T,  TDEL  AY,  (LST  AR  ( 1 1 ) ,I  I = l,LOC) 

3 OB  F0RMAT(1XF6.1,5XIS,9XI2,5XF7.1,2X51A1) 

Lr  = B+10*UMT 

1F(NMUNIT.GT.1)WRITE(LD,308)TMIN,NPS,JFUNIT,TDELAY 
1 ,(LSTAR(II),II=l,LOC) 

If  (PRII.T.EQ.8.AND.LOC.GT.39)LSTAR(39)=STAR 
II=MINO(LOC,39) 

IF (PRINT. EQ.8)TYPE  308,TPIN, NPS,JFUNIT,TDELAY 
1 ,(LSTAP(IJ),1J=1,II) 

LSTAP(39)=BLANK 
LSTAi<(LOC)  = BLANK 
GO  TO  9 

C PROCESS  NEXT  REPAIR 
C 

3 NM(UNIT)=NM(UNIT)+1 

TYPE  898, UNIT, TMIN 

898  FORMAKbX'REPAIR  AT  UNIT',16,',  AT  TIME  ',F8.D 

CALL  BRAKDN(2,0.,UNIT) 

C 

C A MACHINE  LEAVES  THE  BREAKDOWN  QUEUE 
C 


lF(QL(UNIT).LT.NM(UNIT))GO  TO  9 
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CALL  RANDUC 1X3, IY,V) 

IX3=IY 

SF.RVTM=-SRVnAR(UNIT)-*^ALOG(l.“Y) 
dtime=t.vin+si-:rvtm 
D:.LAY  = TMIN-  ARRTM(UMT,  1) 

NRR=PASQ(UNIT,1) 

TDELAY=TDLAYQ(UNIT,1)+DELAY 
fiPASS(Ul'.’IT)  = NPASS(UNIT)  + l 
LOC=MIN1C(DELAY/DLYMAX*50. )+1. 99,51. ) 

LSTAR(LQC)=STAR 

TG'IDLY(UN1T)=T0TDLY(UNIT)  + 0ELAY 

wRITE( 4,301 )NPASS( UNIT), UNIT, ARRTM( UNIT, 1),SERVTM 
1 , DELAY, (LSTAR(II), 11=1, LOC) 

IF(PRI NT.EU.4.AND.L0C.GT.39)  LSTAR(39)=STAR 
il=MlMO(LOC,39) 

IF(PHirJT.  EQ.4)  TYPE  1301 , flP ASS  (UNIT  ),UNIT,ARRTM (UNIT,1 ),  SERVTM, 
1 DELAY, (LSTAR(IJ),IJ=1, IT) 

IF(?P1NT.E(J.4.  AND.LOC.GT.3  9)  LSTAR ( 39 ) = BLANK 
LD=4+10*UN1T 

If (NMUNIT.GT.l) WRTT£(LD,301)NPASS(UNIT),UNIT,ARRTM(UNIT,l) 

I ,SERVTM,1)ELAY,(LSTAR(II),II  = 1,LQC) 

LSTAr<(LOC)  = BLA.‘iK 

C.VLL  DPTIME(1,DTIme,N?R,TDEL  AY,  FUNITQ(  unit,  1),  UNIT) 

P.E?!nVl-:  AN  ARRIVAL  FROM  THE  QUEUE 

CALL  QUEU>=:(2,0.,0,0.,0,UNIT) 

GO  TO  9 
L 

C END-Qr-JOB  PROCESSING 

C (INCLUDING  COMPUTATION  OF  MEAN  QUEUE  LENGTH  AND  MEAN  DELAY  TIME) 

4 LFUC3ST.EQ.0  )G0  TO  41 

C;  LL  CSTMUD(CUUT,LIFE,DTSCRT,SPARES,OPER,HRSMNT,WRATE 
1 ,PASSYR,NM,MCEF,RPRHAR, PROP, COST) 

41  DO  40  UN1T=1, 3 

I *-’(NM('jMT)  .EQ.0)WRTTE(9,9n4)UNIT 
i r ( ! j -1  ( U N I T ) . E Q . 0 . A ri  D . P K I N T . E Q . 9 ) T Y P E 9 0 4 , U N I T 

904  rC)KWAT(//'  UMT',12,'  IS  UNASSIGNED') 

IF(NM(UMT)  .EQ.O)GO  TO  46 
ULBAK(UN1T)=0. 

IF(NP2(UNIT).NE.0)QLBAR(UNTT)=TOTQL(UNIT)/NP2(UNIT) 

DLYBAR(UNIT)=0. 

IF(NPASS(UNIT).NE.O)DLYBAR(UNIT)=TOTDLY(UNIT)/NPASS(UNIT) 

W?1TE(9,303)UNIT,QLBAR(UNIT),DLYBAR(UNIT) 

3C3  FGRMAT(///'  UNIT', I3///4X'  MEAN  QUEUE  LENGTH  IS  ',F8.1 

1 /4X'  -MEAN  DELAY  IS  ',7XF8.1,'  SEC-') 

1F(  IPRINT.EQ.DTYPE  30  3,UNIT,QLBAR(UKIT),DLYBAR(UNIT) 

IMAX=-1 

DO  461  1=0,500 

461  IF(QLD1ST(UNIT, I) .GT.0)IMAX=I 

IF( IMAX.EQ.-l)GO  TO  462 
TPASS=NF2(UNIT) 

KMAX=MIN0(9,IMAX) 

WRITE(9,901)(K,K=0,KMAX) 

WR1TE(9,905) 

901  F0RMAT(//1X'  QUEUE  LENGTH  PROB AB ILITI ES '//2X10I 7) 

905  FOkMAT('  ') 

1F(  IPRINT.EQ.DTYPE  90  1,  ( K,K=0,  KMAX) 

IF(  IPRINT.EQ.DTYPE  905 
DO  4b3  II=0,IMAX,1C 
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III  = .MINO(  II+9,IMAX) 

DO  464  JJ=II,III 

GLPkOb(JJ)=GLDlST(UNIT,JJ)/TPASS 

464  CONTINUE 

WRITE(9,903)II,(QLPROB(JJ),JJ=II,III) 

IF(  IPRINT.EQ.DTyPE  90  3, II,  (gLPROB ( JJ  ),  JJ=I I,  II I ) 

903  F0WMAT(1XI3,1C(1X2FF5.1, '%')) 

463  CONTINUE 

ITc(9,906)(K,K=0,KMAX) 

IF(I?PINT.EG,1)TYPE  906, (K, K=0,KMAX) 

906  FORMATC//'  CUMULATIVE  QUEUE  LENGTH  PROB ABIL ITIES '// 2X1017 ) 

NRIT£(9,905) 

IFdPRINT.EQ.DTYPE  905 
DO  465  II=0,IMAX,10 
II1=mIN0(II+9,IMAX) 

DO  456  JJ=II,III 

I? ( JJ.EQ.OGO  TO  466 

QLPR0B( JJ)=QLPkOB(JJ-l )+QLPROB(JJ) 

466  CONTINUE 

NkITE(9,QC3 )II,(QLPROB(JJ),JJ=II,III) 

IFdpt^INT.EQ.DTYPE  903, II,  ( QLP ROB  ( JJ  ),  JJ=I I, III  ) 

465  CONTINUE 

NSFkV  = NP2 (U  NIT)-NPASS(UNTT) 

^*^ITE(9,203  ) NP2  ( UNIT  ),  NPftSS (UNIT ),NSLRV 
IF(  IPRINT.EQ.DTYPE  203,  NP2  ( UNIT),  NPASS(  UNIT)  , NSERV 
203  1-l)PMAT(/1XI5,'  passengers  ARRIVED  AT  SERVICE  AREA' 

1 /1XT5, • PASSENGERS  SERV ED -/I XI5, ' PASSENGERS  ' 

2 ,'NGT  YET  SERVED') 

462  IF(ICOST.EQ.O)GO  TO  46 

1F(L1FE(UNIT).EQ.0)GD  TC  46 

«vk1T£(  9,9C2)UNIT,  (CQST(UN’IT,I),I  = 1,6) 

!=’(  IPRINT.EQ.DTYPE  902,UNIT,  (COST(  UNIT,  I ),  1=1,  6) 

902  FO^:^AT(//'  COSTS  FOR  UNIT  ',14 

1 /7X'ANNUALIZED  CAPITAL  COST',T35, '$ ',F9.2 

2 /7X'ANNUALI7.ED  SPARES  COST', T35, '$ ',F9. 2 

3 /7x'U?ERATING  COST',T35, '$ ',F9. 2 

4 /7y.'SC:;EDULED  ^-MAINTENANCE  COST',T35, '$',F9.2 

5 /7^ 'CURP.SCTIVE  REPAIR  COST',T35, '$',F9.  2 

6 //7X'T0TAL  COST ',T34, '$ ',F10 . 2) 

46  CONTINUE 

TYPE  399,1X1,1X2,1X3,1X4,1X5 

399  r OR*^AT(//3X'  LAST  RANDOM  HUMBER  • SEEDS '/5  ( 4X115/ ) ) 

C IS  THERE  ANOTHER  RUN  TO  DO? 

C 

CALL  ENDING 
U (LA5T.EQ.DST0P 
GO  TO  20 
END 

SUBROUTINE  DP TI ME( IND EX,DTIME,P ASS,DEL AY,FUNIT, UNIT) 

C 

C DTIME  IS  THE  CURRENT  TIME  PLUS  SERVICE  TIME 
C DEPARTURE  TIME. 

C 

COMMON/S1/UEP(3,500 ), N CUST (3 ),P ( 3,500 ),D ( 3,500 ), U ( 3,500 ) 
INTEGER  UNIT,P,PASS,U,FUNIT 
GO  TO  (1,2), INDEX 

1 1F(NCUST(UNIT).EQ.O)GO  TC  14 

NC=NCUST(UNIT) 

DO  11  1=1, NC 
1SV  = I 
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IF(  Dc:P(  UNIT,  I).LT.DTIME)GO  TO  11 
GO  TO  12 
11  continue; 

14  NCUST(UNIT)=NCl)ST(UNIT)+l 

DLP(UNIT,NCU5T(UNIT) )=DTTME 
P(UN1T,NCUST(UMT))=PASS 
D(UNTT, NCU3T(UNIT) )=DELAY 
U(UNIT,NCUST(UNIT) )=FUNIT 
KSTUR^' 

12  NC^NCUSKUMT) 

DO  13  JJ=1SV,NC 
J = NCi;ST(UNIT)  + ISV-JJ 
DP:P(  UNIT,  J + 1)  = DEP(  UNIT,  J) 

P(UMT,  J + 1)  = P(UNIT,J) 

D(UN1T, J+1)=D(UNIT,J) 

UCUNIT, J+i)=U(UNIT,J) 

13  CONTINUE 

DEPCUNIT, ISV)=DT1ME 

P(UNIT,ISV)=PASS 

D(UNIT,  ISV)  = DF.LAY 

U(UNIT,ISV)=FUMIT 

NCUST(UNIT)=NCUST(UNIT)+1 

KETURN 

2 lF(f’CUST(UNIT).EU.l  )DEP(UNIT,1  ) = 1.£9 

IF( NCUST( UNIT ).EQ.l) RETURN 
N'C=f'!CUST(UNIT) 

DO  21  1=2, NC 

DF.P  ( UNIT,  I-1)=DEP(  UNIT,  I) 

PCUNiT, 1-1)=P(UNIT,I) 

D(U  riIT,I-l)  = D(UNJT,I) 

U(U NIT, I-1)=U(UN1T,I) 

21  CONTINUE 

NCuST(UNIT)  = .\CUST(UNlT)-l 

RETURN 

END 

SUSROUTINE  OUEUEC INDEX, ARTM,?ASS, DELAY, FU NIT, UNIT) 

: ^yRIVAL  QUEUE 

COPMQN/S2/ARRTM(3,500),NQUEUE(3),P(3,500),D(3,500),U(3,500) 
I. .T EGER  UNIT,P,PASS,U,FUNIT 
GO  TO  (1,2), INDEX 

1 IF(NwUEUE(UNIT).EQ.O)  GO  TO  13 
NU=N’QUEUE(UN1T) 

DO  11  1=1, NQ 

IF(ARPTM(UNIT,I).LT.ARTM)  GO  TO  11 
DO  12  JJ=1,NG 
J - NQ+ 1“ JJ 

PCUNIT, J+1 )=P(UNIT,J) 

D(UMT,J  + 1)  = D(UNIT,J) 

U(UN1T,J>1)=U(UNIT,J) 

12  ARRTM(UNIT,J+1  ) = ARRTM ( UNIT, J ) 

F (UNIT, I)=PASS 

DCUNIT, I)=DELAY 
U(UNIT,I)=FUMT 
ARRTM(UNIT,I)=ART^' 

NOUEUE ( UNIT )=NQUEUE( U NIT )+l 
RETURN 

11  CONTINUE 

13  NQU£UE(UNIT)=NQUEUl(UN IT)+1 
ARRTmCUNIT, NGUEUE(UNIT))=APTM 
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P(UNIT,NQUEUE( UNIT ))=F ASS 
D(UNIT,NQUEUE(UNIT))=DELAY 
U(UNIT,NQUEUE(UNIT))=FUNIT 
RETURN 

2 NQUEUE(UNIT  ) = NQUEUE(UNIT)-1 

1F(NQU£UE(UNIT) .EQ.O) ARRTM(UNIT,1)=1.E9 
IF  ( ‘{QUEUE  (U  NIT).  EQ.O)  RETURN 
NQ  = NQUEUE(1.'N1T) 

DO  21  1 = 1, ?JQ 
RCUNIT,I)=P(UNIT,I+1) 

D(UNIT,I)=D(UN1T,I+1 ) 

U(IJNIT,I)=U(UNIT,I+1) 

21  ARRTM(UNIT,I)  = ARRT^](UNIT,I  + l) 

AHR,  TM  ( UNI  T,NQUEDE  (UNIT  ) + l)=0. 

RETURN 
E ND 

SUBROUTINE  PR AKDN (INDEX, RPTIME, UNIT ) 
COMMiJN7S3/REP(3,30),NDOWS(3) 

INTF.UER  UNIT 
GO  TO  (1,2), INDEX 

1 IF( NDOWH(UNIT) .EQ.O)GO  TO  13 
ND=MDQ«‘N(UNIT) 

DO  11  1=1, ND 

IF(REP(UNIT,I) .LT.KPTIME)GO  TO  11 
DO  12  JJ=I,NUOWN(UNIT) 

J = NDOVi'N(UNIT)  + I-JJ 

12  FFP(UNIT,J+1)=RSP(UNIT,J) 

R£P(UNIT, I)=RPTIME 

xN  D 0 vx  N ( U N I T ) = N D 0 W N ( U N I T ) + 1 
RETURN 

11  CONTIMUS 

13  ND0WN(UNIT)=ND0WN(UNIT)+1 

R E P ( U N 1 T , N D 0 W N ( U N I T ) ) = R P T I M E 
RETURN 

2 IF(NDOhN(UNIT).Ea.l)KEP(UNIT,l)=l.£9 
I ? ( -.DOw  N ( UN  I T)  . EQ.  1 ) R ETURN 
N[)=NDG»iU(UNIT) 

l:  0 21  I = 1 , D- 1 

21  RLr(UHIT, I)=REP( UNIT, I+l) 

NDQWN(UNIT)=NDQWN(UNIT)-1 
RETURN 
END 

SUBROUTINE  RA NDU ( IX, I Y,YFL ) 

RANDOM  NUMBER  GENERATOR 

IY=IX*262149 
IF(IY)  5,6,6 

5 IY=IY+34359738337+l 

6 YFL=IY 

YFL=YfL*. 2910383E-10 
RETURN 
E N'D 

SUBROUTINE  BEGI N (FkOMfeH, NM , ARRATE,N ARR,SRVBAR,MCBF,RPRBAR, 

1 TIMEND,DLYMAX,PRTlM,ALT,PR09,ISEED,PRINT,NMUNIT,CUNIT,LIfE 

2 ,DISCRT,SPARES,OPER,HRSMNT,WRATE,PASSYR,lCOST,IPRINT) 
DOUBLE  PRECISION  IFN 1, 1FN2, IFN3, IFN 4, IFN5, IFN6 
INTEGER  FROMisH,  NM(  3 ) , MCPF  ( 3) / YES,NO,ANSR,PRINT 

1 , STORNM(3),UNIT,LIFE(3),PftSSYR 

REAL  SRVBAR(3),  RPRBAR(3),  PROB(3) 

REAL  CUNIT(3),SPARES(3),OPER(3),HRSMNT(3),WRATE(3) 
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DATA  YES,NO/1HY,1HN/,STORNM/0,0,0/,FILE/1HF/ 

c 

c 

C STRUCTURE  OF  THE  INPUT  AND  OUTPUT  FILE  QUESTIONS 
C 

C OUTPUT 

C 

C FIVE  CUES'tIQNS  ARE  ASKED  AT  THE  BEGINNING  OF  EACH  RUN 

C REGARDING  GIVING  FIVE  NAMES  OF  NEW  DISK  OUTPUT  FILES. 

C THESE  DISK  OUTPUT  FILES  CONTAIN  (IN  UNFORMATTED  WORDS), 

C ALL  THE  DATA  THAT  IS  PRINTED  IN  THE  GRAPHS,  ETC.,  AT  THE  FIRST 

C PART  OF  EACH  OF  THE  FIVE  REPORTS.  THE  REASON  FOR  STORING  THIS 

C IS  SO  THEY  CAN  BE  SORTED  LATER  AND  USED  FOR  PRODUCING 

C INDIVIDUAL  REPORTS  BY  UNIT-  AT  THE  END  OF  A RUN,  A QUESTION  IS 

C ASKED  ABOUT  WHETHER  THESE  FILES  SHOULD  BE  KEPT;  UNLESS 

C THEY  ARE  BEING  USED  AS  INPUT  TO  ANOTHER  PROGRAM  OR  NEEDED  FOR 

C MORE  COPIES  OR  COMPARISONS,  THEY  CAN  BE  DELETED. 

C 

C 

C INPUT 

C 

C FOR  THE  FIRST  RUN,  INPUT  MAY  BE  READ  IN  AS  ONE  RECORD  IN  FREE 

C FIELD  FORMAT  FROM  A DISK  FILE,  OR  THE  PROGRAM  CAN  ASK  FOR  DATA 

C TO  BE  TYPED  IN  AS  A RESPONSE  TO  EACH  OF  ITS  INPUT  QUESTIONS. 

C when  the  run  IS  FINISHED,  A QUESTION  IS  ASKED  ABOUT  WHETHER 

C ANOTHER  RUN  SHOULD  FOLLOW,  THE  USER  HAS  A CHOICE  Op  'YES' 

C OR  'NO',  AND  IF  'YES'  HAS  ANOTHER  CHOICE  REGARDING  INPUTTING 

C T^E  data.  if  the  data  IS  ON  DISK,  IT  MUST  BE  THE  NEXT 

C KRCORD  OF  THE  SAME  INPUT  FILE,  OR  THE  USER  CAN  CALL  FOR  THE 

C INPUT  PROMPT  SEQUENCE  TO  CHANGE  A FEW  VARIABLES  OR  ALL  VARIABLES. 

C ANy  TIME  THAT  INPUT  VARIABLES  ARE  CHANGED,  A LIST  OF  THEIR 

C CURRENT  VALUES  IS  PRINTED  BEFORE  EXECUTION. 


DO  600  UNIT=1,NMUNIT 
N M ( U M T ) = S T 0 R N M ( U N I T ) 
oOU  CONTINUE 

OPEN(’JMT=4,FILE='FOR04.DAT',DEVICE  = 'DSKS',ACCESS='S£QINOUT') 

Opt.  N(U  NIT-6, FILE='FQR0  6.  DAT', DEVICE='DSKS',ACCESS='SEQINOUT') 

OPE N(UNIT=7,FILE='FQR07. DAT', DEVI CE='DSKS',ACCESS='SEQINOUT') 

OPRN(UMT=8,FILS='f 0R08.DAT',DEVICE='DSKS',ACCESS='SEQINOUT') 

OPE N(UNIT=9,FILE= 'FOR09.DAT', DEV ICE='DSKS',ACCESS='SEQINODT') 

0PEN(UNIT=14,FILE='U14.DAT',DEVICE='DSKS',ACCESS='SEQIN0UT') 

0?EN(UNIT=17,FILE='U17-DAT',DEVICE='DSKS',ACCESS='SEQIN0UT') 

OP£N(UNIT=18,FILE='U18.DAT',DEVICE= 'DSKS ', ACCESS= 'SEQINOUT ') 

0?EN(UKIT=24,FILE='U24.DAT',DEVICE='DSKS',ACCESS='SEQIN0UT') 

QPEN(UNIT=27,FILE='U27.DAT',DEVICE='DSKS',ACCESS='SEQINOOT') 

UPEN(UNIT=28,FIL£='U28.DAT',DEVICE='DSKS',ACCESS='SEQINOUT') 

0PEN(UNIT=34,FILE='U34.DAT',DEVICE='DSKS',ACCESS='SEQIN0UT') 

0?EN(UNIT=37,FILE='U37.DAT',DEVICE='DSKS',ACCESS='SEQIN0UT') 

OPEN(UNIT=38,FILE='U38.DAT',DEVICE='DSKS',ACCESS='SEQINOUT') 

ICDNT=0 

IFCFROMWH.GT. 0)GO  TO  241 
TYPE  20 

20  FnPMAT(/20X'FARE  COLLECTION  PASSENGER-DELAY  MODEL') 

TYPE  294 

294  F0RMAT(//6X'THE  MODEL  PRODUCES  OUTPUT  FILES  SHOWING  MEAN  ' 

1 , 'QUEUE  LENGTH  AND'/6X'DELAY,  AND  QUEUE  LENGTH  ' 

2 , 'DISTRIBUTION  (AND  COSTS  IF  INCLUDED)-' 

3 //BX'SHOULD  THESE  BE  DISPLAYED  ON  THE  TERMINAL  AS  WELL?' 

4 /6X'("Y"  OR  "N"):  '$) 
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242  ACCEPT  7,ANSR 

IPRINT=-1 

1F(ANSR,EQ. YES) IPRINT=1 
IF(ANSR.EQ.NQ)IPRINT=0 
IF( IPRINT.GE.O)GO  TO  24 
iCONT=ICON'T-»-l 
IF(1C0NT.GT.3)ST0P  77 
TYPE  15,AN'SP. 

GO  TO  242 
24  TYPE  1 

1 F0RMAT(//,6X'THE  GRAPHS  PRODUCED  BY  THE  MODEL  ARE'// 

1 yy'l  — PASSENGER  DELAYS  BY  SERVICE  AREA'/ 

2 9X'2  — QUEUE  LENGTHS  AT  INTERVALS  OF  10  SECONDS'/ 

3 9X'3  — EVENT  LOG'/ 

4 9X'4  — QUEUE  LENGTHS  AT  PASSENGER  ARRIVALS'/ 

5 9X'5  — TOTAL  DELAY  PER  PASSENGER'/) 

TYPE  21 

21  FOkMAT(/,'  ARE  ANY  OF  THE  GRAPHS  TO  BE  ALSO  PRINTED', 

1 ' DIRECTLY  ON  THE  TERMINAL-',/,'  (ZERO  FOR  "NO"', 

2 ' OR  THE  GRAPH  NUMBER  (1,2, 3, 4 OR  5)  FOR  "YES"):  '$) 

25  ACCEPT  203,  PRINT 

PRlNT=PRINT+3 
203  fORMAT(I) 

TF{PHTNT.Eq.3.0R.PPINT.EQ.4.QR.PRINT.EQ.5.0R.PRINT.EQ.6. 

1 0k.PRTNT,EG.7.0R.PRINT,EQ*8)  GO  TO  241 
IC0NT=TC0NT^1 
IF( ICONT.GT.3)  STOP  12 
TYPE  310 

310  FUkV!AT('  THE  PREVIOUS  RESPONSE  WAS  ',15, 

1 /,'  PLEASE  ANSWER  WITH  ONE  INTEGER  VALUE', 

2 'EITHER  0,  1,  2,  3,  4,  5:  '$) 

GO  TO  25 

241  1F(FR0MWH.EQ.2)  GO  TO  13 

IF (FROMftH.EQ.O)GO  TO  19 
TYPE  27 

27  F0FMAT(6X'IS  THE  NEW  DATA  TO  COME  FROM  A NEW  INPUT  FILE?'/ 

1 7X'(THE  ALTERNATIVE  IS  TO  vaKE  CHANGES  IN  THE  CURRENT  INPUT' 

2 DATA)'/oX'("Y"  OR  "N"):  '$) 

271  ACCEPT  7,ANSR 

TF( ANSk.EQ. YES)FPOMWH=0 
IFiANSR.EQ. YES)GO  TO  19 

IF ( ANSR.EQ. NO) CALL  INPUT(FR0MWH,NM,ARRATE,NARR,SRVBAR,MCBF, 

1 rprbar,timend,dlymax,prttm,alt,prob,iseed,print,nmunit, 

2 STORNM, CUNIT, LIFE, DISC RT,SPARES,0PER,HRSMNT 

3 ,WRATE,PASSYR,ICOST,IPRINT) 

IF( ANSR.EQ. NO)RETURN 
ICONT=ICONT+l 
1F(IC0NI.GT.3)ST0P  1 

TYPE  15,ANSR 
GO  TO  271 
1<3  ICONT  = 0 

TYPE  6 

6 FORMAT!/'  ACCEPT  INPUT  FROM  THE  TERMINAL-',/, 

1 ("Y"  UR  "N"):  '5) 

119  ACCEPT  7,  ANSR 

7 FORMAT(Al) 

IF(ANSR.EQ.YES)  GO  TO  14 
IF( ANSR.EQ. NO)  GO  TO  30 
ICONT=ICONT+l 
IF(ICONT.GT.3)  STOP  1 
TYPE  15,ANSP 
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15  FOKMATC/'  THK  PREVIOUS  RESPONSE  WAS 

1 Al,"',  PLEASE  ANSWER  "Y " OR  "N":  '$) 

GO  TO  119 

14  CALL  INPUT(FROMWH,NM,ARRATE,NARR,SRVBAR^MCBF, 

1 RPRLAR,TlMEND,DLYMAX,PRTTMr ALT, P ROB, ISEEU, PRINT, NMONIT, 

2 STORNM,CUHTT,LIFE,DISCRT,SPARES,OPER,HRSMNT 

0 ,WHATE,PASSYR,ICOST,IPfiINT) 

FROMKH=l  ’ 

RETURN 

13  ICUNT=0 

IF(FROMWH.EQ.O)  GO  TO  30 
TYPE  28,IFN5 

28  FORMAT!/,'  IS  THIS  NEW  DATA  EITHER  A NEW  FILE  OR' 

1 ,'  THE  NEXT  RECORD  OF  FILE  ', A10/6X '( "F " FOR  NEW  FILE', 

2 ',  "Y"  FOP  NEW  RECORD,  OR  "N " FOR  NEITHER);  '$) 

40  ACCEPT  7,  ANSR 

IFCANSR.EG. YES)  GO  TO  29 

ir(ANSK.EQ.Na)  CALL  1 NPUT( FROMWH, NM, ARR ATE, NARR,SRVBAR, MCBF, 
kPRBAR,TlMEND,DLYMAX,PRTIM,ALT,PROB,ISEED,PRINT,NMUNIT, 

2 STORNM,CUNlT,LIFE,DISCRT,SPARES,OPER,HRSMNT 

3 , WRATE,PASSYR,ICOST,IPRINT) 

IF(ANSR.EQ.NU)  RETURN 

IF( ANSR.EQ. FILE)FROMWH=0 

I^( ANSR.EQ.FILE)GO  TO  19 

ICQNT=ICONT+l 

IF( 1CONT.GT.3)  STOP  4 

TYPE  15,  ANSR 

GO  TU  40 

: NOTE  THAT  IF  A BAD  RESPONSE  IS  GIVEN  TO  THIS  REQUEST,  THAT  IS, 

C IF  ThL  PROGRAM  CANNOT  FIND  THF  SPECIFIC  DATA  FILE  ON  DISK,  THE 
: SYSTE*^  i^ILL  GENERATE  A REQUEST  FOR  A NEW  INPUT  FILE.  HOWEVER, 

C THE  InCURkECT  FILE  NAME  WILL  STILL  BE  PRINTED  IN  THE  SECTION 
: REGARDING  THE  CURRFNT  INPUT  FILES,  BECAUSE  THE  PROGRAM 

: DIDN'T  GENERATE  THE  NEW  REQUEST,  THE  SYSTEM  DID. 

30  TYPE  22 

22  FOP''AT(/'  ENTER  INPUT  FILE  NAME  ',/, 

1 ' (UP  TU  10  CHARACTERS);  ',$) 

ACCEPT  23,IFN5 

23  FORFAT(AIO) 

IF  (FROMWH.  EQ.2)  CLOSE (UMT=55) 

0PEN(UN1T=55,FILE=1FN5,DEVICE='DSK',ACCESS='SEQIN') 

TYPE  3,  IFN5 

3  F0RMAT(6X'THE  INPUT  FILE  IS;  ',A10) 

29  ICU3T=0 
ICQNT1=0 
TYPE  401 

401  FOPXAT(//6X'ARE  COSTS  INCLUDED  IN  THIS  RECORD?' 

1 /6X'("Y"  OR  "N");  '$) 

132  ACCEPT  7, ANSR 

IF(ANSR.EQ. NO)GO  TO  291 

1F(ANSR.EQ.YES)G0  TO  292 

ICO»Tl=TCONTl+l 

IF(  ICONTl.GT.3)STOP  2 

TYPE  15, ANSR 

GO  TO  132 

291  READ(55,*)ST0RNM,ARRATE,NARR,SRVBAR,MCBF,RPRBAR,TIMEND, 

1 DLYMAX,PRTlM,ALT,PROB,ISEED 

GO  TO  293 

292  READ(55,*)STORNM,ARRATE,NAPR,SRVBAR,MCBF,RPRBAR,TIMEND, 
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1 DLYMAX,PRTI.M,ALT,PRQB,ISEED,CUNIT,LIFE,DISCRT^SPARES,OPER, 

2 HRSMNT,WRATE,PASSYR 
ICQST-1 

293  CLDSE(UNIT=55) 

DO  16  UNIT=lvNMUNIT 
NM(UNIT)=STaRNrt(UNIT) 

16  CONTINUE 
FPaMWH=2 

c 

C ARRATt  (MEAN  ARRIVAL  RATE)  IS  NOT  AN  ATTRIBUTE  OF  EACH 
C MACHINE  OR  SERVICE  UNIT,  THEREFORE  IT  IS  NOT 
C DIMENSIONED  BY  UNIT 
C 

C INPUT  IS  FREE  FORMAT 
C 


c 

Cm, 

VARIABLE 

NAME 

TYPE 

DESCRIPTION 

C 

1-3 

im 

INT 

THE  NUMBER  OF  MACHINES  FOR 

C 

c 

EACH  SERVICE  UNIT 

c 

c 

4 

ARRATE 

REAL 

MEAN  ARRIVAL  RATE  PER  HOUR 

c 

B 

NARR 

INT 

THE  NUMBER  OF  ARRIVALS  IN  A 

c 

r 

GROUP 

c 

r 

SRVBAR 

REAL 

MEAN  SERVICE  TIME  IN  SECONDS 

c 

9-11 

MCBF 

INT 

MEAN  CYCLES  BETWEEN  FAILURES 

c 

r 

FOR  EACH  SERVICE  UNIT 

c 

r 

12-14 

RPRBAR 

REAL 

MEAN  TIME  TO  REPAIR  IN  SECONDS 

Lr 

C 

15 

TIMEND 

REAL 

TIME  THE  SIMULATION  ENDS  IN 

SECONDS 

U 

16 

DLYMAX 

REAL 

THE  MAXIMUM  DELAY  WHICH  CAN  BE 

r 

REPRESENTED  BY  THE  GRAPH 

o 

w 

17 

PRTIM 

REAL 

THE  TIME  IN  SECONDS  WHEN  THE 

n 

NEW  ARRIVAL  RATE  WILL  OCCUR 

o 

IB 

ALT 

REAL 

THE  ALTERNATE  MEAN  ARRIVAL 

U 

r 

TIME  IN  SECONDS 

C 

19-21 

PROP 

REAL 

THREE  PROBABILITIES  OF  ARRIVAL 

n 

c 

FOR  THE  THREE  SERVICE  UNITS 

22 

I SEED 

INT 

RANDOM  NUMBER  SEED 

CALL  INPUT(PROMWH,NM,ARPATE,NARR,SRVBAR,MCBF,RPRBAR, 

1 1IMEND/DLYMAX,PRTIM,ALT,PR0B,ISEED,PRINT,NMUNIT,ST0RNM,CUNIT 

2 ,LIFE,DISCRT,SPARES,OP£R,HRSMNT,WRATE,PASSYR,ICOST,IPRINT) 

RETURN 

E'ND 

SUBROUTINE  INPUT(FROMWH,NM,ARRATE/NARR,SRVBAR,MCBF,RPRBAR, 

1 TImEN'D,DLYMAX,PRTIM,ALT,PROB,ISEED,PRINT,NMUNIT,STORNM, 

2 CUNIT,LIFE,DISCRT,SPARES,OPER,HRSMNT, WRATE,PASSYR,IC0ST,IPRINT) 
INTEGER  PRINT,FR0MWH,CHANGE,YES,N0,ANSR,NM(3),MCBF(3) 

1 ,S1QRNM(3),UNIT,LIFE(3),PASSYR 

PEAL  SPVBAR(3),RPRBAR(3),PR0B(3),CUNIT(3),SPARES(3),0PER(3) 

1 ,HRSMNT(3), WRATEC3) 
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1 

101 

1 

5U 

2 

102 

1 

3 

1C3 

1 

4 

104 

1 

3 

105 

1 

t 

lOn 

1 

7 

107 

1 

B 

lOd 

2 

1 

9 

109 

1 

2 

3 

10 

110 

1 

2 


DATA  YES,  NQ/IHY/IHN/ 

CHANGE,=0 

1CUNT=0 

1 F(fROMWH.GT-.0)  GO  TO  667 

1CIJNT=G 

TYPE  101 

FORMAT(/'  1.  ENTER  NUMBER  OF  MACHINES  FOR  EACH', 

' SERVICE  UNIT-  (THREE  INTEGER  VALUES):  '$) 

ACCEPT  201,  STORNM 
DO  50  UN1T=1, NMUNIT 
NM( UNIT)=STOPNM(UNIT) 

CONTINUE 

IF(ICONT.GT.O)  GO  TO  300 
TYPE  102 

FORMATC/'  2.  ENTER  MEAN  ARRIVAL  RATE  PER  HOUR-',/, 

' (ONE  REAL  VALUE):  '$) 

ACCEPT  202,ARRAT£ 

IF(ICONT.GT.O)  GO  TO  300 
TYPE  103 

FORMAT(/'  3.  ENTER  NUMBER  OF  ARRIVALS  IN  EACH  ', 

'GROUP-',/,'  (ONE  INTEGER  VALUE):  ',$) 

ACCEPT  203,  NARR 

IFC  ICONT.GT.O)  GO  TO  300 

'^YPE  104 

FORMAT(/'  4.  ENTER  MEAN  SERVICE  TIME  IN  SECONDS-', 

/,'  (THREE  REAL  VALUES):  '$) 

ACCEPT  202,  SRV8AR 

IFC ICONT.GT.O)  GO  TO  300 

TYPE  105 

FORMATC/'  5.  ENTER  MEAN  CYCLES  BETWEEN  FAILURES-', 

/,'  (THREE  INTEGER  VALUES):  ',$) 

ACCEPT  201,MCBF 

IFC ICONT.GT.O)  GO  TO  300 

TYPE  105 

fORMATC/'  6.  ENTER  MEAN  TIME  TO  REPAIR  IN  SECONDS-', 

/,'  (THREE  REAL  VALUES):  ',$) 

ACCEPT  202,  RPR3AR 
IFCICONT.GT.O)  GO  TO  300 
TYPE  107 

FORMATC/'  7.  ENTER  TIME  TO  END  THE  SIMULATION  IN  ', 
'SbCuNDS-',/, ' (ONE  REAL  VALUE):  ',$) 

ACCEPT  204,TIMEND 

IFC ICONT.GT.O)  GO  TO  300 

TYPE  108 

FORMATC/'  8.  ENTER  THE  MAXIMUM  DELAY  TIME  WHICH  CAN  BE  ' 

, 'REPRESENTED'/6X'1N  THE  DELAY-TIME  GRAPH,  IN  SECONDS-' 

,/,  ' (ONE  REAL  VALUE):  ',$) 

ACCEPT  204,  DLYMAX 
IFCICONT.GT.O)  GO  TO  300 
TYPE  lOy 

FORMATC/'  9.  ENTER  THE  TIME  IN  SECONDS  WHEN  THE  NEW', 

' MEAN  ARRIVAL  RATE  OCCURS'/ 

eX'CIF  THE  MEAN  ARRIVAL  PATE  DOES  NOT  CHANGE,  ENTER  99999.)-'/ 
' (ONE  REAL  VALUE):  ',$) 

ACCEPT  204,  PRTIM 
IFCICONT.GT.O)  GO  TO  300 
TYPE  110 

FORMATC/'  10.  ENTER  THE  NEW  MEAN  ARRIVAL  RATE', 

' IN  SECONOS'/6X'(IF  THE  MEAN  ARRIVAL  RATE  DOES  NOT  CHANGE,  ' 
,'PRESS  RETURN)-'/,'  (ONE  REAL  VALUE):  ',$) 

ACCEPT  204, ALT 
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11 

111 


12 

112 


A 

1 


401 

132 


131 

13 

113 


14 

114 


15 

115 


16 

116 


17 

117 


18 

118 


19 

119 


20 

120 


IF( ICOKT. GT.O)  GO  TO  300 
TVPE  111 

FORMA'K/'  11.  ENTER  THE  ARRIVAL-SPLIT  PROBABILITIES-', 

/,'  (THREE  REAL  VALUES):  ',$) 

ACCEPT  202,PRQb 

IF( ICONT.GT.O)  GO  TO  300 

TYPE  112 

FOPHATC/'  12.  ENTER  THE  RANDOM  NUMBER  SEED-', 

/6X'(USE  "9999«^"  TO  INDICATE  LAST  SEEDS  FROM  PRIOR  RUN)', 
/,'  (ONE  I,y|TEGER  VALUE):  ',$) 

ACCEPT  203,  ISEED 

IF( ICONT.GT.O)  GO  TO  300 

TYPE  401 

FORVAT(//7X'DO  YOU  WISH  TO  INCLUDE  COSTS  IN  THIS  RUN?' 
/7X'("Y"  OR  "K"):  '$) 

ACCEPT  302,ANSR 
1COST=0 

IF(AN3R.EQ. NQ)GO  TO  667 

IF (ANSR.EQ. Y£S)GO  TO  131 

ICn^!Tl  = ICONTl  + l 

1F( ICONT1.GT.3) STOP  2 

TYPE  303,ANSR 

CO  TO  132 

1CQST=1 

TYPE  113 

FORMAT!/'  13.  ENTER  THE  CAPITAL  COST  PER  UNIT' 

/7X'(T}IREE  INTEGER  VALUES):  ',$) 

ACCEPT  202,CUN1T 
if(tcont.gt.u)go  to  300 

TVpE  114 

F0kMAT(/'  14.  ENTER  THE  USEFUL  LIFE  OF  THE  UNIT' 

/7X'(THREE  INTEGER  VALUES):  ',3) 

ACCEPT  201, LIFE 
IF(ICUNT.GT.O)GO  TO  300 
TYPE  115 


FOKMAT!/'. 15.  ENTL^  THE  DISCOUNT  RATE,  IN  PERCENTAGE  TERMS 
/7X'(0NE  REAL  VALUE):  ',$) 

ACCEPT  204,D1SCRT 

IF(  ICUNT.GT.O)GO  TO  300 

TYPE  116 


FORMAT!/'  16.  ENTER  THE  SPARES 
/7X'(THREE  REAL  VALUES:)  ',$) 

ACCEPT  202,SPARES 
IF(ICONT.GT.O)GO  TO  300 
TYPE  117 

FORMAT!/'  17.  ENTER  THE 
/7X'!THREE  REAL  VALUES:) 

ACCEPT  202,OPER 
IF!ICONT.GT.O)GO  TO  300 
TYPE  118 

FORMAT!/'  18.  ENTER  THE  ANNUAL  SCHEDULED 

,'P£R  UNIT'/7X'(THREE  REAL  VALUES):  ',$) 

ACCEPT  202,HRSMNT 

IF! ICONT.GT.O)GO  TO  300 

TYPE  119 


RATIO,  IN  PERCENTAGE  TERMS 


ANNUAL 

',$) 


OPERATING  COST  PER  UNIT' 


MAINTENANCE  HOURS 


FORMAT!/'  19.  ENTER  THE  REPAIR  WAGE  RATE' 

/7X'!THREE  REAL  VALUES):  ',$) 

ACCEPT  202,WRATE 
IF!ICONT.GT.O)GO  TO  300 
TYPE  120 

FORMAT!/'  20.  ENTER  ANNUAL  PASSENGER  VOLUME  AT  STATION' 
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1 


21 

22 

1 

2 

25 

403 

1 


410 

1 

2 

24 

294 

1 

3 

4 

242 


2 41 


201 

202 

20J 

204 

667 

671 


672 


673 

1 

666 

1 

2 

3 

4 

5 

6 


/7X'(0NE  integer  VALUE):  ',5) 

ACCEPT  203,PASSYR 

1F( ICGNT.GT.0)GO  TO  300 

GO  TO  6t)7 

1C0NT3=C 

TYPE  22 

EQRMAT( /6X' ARE  ANY  Of  THE  GRAPHS  TO  BE  ALSO  PRINTED', 

' DIRECTLY'  ON  THE  T ERM  IN  AL- '/6X  ' ( ZERO  FOR  "NO"', 

' OP  THE  GRAPH  NUMBER  (1,2,3, 4 OR  5)  FOR  "YES"):  '$) 
ACCEPT  403,PRINT 

format  (I) 

IF  (PR  INT.  EQ.O.OR.PRI  NT.  EQ.  4.  OR  . PRI  NT.  ECl.  5.  OR.  PRI  NT.  EQ.  6 
. OR. PRINT. E^. 7. OR. PRINT. EQ. 8)  GO  TO  24 
1C0NT3=IC0NT3+1 
IF(IC0NT3.GT.3)  STOP  12 
TYPE  410 

FORMAT(5X'THE  PREVIOUS  RESPONSE  WAS  ',15, 

/6X'PLEASE  ANSWER  WITH  ONE  INTEGER  VALUE', 

'EITHER  0,  1,  2,  3,  4,  5:  '$) 

GO  TO  25 
TYPE  294 

FOR^^AT(6X'TliE  MODEL  PRODUCES  OUTPUT  FILES  SHOWING  MEAN  ' 
,'GUEUE  length  AND  DELAY, '/6X' AND  QUEUE  LENGTH  ' 

, 'DISTRIBUTICN  (AND  COSTS  IF  INCLUDED).' 

/oX'SHOULD  THESE  BE  DISPLAYED  ON  THE  TERMINAL  AS  WELL?' 
/6X'("Y"  OR  "N"):  '$) 

ACCEPT  302,ANSR 
IPk  1NT=-1 

1F( ANSR.EQ.YES)IPRINT=1 
IF( ANSR.EQ. N0)IPRINT=0 
IFCIPRINT.GL.OGO  TO  241 
IC01jT=ICDM  + l 
IF{ ICONT.GT.3)STOP  77 
type  303,ANSR 
GO  TO  242 
ICuNT3=0 

IFCICONT.GT.OGO  TO  300 
h OR^YATOl  ) 

F0PMAT(3F ) 

FORMAT(T) 

FOPMAT(F) 

TYPE  671 

F 0kMAT(6X'W0ULD  YOU  LIKE  TO  SEE  THE  CURRENT  INPUT  FILE?'/ 
6X'("Y"  OR  "N"):  '$) 

ICONT3=0 


ACCEPT  302,ANSR 
IF(ANSR.EQ.NO)GO  TO  674 
IF ( ANSR.EQ.YES)GO  TO  673 
1 CONT3=ICONT3+1 
IF(ICONT3.GT.3)STOP  16 
TYPE  303, ANSR 
GO  TO  672 

TYPE  666,  NM,  ARRATE, NARR, SRVBAR, MCBF,RPRBAR,TIMEND, 
DLYMAX,PRTIM, ALT,PROB,ISEED 

FURMAT(/,'  1.  NUMBER  OF  MACHINES  IS  ',  318, 

/,'  2.  MEAN  ARRIVAL  RATS  IS  ',  FlO.2, 

/,'  3.  NUMBER  OF  ARRIVALS  IN  GROUPS  IS  ',  18, 

/,'  4.  MEAN  SERVICE  TIME  IS  ',3F10.2, 

/,'  5.  mean  CYCLES  BETWEEN  FAILURES  IS  ',318, 

/,'  6.  MEAN  TIME  TO  REPAIR  IS  ',3F10.2, 

/,'  7.  TIME  TO  END  THE  SIMULATION  ',F10.2, 
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670 


669 

664 

665 
668 

674 


300 

301 


309 

302 


303 

30o 

304 

311 

312 
306 


675 

676 


7 

//' 

8. 

MAXIMUM  DELAY  TIME  IS 

',F10.2, 

8 

9. 

TIME  AT  WHICH  ARRIVAL 

RATE  CHANGES  ',F15.2, 

9 

10. 

NEW  MEAN  ARRIVAL  RATE 

',F10.2, 

A 

tr 

11. 

ARRIVAL  PROBABILITIES 

ARE  ',3F10.2, 

b 

//' 

12. 

RANDOM  NUMBER  SEED  IS 

',111/) 

IF(ICOST.EU.O)GO  TO  669 

TYPE  670, CUNIT, LIF £,DI SORT, SP ARES, OPER,HRSMNT,WRATE/PASSYR 
FORMATC/'  13.  CAPITAL  COST  PER  UNIT  IS  ',3F10.2 


r 

r 

r 

/' 

r 

r 


14. 

15. 

16. 

17. 

18. 

19. 

20. 


USEFUL  LIFE  IS  *,318 
DISCOUNT  RATE  IS  ',F10.2,'%' 

SPARES  RATIO  IS  ',3(F10.2, '%') 

UNIT  OPERATING  COST  IS  ',3F10.2 
HOURS  OF  SCHEDULED  MAINTENANCE  IS 
kAGE  RATE  IS  ',3F10.2 
STATION  PASSENGER  VOLUME  IS  ',19) 


',3F10.2 


ANNUAL 
ANNUAL 
REPAIR 
ANNUAL 

lGRAPH=PRINT-3 
IFdPRINT.EQ.OTYPE  664 
TFdPRINT.EQ.DTYPE  665 
F0RMAT(/6X'THE  OUTPUT  FILES 
, 'TERMINAL') 

FQRMAT(/6X'THE  OUTPUT  FILES  WILL  BE  DISPLAYED  ON  THE  TERMINAL') 
TYPE  668,  IGRAPH 
FORMAT(6X,'THE  GRAPH  DISPLAYED 
'TERMINAL  IS  ("0"  IF  NONE);  ' 

IF (CHANGE. EQ.1)CHANGE=0 
If  (CHANGE. EU.DGO  TO  888 
ICONT=l 


WILL  NOT  BE  DISPLAYED  ON  THE 


ON  THE 
',15) 


ICONT1=0 
TYPE  301 

FORMAT(//'  DO  YOU  WISH  TO  MAKE  ANY  CHANGES  ?', 

1 /,'  ("Y"  OR  "N"):  '$) 

ACCEPT  302,  ANSR 
FORMAT! Al) 

1F(aNSR.£Q. NO. AND. CHANGE. EG. l)GO  TO  667 

IF(ANSR.EQ.NC)  GO  TO  888 

IF ( ANSR. EQ. YES)  GO  TO  306 

ICO NT1=ICGNT1+1 

lFdCONTl.GT.3)  STOP  2 

TYPE  303,  ANSR 

FORMAT!/,'  THE  PREVIOUS  RESPONSE  WAS 

1 AI,"'.  PLEASE  ANSWER  "Y"  OR  "N":  '$) 

GO  TO  309 
ICONT2=C 
TYPE  304 

FaPMAT(6X'ENTER  THE  NUMBER  OF  THE  QUESTION  YOU  WISH  ' 

1 ,'TQ  CHANGE') 

IF(CHANGE.EQ.O)TYPE  311 

FORMAT(6X'(  USE  "21"  10  CHANGE  REPORTS  DISPLAYED  ON  ' 

1 , 'TERMINAL)-') 

TYPE  312 

FORMAT!'  (ONE  INTEGER  VALUE):  '$) 

CHANGE=1 
ACCEPT  203,1 

IF!1.LT.1.0R.I.GT.21)  GO  TO  305 

If (.NOT. ((I. GE. 13  .AND.  I.LT.21)  .AND.  ICOST. EQ. 0) ) GO  TO  678 

CHANGE=0 

TYPE  675 

FOkMAT('  COSTS  ARE  PRESENTLY  NOT  INCLUDED  IN  THIS  MODEL.' 

1 ,'  DO  YOU  WISH  TO  ADD  COST  DATA? '/7X'(  "Y"  OR  "K"):',$) 

ICONT3=0 
ACCEPT  302, ANSR 
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IF( ANSR.EQ. NO)GO  TO  300 
1F(ANSR.EQ. Yfc:S)GO  TO  677 
ICa^•T3  = TCON^3+l 
IF(ICONT3.GT.3)STOP  15 
TYPE  303,ANSR 
GO-  TO  67b 

677  TCONT=0 
1C0ST=1 
CPANGF=0 
GO  TO  13 

678  GO  TO  Cl/2,3,4,5, 6,7, 8,9, 10,ll,12,13,14,15<rl6,17,18,19^20,21)^I 

30b  1C0NT2=IC0NT2+1 

IF(ICONT2.GT.3)  STOP  3 
TYPE  310 

310  FQPMATC'  THE  PREVIOUS  RESPONSE  WAS  ',15, 

1 /,'  PLEASE  ANSWER  WITH  ONE  INTEGER  VALUE', 

2 'FROM  1 TO  12:  '$) 

GO  TO  308 

888  WRITE(9,666)NM, ARRATE,NARR,SRVBAR,MCBF,RPRBAR 

1 ,T1MEND,DLYMAX,PRT1M,ALT,PR0B,ISEED 

IF( ICOST.EQ.1)WR1TE(9,670)CUNIT,LIFE,DISCRT,SPARES,OPER 
1 , HRS.MNT,IvRATE,PASSYR 

RETURN 
END 

SUBROUTINE  ENDING 

CORMUN7S4/I5,IMULT(1000),IQUE(1000),IUNIT(1000),LOC2(1000) 
COMMON/ 35/ NM, ARP ATE, N ARP, SRVBAR,MCBF,RPRBAR,TIMEND,DLYMAX 

1 ,PRT1M,ALT,PR0B,IS£ED,LAST,IC0ST,IPRINT 

2 ,conit,life,discrt,spares,cper,hrsmnt,wrate,passyr 
dimension  NM(3),SPVBAR(3),MCBF(3),RPRBAP(3),PR0P(3),CUNIT(3) 

i , LIFE (3 ),0PEP(3),HRSMNT(3),WRATE(3) 

DOUBLE  PRECISION  IF NX 

INTEGER  ANSR,NO,YES,(jUE,QL,PRINT,PUNIT 

I^'ISGER  LSTAK(51),STAR,BLANK,QSTAR,QLSTAR(101),  UNIT 

DATA  YES,L0/1HY,1HN/ 

DATA  STAR/'*'/,  BLANK/'  '/,  QSTAR/'#'/ 

I^'r  2 = 0 


C WHEN  THE  PROGRAM  ENDS,  THE  STATEMENT  'STOP'  APPEARS  AT  THE 
: END  OF  THE  TERMINAL  PRINTOUT,  IF  THE  PROGRAM  ENDS 

: normally,  THE  'STOP'  WILL  BE  UNACCOMPANIED  BY  ANY 

C number.  however,  after  SOME  TYPES  OF  ABNORMAL  ENDINGS,  THE  'STOP' 
: STATEMENTS  WILL  HAVE  AN  OCTAL  NUMBER  AFTER  THEM. 

C THESE  NUMBERS  ARE  AS  FOLLOWS: 


c 

(STOP 

— 

NORMAL  END 

i-OF-JOB) 

r* 

STOP 

1 - 

MORE  THAN 
REGARDING 

THREE  BAD  RESPONSES  TO  THE  QUESTION 
ACCEPTING  INPUT  FROM  THE  TERMINAL 

n 

STOP 

2 - 

MORE  THAN 
REGARDING 

THREE  BAD  RESPONSES  TO  THE  QUESTION 
MAKING  CHANGES  TO  THE  INPUT  DATA 

n 

c 

STOP 

3 - 

MORE  THAN 
REGARDING 
CHANGED 

THREE  BAD  RESPONSES  TO  THE  QUESTION 
THE  NUMBER  OF  THE  QUESTION  TO  BE 

V, 

c 

c 

STOP 

4 - 

MORE  THAN 
REGARDING 
DATA  FILE 

THREE  BAD  RESPONSES  TO  QUESTION 
NEXT  RECORD  OF  CURRENT  DISK 
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O O 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C' 


2 


302 

301 
4 80 
520 
202 


201 

400 

401 


STOP 

5 

STOP 

6 

STOP 

7 

STOP 

10 

STOP 

11 

STOP 

12 

STOP 

13 

STOP 

14 

MORE  THAN  THREE  BAD  RESPONSES  TO  THE  QUESTION 
REGARDING  MAKING  ANOTHER  RUN 

THE  ARRAY  'BRK'  CONTAINING  2500  RANDOM 
NUMBERS  IS  BEING  EXCEEDED 

THE  ARRAY  'RPR'  CONTAINING  90  RANDOM  NUMBERS 
IS  BEING  EXCEEDED 

THE  NUMBER  OF  MACHINES  FOR  A SERVICE  UNIT 
(N.M(UNIT))  HAS  GONE  NEGATIVE 

MORE  THAN  THREE  BAD  RESPONSES  TO  THE  QUESTION 
REGARDING  DELETION  OF  DISK  STORAGE  FILES 

MORE  THAN  THREE  BAD  RESPONSES  TO  THE  QUESTION 
REGARDING  CHOICE  OF  OUTPUT  GRAPHING  ROUTINES 

MORE  THAN  THREE  BAD  RESPONSES  TO  THE  QUESTION 
REGARDING  ADDITIONAL  PRINTOUT  ON  THE  TERMINAL 

MORE  THAN  THREE  BAD  RESPONSES  TO  THE  QUESTION 
REGARDING  PRINT  NUMBER  REPORT  AND/OR  UNIT  NUMBER 


LAST=0 

DO  2 1=1,51 

LSTAkCl )=BLANK 

CONTINUE 

DO  3 1=1,101 

QLSTAR(I)=BLANK 

CONTINUE 

ICONT=l 

1CONTO=0 

GO  TO  33 

FQkMATUHI,'  PASSMGR','  at  ARRIVAL  SERVI CE', 3X 'DELA Y ', 
1 /,9X, 'UNIT', 4X,'TIME',5X,' TIME', 4X'TIME'/) 

FOkMAT(1XI5,2X,I4,2X,3F8.1,2X51A1) 

FGPVATdHl,'  TIME',2X'QUFUE',4X,'AT',/7X'LENGTH  UNIT'//) 
FnRMAT(lXI4,2(3X,I4),3X101Al) 

F0<^MAT(1H1,2X'TIME','  EVENT ',6 X 'N EXT ',7X'NEXT ',5X'NEXT  ' 

1 ,4X'QU£UE',3X'MACHINES  AT '/19X'ARRIVAL',2X'DEPARTURE' 

2 ,2X'REPAIR',2X'LENGTH',2X,'IN  SERVICE  UNIT'/) 

FORMATdXF  6,1,2XA5,3XF7.1,4XF7.1,2XF7.1,4XI4,5XI4,5X,I4) 
FORMATdHl,  ' PASS  UNIT  ',' ARRIVAL', 2X'QUEUE'/15X'TIME', 

1 4X, 'LENGTH'/) 

F QKMAT (1X14, IX, I4,3XF7. 1, 3X14, 3X101 Al) 


SEE  IF  ANY  OTHER  REPORTS  ARE  TO  BE  PRINTED  ON  THE  TERMINAL 


33  ICONT1=0 

IF( ICONT.EQ.O)TYPE  2000,PUNIT 
2000  f ORMATOX'UNIT  ',11,'  INACTIVE  OR  NON-EXISTENT') 

TYPE  10 

10  FORMATC/,'  DO  YOU  WISH  TO  SEE  ANY  OTHER  GRAPHS-',/, 

1 t)X,'(ANSWER  "Y"  OR  "N"):  ',$) 

12  ACCEPT  1,ANSR 

1 FQRmaT(AI) 

IF(ANSR.EQ.YES)  GO  TO  19 
IF(ANSR.EQ. NO)  GO  TO  9 
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lCQMTl=ICnNTl+l 
IF (ICGNTl.GT. 3)  STOP  13 
TYPE  308,ANSR 
GO  TO  12 

19  ICONT2=0 

IF(  ICONTO.EQ.OTYPE  901 

901  FOPMATC'  GRAPHS  CAN  BE  DISPLAYED  ON  THE  TERMINAL  IN  ' 

1 , 'THEIR  ENfiRETY'/'  OR  IN  A SHORT  FORM  CONSISTING  OF  EVERY' 

2 TENTH  ENTRY.  TO  REQUEST'/'  THE  SHORT  FORM,  ADD  10  TO  THE' 

3 ,'  GRAPH  NUMBER  CE.G.,  FOR  THE  SHORT  FORM'/'  OF  GRAPH  4,' 

4 ,'  TYPE  "14")'// 

5 ' TO  STOP  DISPLAY  OF  ANY  PRINTOUT,  TYPE  CONTROL-O'/ 

6 ' C"0",  NOT  "ZERO"),  THEN  RESPOND  "Y"  OR  "N"  TO  THE'/ 

7 ,'  QUESTION  "DO  YOU  WISH  TO  SEE  ANY  FURTHER  PRINTOUTS?"'/ 

8 ' (THIS  QUESTION  WILL  NOT  BE  DISPLAYED  ON  THE  TERMINAL)') 
ICGNT0=1 

ICONT=0 
TYPE  20 

20  FO»MAT(/,'  ENTER  THE  GRAPH  NUMBER  (1,2,3, 4 OR  5)' 

1 AND  THE  UNIT  NUMBER  (0=ALL,1,2  OR  3)-',/, 

2 ' (TWO  INTEGER  VALUES):  ',$) 

C 

C NOTE  THAT  THE  VARIABLE  'PRINT'  DOESN'T  GET  TRANSFERRED  BACK  TO 
C THE  MAIN  ROUTINE  SO  THERE  IS  NO  WORRY  ABOUT  LOSING  ORIGINAL 
C VALUES 
C 

255  ACCEPT  21, PRINT, PUMT 

21  FORMAT(2I) 

IGkAPH=PRINT 

PRlNT=PRINT+3 

TF(((PRINT.GE.4  -and.  PRINT. LE. 8)  .OR. 

A (PRINT. GE. 14  .AND.  P RINT.LE. 18 ) ) .AND. 

1 (PUNIT.EQ.O.QR.PUNIT.EQ.1 .0R.PUNIT.EQ.2.0R.PUNIT.EQ.3)) 

i G D To  2 3 


1C0N’»’2  = IC0NT2+1 
IF  (ICONT2.GT.3)  STOP  14 
TYPE  29,1GR APH,PUNIT 

29  FORMAT(/, ' THE  PREVIOUS  RESPONSE  FOR  PRINT  WAS  ',13, 

1 ' A^JD  UNIT  WAS  ',13,/,'  PLEASE  ANSWER  1,2, 3,4', 

2 ' OR  5 FOR  PRINT  (OR  11,12,13,14,  OR  15  FOR  THE  ' 

3 , 'SHORT-FORM  PRINTOUTS)  '/7X'AND  0,1,2, 3 FOR  UNIT:  ',$) 

GO  TO  255 


23  ISHORT=0 

IF(PRINT.GT.IO) ISHORT=l 
IF(PRINT.GT.10)PRI NT=P  R I NT- 1 0 
PRINT=PRINT-3 

GO  TO  (24,25,26,27,28)  PRINT 


: REPORT  1 ON  THE  TERMINAL 

24  LD=4+10*PUNIT 

REWIND  LD 
PEAD(LD,302) 

1302  FORMATUHl,  ' PASS  AT  ARRIVAL  SERVICE  DELAY  ',/, 

1 ' UNIT  TIME  TIME  TIME  '/) 

K=-l 

342  READ(LD,301,END=33)NPASS,UNIT,TIME,SERVTM,DELAY,LSTAR 
1F(  ICONT.EQ.OTYPE  1302 
ICUNT=1 

IFdSHORT.EQ.OGO  TO  241 
K=K+1 
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241 


IF(K/10*10.  ?JE.K)GO  TO  342 
DO  34  1=1,38 

TFCLSTARd)  .EQ.STAR)GO  TO  341 
34  CONTINUE 

1 = 38 

LSTAR(38)=STAR 

341  TYPE  1301,  NPASS,UNIT,TIME,SERVTM,DELAY,(LSTAR(IJ),IJ=1,I) 

1301  F0RMAT(1XI5,1XI2,3F8.1,1X,38A1) 

GO  TO  342 
C 

C REPORT  2 ON  THE  TERMINAL 
C 

25  IFtl5.EQ.0)G0  TO  33 

K=-l 

DO  30  1=1,15 

IF(.NOT,(PUNIT.EQ.O  .OR.  PUNIT. EQ. lUNIT ( I ) ) )GO  TO  30 

IF(1CONT.EQ.O)TYPE  1480 

ICQNT=1 

IFC1SHORT.EQ.O)GO  TO  251 
K=K+1 

IF(K/10*1G.NE.K)GO  TO  30 
251  IJDEX=MINO(LOC2(I),55) 

QLSTAR(IJDEX)=QSTAR 

TYPE  1520,1MULT(I),IQU£(I),TUNIT(I),(QLSTAR(II),II=1,IJDEX) 
GLSTAR(IJDEX)=BLANK 
30  CONTINUE 

GO  TO  33 

1400  FOPMATCIHI,'  TIME  QUEUE  AT',/, 7X, 'L£ NGTH  UNIT'/) 

1520  F0RMAT(1X,I5,1X,I4,2X,I2,2X,55A1) 

u 

C R;,?ORT  3 ON  THE  TERMINAL 
C 

2b  REWIND  6 
TYPE  903 

903  fOKMATC'  WARNING  — GRAPH  3 CAN  BE  VERY  LONG.  TO  STOP' 

1 ,'  DISPLAY,  TYPE  CONTROL-0') 

KF.AD(6,202) 

K=-l 

31  READ(6,2C1,END=33)  TIME, LABEL, AR,DP,RP,QL, NMACH,UNIT 

IF(. NOT. (PUNIT. EQ.O  .OR.  PUNIT. EQ. UNIT) )GO  TO  31 
IF( ICONT. EQ.G)TYPE  202 
ICDM'=1 

IF(ISHORT.EQ.O)GO  TO  261 
K=K  + 1 

IF(K/10*10.NE.K)GO  TO  31 

261  TYPE  201, TIME, LABEL, AP, DP, RP,QL,NMACH, UNIT 

GO  TO  31 

C REPORT  4 on  THE  TERMINAL 
C 

27  LD=7+10*PUNIT 

REWIND  LD 
K=-l 

READ(LD,400) 

32  PEA9(LO,401,END=33)NP2,UNIT,TMIN,QL,QLSTAR 

IF ( ICONT. EQ.0)TYPE  1400 
ICONT=l 

IF( ISHORT.EQ.O)GO  TO  271 
K=K+1 

IF(K/10*10.N£.K)GU  TO  32 
271  DO  37  1=1,49 
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IFCQLSTARd  ).EQ.QSTAK)GQ  TO  371 
37  CONTINUE 

1-49 

QLSTAR(49)=QSTAR 

371  TYPE  1401, NP2, UNIT, TMIN,QL,(QLSTAR(IJ)^IJ=1^I) 

GO  TO  32 
C 

C REPORT  5 ON  THE, TERMINAL 
C 

28  LD=8+10»PUNIT 

REWIND  LD 
K=-l 

READ(LD,490) 

490  FORMAT('  DEPARTURE', 13X'INITIAL  *’,3X'TOTAL'/ 

1  ,3X'TIME',5X'PASSENGER',4X'UNIT',5X'DELAY'/) 

291  READCLD,309,END=33)TMIN,NPS,JFUNIT,TDELAY,LSTAR 

309  FORMAT (1XF6, 1,5X15, 9X12, 5XF7.1,2X51A1) 

IF(  ICUNT.EQ.OTYPE  490 
ICGNT=1 

1F(ISHORT.EQ.O)GO  TO  281 
K=K+1 

IF(K/10*10.N£.K)GO  TO  291 
281  DO  292  1=1,30 

IF(LSTAR(I) .EQ. STAR)GO  TO  293 

292  CONTINUE 
1 = 30 

LST AR(30)=STAR 

293  TYPE  130b,TMIN,NPS,JFUNIT,TDELAY,(LSTAR(IJ),IJ=l,I) 

130  6 FDr:MAT(lXF6.1,5XI5,9XI2,5XF7.1,2X30Al) 

GO  TO  291 

14DU  FQRMATCIHI, ' PASS  UN  ARRIVAL  QUEUE', /7X'IT',3X, 'TIME  LENGTH'/) 
1401  FQ^MaT(1X,I5,1X,I2,F7.1,1X,I4,1X,49A1) 

9 iCONT2=0 

TYPE  4 

4 FORMATC/'  DO  YOU  WISH  TO  MAKE  ANOTHER  RUN?'/ 

1 7X'(T:ITS  WOULD  DELETE  ALL  REMAINING  CURRENT  OUTPUT  REPORTS)'/ 

2 6X'(AN3WER  "Y"  OR  "N"):  '$) 

ISi|R£  = IPPINT 

9,^  ACCEPT  l,ANSk 

1F(ANSR,EQ. YES)  GO  TO  5 
IF (ANSR.EQ.no)  go  to  88 
ICOMT2=TCONT2+l 
1F(ICQNT2.GT.3)  STOP  11 
TYPE  308,ANSR 

308  FORMATC/,'  THE  PREVIOUS  RESPONSE  WAS 

1 Al,"'.  PLEASE  ANSWER  "Y"  OR  "N";  '$) 

GO  TO  98 

5 1F(ISURE.EQ«1)GG  TO  51 
TYPE  294 

294  FORMAT(6X'IF  YOU  MAKE  ANOTHER  RUN,  YOU  WILL  HAVE  NO  ' 

1 , 'QUEUE  LENGTH  OR  '/6X'COST  OUTPUTS  FROM  ' 

2 ,'THIS  RUN  (APART  FROM  THE  GRAPHS) . '/6X 'DO  YOU  STILL  ' 

3 ,'WISH  TO  MAKE  ANOTHER  RU N? '/7X '( "Y"  OR  "N"):  '$) 

ISURE=1 

GO  TO  96 

51  CLOSE(UNTT=4,DISPOSE='DELETE') 

CL03E(UNIT=6,DISP0SE= 'DELETE') 

CLOSEC UNI T=7,DISPOSE= 'DELETE') 

CL  CJSE(UNIT=a,DlSPOSE= 'DELETE') 

CLOSEC UNI T=9,DISP0SE='DELETE') 

CLQSE(UNIT=14,DISPOSE='DELETE') 
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CLCjSEC U MI  T=  17, DISPOSL= 'DELETE') 

CLGSE(UNIT=18,DISP0SE= 'DELETE') 

CL GS E( U MI T=24,DISP0SE= 'DELETE') 

CLUSK(UNIT=27,DISPOSE= 'DELETE') 
CLCSE(UNIT=26,DISP03£='DELETE') 
CLOSt'(UNIT=34,DISPOSE='DELETE') 

CLGSE(U MI T=37,DISP0SE= 'DELETE') 

CL GSE(UNIT=38,DISPUSE= 'DELETE') 

RETURN 
83  LAST=1 

CLGSE(UNIT=4) 

CL03E(UNIT=6) 

CLQ3£(UNIT=7) 

CLU3E{UNIT=8) 

CLGSE(UNIT=9) 

CLGSE(UNIT=14) 

CLC3E(UNIT=17) 

CLQS£(UNIT=18) 

CLGSE(UNIT=24) 

CL0SL'(UNIT=27) 

CL0SE(UNIT=28) 

CL0SE(UNIT=34) 

CLGSE(UNIT=37) 

CL0SE(UNIT=3fa) 

ICGNT1=0 

1F(I5.EQ.0)GQ  TO  40 

nP£N(UNIT=5,KTLE='FOR05.nAT',DEVICE='DSKS'/ACCESS='SEQINOUT') 
UPEN(UNIT=15,FILE='U15.DAT',DEVICE='DSKS'/ACCESS='SEQIN0UT') 
03FN(UNIT=25,FILE='U25.DAT',DEVICE='DSKS',ACCESS='SEQIN0UT') 
0?LN(UNIT=35,FIL£='U35.D AT', DEVICES 'DSKS',ACCESS='SEQINOUT') 
DG  90  L 0=5, 35,10 

90  ITE(LD,46C) 

DO  <^1  1 = 1,15 

Lj  = 5+10*IUhiIT(l) 

GLSTAR(LOC2(I))=STAR 

Ft;  ITE(LD,52C)  TMULT(I),IQUE(I),IUNIT(I  ), 

1 (QLSTAP(II),I1=1,LUC2(I)) 

TE(  5,5  20)1MULT(  I)  , IQU E C I ), lU N IT ( I ), 

1 (qLSTAR{II),II=l,LGC2(I)) 

GLSTAR(L0C2(I))=BLAHK 

91  CONTINUE 
CL03E(UNIT=5) 

CLG3r:(uMIT=15) 

CLOS£(UNIT=25) 

CLOSE(UNIT=35) 

40  irrjNTl=0 
TYPE  940 

940  f 0RMAT(6J('DG  YOU  WISH  TO  SAVE  THE  CURRENT  INPUT  DATA?'/ 

1 6X'(ANSWER  "Y"  OR  "N"):  ',$) 

41  ACCEPT  1,ANSR 
IF(ANSR.EQ. NO)RETUKN 
IF(ANSR.EQ.YES)GO  TO  42 
IC0NT1=ICQNT1+1 
IFdCONTl.GT.  3)STOP  15 
TYPE  308,ANSR 

GO  TU  41 

42  TYPE  942 

942  FQRYAT(6X'ENTER  FILE  NAME  (DP  TO  10  CHARACTERS):  ',$) 

ACCEPT  943, IFNX 

943  FQhMAT(AlO) 

0P£N(UNIT=56,FILE=1FNX,DEVICE='DSKS',ACCESS='SEQIN0UT') 
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WPITE(56,9  44)fJM,ARKATE,NARR,SRVBAR,  MCBF,RPRBAR,TIMEND 
1 ,DLYMAX,PRTIM,ALT/PROB,ISEED 

944  F0kMAT(lX3I3/lXFb.2/lXI3/lX3F8.2/lX3I6/lX3F8.2/lXF8.1 
1 /lXF8,l/lXFfl.l/lXF8.2/lX3F6.2/lXIll) 

IF  ( I COST.  EQ.l  )»JRITE(  56,  945  )CUNIT, LIFE, DISCRT,  SPARES,  OPER 
1 , HRSMNT, WRATE,PASSYR 

945  FORM AT (1X3F 9. 1/1X31 4/1 XF4. 1/1 X3F5. 2/1 X3F9.1/1X3F5.1 

1 /1X3F6. 2/1X112) 

CL0Sb:(UNIT=56) 

RETURN 

EED 

SUBROUTINE  CSTMOD(CUNIT,LIFE,DISCRT,SPARES,OPER, HRSMNT, WRATE 
1 ,?ASSYR,NM,MCBF,RPRBAR,PROB,CQST) 

DIMENSION  COST(3,6),CUNIT(3),LIFE(3),SPARES(3),OPER(3),HRSMNT(3) 
1 , ^RATE( 3),NM(3),MCBF(3 ),RPRBAR(3 ),PROB(3) 

INTEGER  PASSYR 
RATE=DISCRT/100. 

RFRAC=0. 

DO  1 1-1,3 
PFRAC=PFRAC+PROB(I) 

IF{LTFL(I).EQ.C)GO  TO  1 
S?KPCT=SPARES(I)/100. 

CAPINV=NM(I )*CUNIT(I) 

SPR INV  = C API  N'V* SPRPCT 
PV=1. 

IF  (RATE.LT..0DGO  TO  2 

PV=R AT£/(l.-( l./(l.+RATE)**LIFE(I))) 

2 CnST(I,l)=CAPINV*PV 
C03T(I,2)=SPRINV*PV 
COST(I,3)=OPER(I)*NM(I) 

COST(I,4)=HRSMNT(I)*taRATE(I)*NM(I) 

CrjST(I,5)=*IKAT£(I)*RPR6AR(I)/3600.*(PASSYR*PFRAC/MCBF(I)  ) 
CnST(I,6)=C. 

DO  3 J=l,5 

3 C0ST(I,6)=C0ST(I,6)+CQST(I,J) 

i corriMUE 

RETURN 

T?  »i’'. 

1,.  ■ L_/ 
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Main  Program 
MAIN. 

Subroutines 

QMATRX 

RMATRX 

CMTRX1 

CMTRX2 

PIVECT 

MATMLT 

VECMLT 

BEGIN 

INPUT 

ENDING 

CSTMOD 

MINV 

GELG 
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DIMENSION  Q( 0: 2,0 :2),R( 0:2,0 ;2)/LAMDA(0:2),PIX(0: 2)^ XX( 0:500) 

1 ,TEMP(0:2,0:2),ZERO(0:2,0;2),C1(0:2,0:2),C2(0:2,0:2),SUM(0:2) 

2 ,TEMP1( 0:2,0: 2 ),R1INV(0:2,0;2),LWV(3),M«V(3),A(0:2),X(0: 500,0:2) 

3 ,U1(0:2),U2(0:2),PI(0:2),XC(0:500,0 :2),QBARJ(0:5) 

DIMENSION  CC1(0:2,0:2),CC2(0;2,0:2),F1(0:2),F2(0:2),H(0;2) 

1 ,RMU(0:2),R2INV(0:2,0:.2),TEMPa‘(0:2),TEMPB(0:2),TEMPC(0:2) 

2 ,TEMPZ(0:2,0:2),UU(0:2,0:2),UU1(0:2,0:2),UU1MO(0:2) 

3 ,Y(0:2,0:2),Y1(0:2,0:2) 

DIMENSION  COSTC6) 

COMMON  IFILNM,FRDHWH,N,ARR,SERVTM,MCBF,MTTR 
1 ,CUNIT,LIFE,DISCRT,SPARES,OP£R,HRSMNT,KRATE,PASSYR,ICOST 
DOUBLE  PRECISION  IFTLNN 
REAL  LAMDA,MCBF,MTTP. 

INTEGER  C,FROMWH,PASSYR 
DATA  TOL,XTOL/.O0001,.00O5/ 

DATA  LIMIT,.MAXIT/500, 50000/ 

FROMWH=0 
CALL  BEGIN 

201  FORMAT(1XI3,4F6.0) 

202  F0RMAT(1X2F9. 6,215) 

34  WR1TE(6,201)N,ARR,SERVTM,MCBF,MTTP 

WRITE(6,202)XTOL,TOL,LIMTT,MAXIT 
SERV=3600./3ERVTM 
C=N 

SIGMA=(3600./MTTR)/SEKV 

THETA=1-/MCBF 

DO  1 1=0, N 

1 {.AMDA(I)  = ARR/SERV 

«RI TEC  6, 1 10 ) SIGMA, THETA, LAMDA 
110  F0RMAT(1X2F9.5,5X(F9.5)) 

DO  31  IC=1,C 
JC=IC 

CALL  P1VECT(N,0C,SIGMA,THETA,PI) 

RH0NUM=0. 

RHODEM=0. 

DC  32  J=0,N 

hHaNllM=RHnNUM+PI(J)*LAMDA(  J) 

32  PHODE'^=RHODEM  + PI  (J)*J 

k h 0=R  H 0 N n M/  RK  0 D EM 

31  »^RITE(6,108)IC,RHQ,(PI(J),J=0,N) 

109  F0RMATCI5,F7. 4,5X7F7.4) 

y^RITE(4,301)N,ARK,SERV,MCBF,MTTR 
301  FORMATC 'IFARE  COLLECTION  QUEUING  ANALYSIS'// 

* ' NUMBER  OF  UNITS  = ',12/'  ARRIVAL  RATE  = ' 

1 ,F6.0, ' /HR'/'  SERVICE  RATE  = ',F6,0,'  /HR'/'  MCBF  = ' 

2 ,F6.0/'  MTTR  = ',F6.0,'  SEC') 

WR1TE(7,302) 

DO  351  J=0,N 

351  lfcRITE(7,303)J,PI(J) 

kP.ITE(7,304)RHO 

3C4  FORMATC//'  THE  TRAFFIC  INTENSITY  IS  ',F7.4) 

IF(RHO.GE.1.0)GO  TO  30 

CALL  QMATRX(Q,N,C,SIGMA,THETA) 

DO  2 1=0, N 

2 WR1TE(6,101){Q(I,J),J=0,N) 

101  F0HMAT((1X7F7.4)) 

CALL  RMATRX(R,LAMDA, N,T0L,MAXIT,Q,H,F1,F2,UU,UU1,UU1MU, RMU) 
WRITEC6,102) 

102  FORMAT('O') 

DO  3 1=0, N 

3 V;RIT£(6,101)(R(I,J),J  = 0,N) 
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5 
4 

6 

7 

26 

8 

12 

10 

9 

27 

99 

999 

28 

11 

13 

103 

98 

29 


DO  4 1=0, N 
DO  ? J=0,N 

TEMP(I,J)=R(I,J)*J+Q(I,J) 

TEMP(I,I)=TEMF( I, I)-l-L AMD A( I ) 

CALL  MATMLTCTEMP, N,N,R,N,ZEKO) 

DC  6 1=0, N 

ZERO(I,I)=ZERO( I,I)^LAMDA(T) 
h'PlTE(6,102) 

DO  7 1=0, N 

WR1TE(6,101)(ZERO( 1,J),J=0,N) 

CALL  CMTRX1CC1,C2,H,U,R,SUM,LAMDA,Y,Y1) 

WP1TE(6,102) 

DO  26  1=0, N 

»vPlTE(6,101)(Cl(I,J),J  = 0,N) 

WR1TE(6,102 ) 

WK17E( 6,101 )SUM 
DO  3 1=0, N 
C2(I,0)=C. 

DO  6 J=0,M 
TEMP(I,J)=0. 

IF(J.EQ.I  .OR.  J.EQ.1+1)TEMP(I,J)=Q(I,J) 

IFCJ.EQ.l  .AND.  I.GT.0)TEMP(I,J)=TEMP(I,J)+Q(I,I-1) 
IFd.EQ.  J)TEMP(1,J)  = TEMP(I,J)-LAMDA(I) 

CALL  M ATMLTCC 1,N,N, TEMP, N, TEMPI) 

DO  12  1=0, N 
DO  12  J=0,N 

TEMPCI,J)=C2(1,J)+TEMP1(I,J) 

DO  9 1=0, N 
DO  10  J=0,N 
R1INV(1,J  )=-R(I,J) 

R1 1NV(1, 1 )=P1INV(I,I)+1. 

S^PITE(  6,102) 

DO  27  !=0,N 

« 1 TE ( 6 , 1 01 ) ( R1 1 NV  ( I,  J ) , J = 0 , N ) 

CALL  MINV(R11NV,N+1,D,LWV,MWV) 
wRlTE(6,102) 

DO  99  1 = 0,  :j 

ftPlTE(6,999)(kllNV(I,J),J=C,N) 

F0PVAT(1X7E11 .4) 
rtPTTE(6,102) 

DO  28  1=0, N 

WR1TE(6,101)(R11NV(I,J),J=0,N) 

DO  11  1=0, N 
TEMP(1,N)=SUM(I ) 

A(1)=0. 

DO  11  J=0,N 

TEMP(I,N)=TEMP(I,N)+R11NY(I,J) 

DO  13  1=0, N 
DO  13  J=0,N 
TEMP1(1,J)=TEMP(J,I) 

A(N)=1. 

EPS=5.E-4 

k»RlT£(6,102) 

F0PMAT(1X7F9. 4) 
kRlT£(6,lC2) 

DO  98  I=0,N 

WRITE(6,999)(TEMP(I,J),J=0,N) 

ksRlT£(6,102) 

DO  29  1=0, N 

WPITE(6,103)(TEMP(1,J),J=0,N) 

CALL  GELGCA, TEMPI, N-*-l,l, EPS, lER) 
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DO  14  1=0, N 
X(N-1,I)=ACI) 

WRITE(6,1C2) 

WPITE(6,101)(X(N-1,I),I=0,N) 

call  CHTRX2 (N, Q, R, R11NV,X,U1,U2,LAMDA, limit, CC1,CC2,Y,Y1 
1 , TEMR2, TEMP  A, TEMPB, TEMPO, R2INV) 

DO  15  JJ=2,N 

15  WRIT£(6,101)(X(N-JJ,I),I=0,N) 

WRIT£(6,102) 

WRITE(6,104)(U1(I),I=0,N) 

kRITE(6,104)(U2(I),I=0,N) 

104  F0RMAT(1X7F9.41 
XTUT=0. 

DO  16  J=0,N 
PIX(J)=0. 

DO  16  JJ=0, N-1 
PIX(J)=PIX(J)+X(JJ,J) 

16  XTOT=XTOT+X(JJ, J) 

JJ=N-1 

17  IFd-XTOT.LT.XTODGO  TC  19 

IF( JJ.GE.LIMIT)GO  TO  18 
JJ=JJ+1 

DO  20  J=0,N 
X(JJ,J)=C. 

DO  21  1=0, N 

21  X(JJ,J)=X(JJ,J)+X(JJ-1,I)*R(I,J) 

PIX(J)=P1X(J)+X(JJ,J) 

20  XTOT=XTOT+X{JJ,J) 

GO  TO  17 

16  WR1TE'(7,105)LIMIT 

105  FOk.MATC'  X DOES  NOT  5UM  10  ONE  AFTER",  15, ' QUEUE  LENGTH  POINTS') 

19  PRITE(7,305 )JJ,XTOT 

3C5  F0RMAT(///1X14,'  POINTS  OF  THE  QUEUE  LENGTH  DISTRIBUTION' 

1 /'  HAVE  BEEN  COMPUTED,  FOR  A TOTAL  PROBABILITY  OF  ',F7.4) 

106  F0RMAT(1XI4,3XF7.4) 

QLB  AR  = 0 • 

GLSTDV=0. 

DO  36  1=0, N 

OLE AR=QLPAR+U1(I) 

IFCPIX(i) .GE.XTQL)QLSTDV=QLSTDV-fU2(I)-(Ul(I)**2/PIX(I)) 

36  CONTINUE 

QLSTDV=SQRT(QLSTDV) 

WRITE! 4, 306)QLBAR,QLSTDV 

306  FORMAT!/'  THE  MEAN  QUEUE  LENGTH  IS  ',F7.2/ 

1 SX'wITH  A STANDARD  DEVIATION  OF  ',F7.2) 

DLYBAR=QLBAR/ARR*3600. 

kRITE!4,310)DLYBAR 

310  FORMAT!/'  THE  MEAN  DELAY  IS  ',F7.2,'  SEC.') 

WRIT£!7,101)!PIX!I),I=0,N) 

WRITE!7,102) 

WRITE! 4,302) 

302  FORMAT!/'OPROBABIL1TY  OF  NUMBER  OF  UNITS  AVAILABLE '/SX'NUMBER' 

1 ,'  OF  UNITS  AVAILABLE', lOX'PROBABILITY') 

DO  35  J=0,N 

35  WRITE!4,303)J,PIX!J) 

303  F0PMAT!13XI2,25XF7.4) 

WRITE!4,102) 

RHONUM=0. 

RHGDEN=0. 

DO  37  J=0,N 

RHONUM=RHONUM+PIX!J)*LAMDA(J) 
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37  RHCJDf;N  = RMODKN+PIX(  J)*J 

kho=rhonum/rhoden 

WR1TS( 4,304)RHQ 
JMAX=MIN0(9,JJ) 
fcRlTE(4,307)( J,J=0,JMAX) 

307  FORVAT(//20X'  QUEUE  LENGTH  — PROBABILITY  DENSITY '//4X1 0 (5X12 ) ) 

vJRlTEC  4,311) 

311  FORMAT(' 

DO  22  1=0, JJ 
XX(T )=0. 

DO  23  J=0,N 
XC(I,J)=0. 

IF(PIX(J),LT.XTOL)GO  TO  23 
XC(I,J)=X(I,J)/PIX(J) 

XX(I)=XXCI)+X(1,J) 

23  CONTINUE 

IF((1+1)/10*10.NE.(I+1))GO  TO  22 
lM9=I-9 

V.R1TE(4,112)IM9,  (XX(J  ),J  = IM9,I ) 

112  FORMAT (3X 13, 10(1X2PF5.1,'%')) 

30B  F0KMAT(1XI5,7XF6.4) 

22  wRlTL(7,107  )I,XX(I), (X(I,J),J=0,N) 

IP1=IM9+10 

IF(JJ-IM9.GT,9)WRTTE(4,112)IP1, (XX(J),J=IP1,JJ) 

107  F0HMAT(1XI4,4X2PF6.2,'%',5X7(1X2PF6.2,'%')) 

WRITE(4,3C9)( J,J=0,JMAX) 

309  F0RRAT(//21X'QUEUE  LENGTH  — CUMULATIVE  PROBABILITY '//4X10(5XI2 ) ) 
WR1TE(4,311) 

WR1TE(7,102) 

NULL=0 

VJRlTE(7,107)i^ULL,XX(0),(XC(0,J),J=0,N) 

QBAR=1,-XX(0) 
uO  39'J=0,N 

39  QP  AP,J(  J)  = 1.-XC(0,J) 

DO  24  1=1, JJ 
DO  25  J=0,N 

Acn,  j)=xc(  i,j)+xc(  i-i,j) 

25  y5APJ(J)=QBARJ(J)+Cl.-XC(I,J)) 

XX(I  ) = XX(I  )+XX(I-l) 

QBAS  = QBAR  + ( l.O-XX(D) 

IF((1+1)/10*10.NE.(I+1))GO  TO  24 
lM9=I-9 

^RiT£(4,112)lM9,(XX(J),J=IM9,I) 

2 4 ksRITE(7,lC7  )I,XX{I),(XC(I,J),J=0,N) 

IP1=IM9+10 

IF( JJ-IM9.GT,9)WRITE(4,112)IP1,(XX(J),J=IP1,JJ) 

^RIT£(7,4C1 )QBAR,QBARJ 

401  FORMATC'  THE  COMPUTED  MEAN  QUEUE  LENGTH  IS',F7.2/5X7F8. 2 ) 

GO  TO  33 

30  WRITE(6,109) 

V;RITE(4,109) 

109  FORMAT('  THF  QUEUE  IS  NOT  STABLE') 

WRITE(4,304)kHO 

33  IF( ICOST.EQ,0)GO  TO  38 

CALL  CSTMODC  CUMT,LIFE,DISCRT,  SPARES, OP  ER,HRSMNT,WRATE 
1 ,°ASSYR,N,MCBF,MTTR,C0ST) 

WRITE(4,111)C0ST 

111  F0kMAT(/15X'—  COSTS  — '/ 

1 /7X'ANNUALIZED  CAPITAL  COST ',T35, '$ ',F9. 2 

2 y7X'ANNUALIZED  SPARES  COST', T35, '$ ',F 9. 2 

3 /7X'UPERATING  COST',T35, '$ ',F9. 2 
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4 /7X'SCHEDULED  MAINTENANCE  COST T35, F9. 2 

5 /7X'COHRECTIVE  REPAIR  COST T35, '$ ',F9. 2 

6 //7X'TOTAL  COST',T34,'$'/F10.2) 

38  CALL  ENDING 

1F(FR0M«,H.LT.2)G0  TO  34 

STOP 

END 

SU3R0UTINE  QM ATPX( G, N, C, SI GMA^ THETA) 

DIMENSION  Q(0:N,0:N) 

INTEGER  C 
DO  1 1=0, N-1 

G(I,I+1)=MIN0(C,N-I)»SIGMA 

1 G(I+1,I )=(I+1)*THETA 
DO  2 1 = 1, N-1 

2 Q(I,I)  = -(Q(I,I-1)+Q(I,I  + D) 

Q(0,0)=-G(0,1) 

Q(N,N)=-Q(N,N-1) 

RETURN 

END 

SUPRGUTINE  RMATRX(R,LAMDA,N,T0L,MAXIT,Q,H,F1,F2,U,U1,U1MU,KMU) 
DIMENSION  r(0:n,0:n),lamda(0:n),h(0;N),f1(0;n ),F2(0:n) 

1 ,u(0:n,0:n),ui(0:n,0:n),uimu(0:n),rmu(0:n),q(0:n,0;n) 

REAL  LAMDA 
ITS  = 0 

DO  16  1=0, N 
DO  16  J=0,N 
16  R(I,J)=0. 

DO  1 1=0, N 

1 H(1)=LAMDA(I)+I-Q(I,I) 

DO  2 1=1, N 

F1(I)=Q(I,1-1)/H(I-1) 

2 F2(I)=Q(I-1,I)/H(I) 

3 DO  4 1=0, N 
DO  4 J=0,N 

4 U(I,J)=R(I,J)*J/H(J) 

DO  5 1=1, w 

U(T,I-1)=U(I,T-1)+F1(I) 

5 U(I-i,I )=U(I-1, I)+F2(I ) 

DO  6 1=0, N 

DO  6 J=0,N 
U1(I,J)=0. 

DO  6 K=0,N 

6 U1(T,J)=U1(I,J)+R(1,K)*U(K,J) 

DO  7 1=0, N 

7 U1(I,I)=U1(I,I)+LAMDA(I)/H(I) 

ITS=ITS+1 

ERR=0. 

DO  8 1=0, N 
DO  8 J=0,N 

8 ERR=AyAXl(ERR,ABS(Ul(I,J)-R(I,J))) 

IF(ERR.LE.T0L)GQ  TO  11 
IF(ITS,LT.MAXIT)GO  TO  9 
WRITE(6,101)MAXIT 

101  FORMAT('  NO  CONVERGENCE  AFTER', 14,'  ITERATIONS.') 

GO  TO  11 

9 DO  10  1=0, N 
DO  10  J=0,N 

10  R(I,J)=U1 (I,J) 

GO  TO  3 

11  kRITE(6,102)ITS 

102  FQRMAT( 'O', 10X15) 
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GO  TO  14 
DO  20  1=0, N 

20  WRITE(6,103)(R( I,J),J=0,N) 

103  FQRMAK  1X7F7,  4) 

DO  12  1=0, N 
RMU(I)=0. 

UlMUd  )=0. 

DO  12  J=1,N 

RMU(I)=PMU(I)+R(I,J)*J 

12  UlMUd  )=U1MU(I)  + U1CI,J)*J 
ISSAME=0 

DO  13  1=0, N 

lF(UlMUd),EQ.RMUd))GO  TO  13 
ISSAME=1 

TEMP=(LAMDA(1)-U1MU(I)  )/(RMU(I)-UlMU(D) 

DO  13  J=0,N 

R(I,J)  = Uld,J)  + TEMP*(Rd,J)-Ul(I,J)) 

13  CONTINUE 

IFdSSAME.EQ.  0)GO  TO  14 

RETURN 

14  DO  15  1=0, N 
DO  15  J=0,N 

15  R(I,J)=U1(I,J) 

RETURN 

END 

SUBROUTINE  CMTRXl ( Cl, C2, N, Q, R, SUM,L AMDA, Y, Y1 ) 

DIMENSION  Cl(0;N,0;N),C2(0:N,0:N),Q(OrN,0;N),R(0;N,0:N),SUM(0;N) 
1 ,Y1(0:n,0:n),Y(0:N,0:N) 

REAL  LA''DA(0:N) 

DO  1 1=0, N 
SUMd)  = 0. 

DO  1 J=0,N 
C2( 1,J)=C. 

QIJ=U(I,0) 

IFCI.EQ.N  .AND.  J. EQ. N-1  )QI J=Q I J-THET A 
IFd.EQ.N  .AND.  J . EQ.  N ) QT  J=QI  J+THET  A 
Cl  d,J)=-(QIJ  + R(I,vj)*J)/LAMDA(J) 

IF(T.NE.J)GU  TO  1 
C2(I,J)=1. 

Cld,J)  = Cl(I,J)>(l.  + AMINO(N-l,J)/LAMDA(J)) 

1 SUMd)=SUMd)-^Cl(I,J) 

1F(N.EC4.2  )RETURN 
JJ=3 

2 kRITE(6,101)JJ 

101  F0RMAT(1XI3) 

DO  6 1=0, N 

6 WRITE(6,102)(C1(I,J),J=0,N) 

102  FORMAT( (1X6F7.4)) 

DO  7 1=0, N 

DO  7 K=0,N 
Y( I,K)=0. 

DO  7 J=0,N 
QJK=Q(J,K) 

IF( J.LE.N-JJ+DGO  TO  7 

IF(J.EQ.K)QJK=QJK+(J-(N-JJ+1))*THETA 

1F(J.EQ.K+1)QJK=QJK-(J-(N-JJ+1))*THETA 

7 Y(I,K)=Y(I,K)-t-Cld,J)*QJK 
DO  3 1=0, N 

DO  3 J=0,N 

3 Y(I,J)  = Yd,J)-ClCI,J)*(LAMDA(J)-»-AMINO(N-JJ-*-l,J)) 

DO  4 1=0, N 
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DO  4 J=0,N 

Y1(T,J)=(-Y(I,J)-C2(1,J)*AKIN0(N-JJ+2,J))/LAMDA(J) 
SUM(1)=SUM( I)+Y1(I,J) 

DO  5 I=C,N 
DO  5 J=0,N 
C2(T,J)=C1(I,J) 

Cl(l,J)=Yl(I,J) 

if(jj.le:.n)gq  to  2 

RETURN 

END 

SUBROUTINE  CMTRX2 ( K,Q, R,R1 IN V,X, Ul, U2,L AMDA,LIMIT/Cl,C2x Y1 
r TEMPI, TEMP  A, TEMPB, TEMPO, R2 1 NV) 

DIMENSION  Q(0;n,0;n),R(0:n,0;N),riinv(0:n,0:n),u1(0:n) 

,U2(0: n),x(0:limit,0:N),lamda(0:n) 
dimension  CKO:  n,0;n),y(0:  n,0:n),y1(0;n,0:n) 
,temp1(0:n,0:N),tempa(0;N),tempb(0;n),tempc(0:n) 
dimension  r2Inv(0:n,0:n),c2(0:n,0:N) 

REAL  LAMDA 
DO  1 1=0, N 
DO  1 J=0,N 
C2(I,J)=0. 

QIJ=Q(I,J) 

Ifd.EQ.N  .AND.  J.EQ.N-1)QTJ=QIJ-THETA 
IFCI.EQ.N  .AND.  J. £Q. N )QI J=QIJ+THETA 
Cl(I,J)=-(aiJ+R(I,J)*J)/LAMDA(J) 

IF(I.NE.J)GO  TO  1 
C2CI,J)=1. 

C1(I,J)=C1(I,J)+(1.+AMIN0(N-1,J)/LAMDA(J)) 

CONTINUE 
DO  2 J=0,N 
X(N-2,J)=0. 

DO  3 1=0,N 

X(N-2,J)=X(N-2,J)+X(N-1,I)*C1(I,J) 

U1(J)=(N-2)*X(N-2,J) 

U2(J)=(N-2)’‘U1(J) 

IF(  n'.EQ.2)GQ  TO  9 
JJ=3 

DO  16  1=0,1. 

DO  16  K=0,N 
Y(I,K)=0. 

DO  16  J=0,N 
QJK=Q( J,K) 

IF(J.LE.N-JJ+1 )G0  TO  16 

IF(J.Ea.K)«JK=QJK+(J-(N-JJ*l))'*THETA 

1F(J.EQ.  K+1  )QJK=QJK-(J-(N-JJ-i-l))*THETA 

Y(1,K)=Y(I,K)+C1(I,J)*QJK 

DO  5 I=C,N 

DO  5 J=0, N 

Y(I,J)=Y(I,J)-C1(T,J)*(LAMDA(J)+AMIN0(N-JJ+1,J)) 

DO  6 J=0,N 
X(N-JJ,J)=0. 

DO  7 1=0, N 

Y1(I,J)=(“Y(I,J)-C2(I,J)*AMIN0(N-JJ+2,J))/LAMDA(J) 

X(N-JJ,J)=X(N-JJ,J)+X(N-1,T)*Y1(I,J) 

U1(J)=U1(J)+(N-JJ)*X(N-JJ,J) 

U2( J)=U2(J)+(N-JJ)*U1(J) 

DO  8 1=0, N 
DO  8 J=0,N 
C2( I,J)=C1( I,J) 

C1(I,J)=Y1(I,J) 


Uj  NJ 


JJ=JJ+1 

IF(JJ,LS. N)GO  TO  4 

9 CALL  MATMLT(R1INV,N,N,R1INV,N,R2INV) 

DO  10  1=0, N 

DO  11  J=0,N 

11  TEMP1(1,J)=-(N-2)*R(I,J) 

TEMPK I,I)=TLMP1(I,I)+(N-1) 

10  TEMPACI )=X(N-1,I) 

CALL  VECHLT(TEMPA,N,R2INV,N,TEMPB) 

CALL  VECMLTC TEMPB,N, TEMPI, N,TEMPA) 

DO  12  1=0, N 

12  U1(I)=U1(I)+TEMPA(I) 

DO  13  1=0, N 

DO  14  J=0,N 

14  TEMP1(I,J)=(N-2)**2*R(I,J) 

13  TEMPK 1,1) =TEMP1 (!,!)-( (2. *N-6. )»N+3.) 

CALL  VECMLT(TEMPB,M,R1INV,N,TEMPA) 

CALL  VECMLTC TEMPA, TEMPI, N,TEMPB) 

CALL  VECMLTCTEMPb,N,R,N,TEMPC) 

DO  15  1=0, N 

15  U2(I)  = U2(I)+(TEMPC(I)-KN-1)**2*TEMPA(I)) 

RETURN 

END 

SUBROUTINE  PI VECT( N,C, SIGMA, THETA, PI) 

DIMENSION  PICOZN) 

INTEGER  C 
PI(0)=1. 

SUM=1. 

CST=C*SIGMA/THETA 
DO  1 J=1,N-C+1 
P1(J)=PI( J-1)*CST/J 
1 SUM=SUM+PICJ) 

IF(C.L£,l)GO  TO  3 
DO  ? J=N’-C  + 2,N 

FI(J)=PI( J-l)*CST*(l.-(J-(N-C+l))/FLOAT(C))/J 
SnM=SUM+?I( J) 

DO  4 J=0,N 
PI(J)=PI(J)/SUM 
RETURN 
END 

SUBROUTINE  MATMLTC A, I, J, P, K, C ) 

DIMENSION  A(0; 1,0 :j),B(0:j,0 :K),C( 0:1,0 ;k) 

DO  1 11=0,1 
DO  1 KK=0,K 
C(II,KK)=0. 

DO  1 JJ=0,J 
1 C(II,KK)=C(I1,KK)+A(II,JJ)*B(JJ,KK) 

RETURN 

END 

SUBROUTINE  VECMLTC A, I,B,J,C) 

DIMENSION  AC0:I),RC0:I,0:J),CC0:J) 

DO  1 JJ=0,J 
CCJJ)=0. 

DO  1 11=0,1 

1 CCJJ)=CCJJ)-^ACII)*BCII/JJ) 

RETURN 

END 

SUBROUTINE  BEGIN 

COMMON  IFILNM,FROMWH,NM,ARRATE,SRVBAR,MCBF,VTTR 
1 ,CUNIT,LIFE,D ISCRT, SPARES, OP ER,HRSMNT,WRATE,P ASS YR,ICOST 
REAL  MCBF,MTTR 
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INTEGER  FROMWH,YES/NO,ANS,PASSYR 
DOUBLE  PRECISIUM  IfILNM 
DATA  Y£S,NQ/1HY,1HN/ 

OPEN(UNIT=4,FILE='FaR04.DAT',DEVICE='DSKS',ACCESS='SEQINQUT') 

OPEN(UNIT=6,FILE='FOR06.DAT',DEVICE='DSKS',ACCESS='SEQINOUT-) 

aPEN(UMT=7^FlLE='FOR07.DAT',DEVICE='DSKS',ACCESS='SEQINQUT'’) 

IF(FROMWH.EQ.l)GO  TO  7 

ICOST=0 

TYPE  101 

101  FORMAT!'  ACCEPT  INPUT  FROM  TERMINAL?'/6X'(  "Y"  OR  "N")  ;'$) 

1 ACCEPT  102, ANS 

102  FQRMAT(Al) 

IF( ANS.EQ.YES)GO  TO  5 
IF(ANS.EQ.NO)  GO  TO  2 
TYPE  103 

103  FORMAT! 'OTYPE  EITHER  Y 0^  N-  TRY  AGAIN.  '$) 

GO  TO  1 

2 TYPE  104 

104  FORMAT!/'  ENTER  INPUT  FILE  NAME  (UP  TO  10  CHARACTERS);  '$) 

ACCEPT  105, IFILNM 

105  FORMAT!A10) 

0P£N!UNIT=55,F1L£=IFILNM,DEVICE='DSKS',ACCESS='SEQIN') 

TYPE  106, IFILNM 

106  FORMAT!/'  THE  INPUT  FILE  IS;  ',A10) 

7 READ!55,*)NM,ARRATE,SRVBAR,MCBF,MTTR 

TYPE  107 

107  FORMAT!/'  APE  COSTS  INCLUDED  IN  THIS  RECORD?' 

1 /6X'  !"Y"  OP  "N");  '$) 

3 ACCEPT  102, ANS 

IF! ANS.EQ.NO)GO  TO  6 
IF(ANS. EQ.YES)GO  TO  4 
TYPE  103 
GO  TO  3 

4 F;£AD!55,*  )CUNI  T,  LIFE,  DISC  RT,  SPARES,  OP  ER,HRSMNT,WRATE,PASSYR 
1C0ST=1 

GO  TO  6 

5 CALL  INPUT 

o FRUMRH=i 

RETURN 

END 

SUBROUTINE  INPUT 

COMMON  IFTLNM,FROMteH,NM,ARRATE,SRVBAR,MCBF,MTTR 
1 ,CUMT,LIFE,DISCRT,SPARES,OPER,IiRSMNT,WRATE,PASSYR,ICOST 
INTEGER  FROMWH, CHANGE, YES, NO, ANSR,PASSYR 
REAL  MCBF,MTTR 
DOUBLE  PRECISION  IFILNM 
DATA  YES,N0/1HY,1HN/ 

CHANGE=0 

ICGNT=0 

IF!FROM«H.GT.O)  GO  TO  667 

1 TYPE  101 

101  FORMAT!/'  1.  ENTER  NUMBER  OF  MACHINES', 

1 /'  !ONE  integer  VALUE);  '$) 

ACCEPT  201,  NM 
IF!ICONT.GT.O)  GO  TO  300 

2 TYPE  102 

102  FORMAT!/'  2.  ENTER  MEAN  ARRIVAL  RATE  PER  HOUR-',/, 

1 ' !ONE  REAL  VALUE);  '$) 

ACCEPT  202,ARRATE 
IF!ICONT.GT.O)  GO  TO  300 

3 TYPE  103 
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103 

1 


4 

104 

1 


5 

105 

1 


401 

1 

132 


131 

6 

106 

1 


7 

107 

1 


8 

1 Ob 

1 


9 

109 

1 


10 

110 

1 


11 

111 

1 


12 

112 

1 


FORMATC/'  3.  ENTER  MEAN  SERVICE  TIME  IN  SECONDS-'^ 

(ONE  REAL  VALUE):  ' $) 

ACCEPT  202,  SRVBAP 
IF(TCONT.GT.O)  GO  TO  300 
TYPE  104 

FORMAT(/'  4.  ENTER  MEAN  CYCLES  BETWEEN  FAILURES-', 

/,'  (ONE  REAL  VALUE):  ',$) 

ACCEPT  202,MCBF 
IF(ICONT.GT.O)  GO  TO  300 
TYPE  105 

FQRMAT(/'  5.  ENTER  MEAN  TIME  TO  REPAIR  IN  SECONDS-', 

/,'  (ONE  REAL  VALUE):  ',$) 

ACCEPT  202,  MTTR 

IF( ICONT.GT.O)  GO  TO  300 

TYPE  401 

F0RMAT(//7X'D0  you  Wish  TO  INCLUDE  COSTS  IN  THIS  RUN?' 
/7X'("Y"  OR  "N"):  '$) 

ACCEPT  302,ANSk 
ICOST=0 

IF( ANSR.EO. NO)GO  TO  667 
IF(ANSR.EQ.YES)GO  TO  131 
IC0NT1=IC0NT1+1 
IF(ICONTl.GT.3)STOP  2 
TYPE  303,ANSR 
GO  TO  132 
ICOST=l 
TYPE  106 

FORMAK/'  6.  ENTER  THE  CAPITAL  COST  PER  UNIT' 

/7X'(ONE  REAL  VALUE):  ',$) 

ACCEPT  202,CUNIT 
IF(ICONT.GT.0)GO  TO  300 
TYPE  107 

FORMAT(/'  7.  ENTER  THE  USEFUL  LIFE  OF  THE  UNIT' 

/7X'(0NE  INTEGER  VALUE):  ',$) 

ACCEPT  201, LIFE 
IF (ICQNT.GT.O)GO  TO  300 

TYPE  loe 

FORMAT(/'  8.  tNTEF  THE  DISCOUNT  RATE,  IN  PERCENTAGE  TERMS' 
/7X'(ONE  REAL  VALUt):  ',$) 

ACCEPT  204,DISCRT 

IF(  ICONT.GT.O)GO  TO  300 

TYPE  109 

FORMAK/'  9.  ENTER  THE  SPARES  RATIO,  IN  PERCENTAGE  TERMS' 
/7X'(QNE  REAL  VALUE:)  ',$) 

ACCEPT  202, SPARES 
IF(ICONT.GT.O)GO  TO  300 
TYPE  110 

FOKMAT(/'  10.  ENTER  THE  ANNUAL  OPERATING  COST  PER  UNIT' 
/7X'(QNE  REAL  VALUE:)  ',$) 

ACCEPT  202,OPER 
IF(ICONT.GT.O)GO  TO  300 
TYPE  111 

FORMAK/'  11.  ENTER  THE  ANNUAL  SCHEDULED  MAINTENANCE  HOURS 
,'PER  UNIT'/7X'(ONE  REAL  VALUE):  ',$) 

ACCEPT  202,HRSMNT 
IF(ICONT.GT.O)GO  TO  300 
TYPE  112 

FORMAK/'  12.  ENTER  THE  REPAIR  WAGE  RATE' 

/7X'(0f'lE  REAL  VALUE):  ',$) 

ACCEPT  202,WRATE 
1F(IC0NT.GT.0)G0  TO  300 
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13 

113 


1 


201 

202 

203 

204 
667 

671 

1 

672 


673 

666 

1 

3 

4 

5 


670 

1 

2 

3 

4 

5 

6 
7 

674 


300 

301 

309 

302 


303 

1 


306 


304 


1 

2 

308 


TYPE  113 

FORMATC/'  13.  ENTER  ANNUAL  PASSENGER  VOLUME  AT  STATION' 
/7X'(0NE  integer  VALUE);  '^$) 

ACCEPT  203,PASSYR 
IF(ICONT.GT.O)GO  TO  300 
GO  TO  667 
FORMAT(I) 

FORMAT(F) 

FORMAT(I) 

FORMAT(F) 

TYPE  671 

F0RMAT(6X'W0ULD  YOU  LIKE  TO  SEE  THE  CURRENT  INPUT  FILE?'/ 
6X'("Y"  OR  "N");  '$) 

IC0NT3=0 
ACCEPT  302,AKSR 
IF( ANSR.EQ. NO)GO  TO  674 
IF(AiiSR.EQ.  YES)GO  TO  673 
ICOKT3=ICONT3+l 
IF(IC0NT3.GT.3)ST0P  16 
TYPE  303/ANSR 
GO  TO  672 

TYPE  666,  NM,  ARRATE, SRV3AR,MCBF,MTTR 
F0RMAT(/,'  1.  NUMBER  OF  MACHINES  IS  ',  18, 

/,'  2.  MEAN  ARRIVAL  RATE  IS  ',  FlO-2, 

/,'  3.  MEAN  SERVICE  TIME  IS  ',fl0.2, 

/,'  4.  MEAN  CYCLES  BETWEEN  FAILURES  IS  ',F9.0, 

/,'  5.  MEAN  TIME  TO  REPAIR  IS  ',F10.2) 

IF(ICOST.Eu.O)GO  TO  674 

TYPE  670, CUNIT,LIFE,DISCRT, SPARES, OPER,HRSMNT,WRATE,P ASS YR 
FOkMAT(/'  6.  CAPITAL  COST  PER  UNIT  IS  ',F10.2 
/'  7.  USEFUL  LIFE  IS  ',18 

/'  8.  DISCOUNT  RATE  IS  ',F10.2,'%' 

/'  9.  SPARES  RATIO  IS  ',F10.2,'%' 

/'  10.  ANNUAL  UNIT  OPERATING  COST  IS  ',F10.2 

/'  11.  ANNUAL  HOURS  OF  SCHEDULED  MAINTENANCE  IS  ',F10.2 

/'  12.  REPAIR  WAGE  RATE  IS  ',F10.2 

/'  13.  ANNUAL  STATION  PASSENGER  VOLUME  IS  ',19) 

IF( CHANGE. EQ.l) RETURN 

1C0NT=1 

1CONT1=0 

TYPE  301 

FORMAT(//'  DO  YOU  WISH  TO  MAKE  ANY  CHANGES  ?', 

/,'  ("V"  OR  "N");  '$) 

ACCEPT  302,  ANSP 
FORMAT(Al) 

IF(ANSR.EQ.  NO.AND.CHANGE.EQ.DGO  TO  667 

IF(ANSR.EQ. NO)  RETURN 

IF(ANSR.EQ.YES)  GO  TO  306 

IC0NT1  = IC0NT1•^1 

IF(ICQNT1.GT.3)  STOP  2 

TYPE  303,  ANSR 

FORMAT!/,'  THE  PREVIOUS  RESPONSE  WAS  '", 

Al,"'.  PLEASE  ANSWER  "Y"  OR  "N";  '$) 

GO  TO  309 
ICONT2=0 
CHANGE=1 
TYPE  304 

FORMAT!'  ENTER  THE  NUMBER  OF  THE  QUESTION  YOU  WISH  ' 

,'TO  CHANGE'/ 

' !ONE  INTEGER  VALUE);  '$ ) 

ACCEPT  203,1 
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1F(.N0T.(I.GE.1.ANL.I.LE.13))G0  TO  305 

IF(.N0T.((I.GE.6  .AND.  I.LE.13)  .AND.  ICOST. EQ. 0 ) ) GO  TO  678 

CHANGE=0 

TYPE  675 

675  FORMAK'  COSTS  ARE  PRESENTLY  NOT  INCLUDED  IN  THIS  MODEL.' 

1 DO  YOU  WISH  TO  ADD  COST  DATA? '/7A '( "Y"  OR  "N");'y$) 

ICDNT3=0 

676  ACCEPT  302,^ANSR 
IF(ANSR.EQ. NO)GO  TO  30C 
IF(ANSR.EO. YES)GO  TO  677 
ICONT3=ICONT3+l 
IF(ICONT3.GT.3)STOP  15 
TYPE  303,ANSR 

GO  TO  676 

677  ICONT=0 
ICOST=l 
CHANGF=0 
GO  TO  13 

678  GO  TO  (1,2,3,4,5,6,7,8,9,10,11,12,13)/! 

305  ICONT2=ICONT2+l 

1F(IC0NT2.GT. 3)  STOP  3 
TYPE  310 

310  FORHAT('  THE  PREVIOUS  RESPONSE  WAS  ',15, 

1 /,'  PLEASE  ANSWER  WITH  ONE  INTEGER  VALUE', 

2 'FROM  1 TO  12:  '$) 

GO  TO  308 

END 

SUBROUTINE  ENDING 

COMMON  IFILNM,FROMWH,NM,*RRATE,SRVBAR,MCBF,MTTR 
1 , CU NIT, LIFE, DISCRT, SPA RES, OP ER/HRSM NT, WRATE/PASSYR, ICOST 

REAL  MCBF,MTTR 

INTEGER  FROMWH,YES,NO,ANS,PASSYR 
DOUBLE  PRECISION  IFILNM 
DATA  VES,NQ/1HY,1HN/ 

FROMWH=2 
TYPE  101 

IC'l  FORN'.ATC 'ODD  YOU  WISH  TC  MAKE  ANY  FURTHER  RUNS?' 

1 /6X'  (ANSWER  "Y"  OR  "N")  ',$) 

1 ACCEPT  102, ANS 

102  FORvaT(AI) 

IF( ANS.EQ.YES)GO  TO  2 
IF( ANS. EU. NO)GO  TO  9 
TYPE  103 

103  F0RMAT( 'OTYPE  EITHER  Y OR  N.  TRY  AGAIN.  ',$) 

GO  TO  1 

2 FROMWH=l 
CLGSE(UN1T=4,DISP0SE='DELETE') 

CLOSE(UNIT=6,DISPOSE= 'DELETE') 

CLGSE(UNIT=7,DISPDSE='DELETE') 

TYPE  104 

104  FORMAT('OIS  THE  NEW  DATA  THE  SAME  AS  THE  CURRENT  DATA' 

1 ,'  (EXCEPT  FOR  SOME  CHANGES  )?'/6X'  (ANSWER  "Y"  OR  "N")  ',$) 

3 ACCEPT  102, ANS 
lF(ANS.Ea.YES)GO  TO  4 
IF(ANS. EQ.NO)GO  TO  5 
TYPE  103 
GO  TO  3 

4 CALL  INPUT 
RETURN 

5 TYPE  105, IFILNM 

105  FORMAT('OIS  THE  NEW  DATA  THE  NEXT  RECORD  OF: 
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,A10,  '?') 


oonoooonoooooorioo 


6 


7 


106 

8 


9 


2 


3 

4 

101 

1 


ACCEPT  102,ANS 

IF( ANS. EQ.YES)G0  TO  8 

1F(ANS.EU.NQ)G0  TO  7 

TYPE  103 

GO  TO  6 

FROMWH=0 

CL0SE(UNIT=55  ) 

TYPE  106 

FORMATCOTHE  NEW  DATA  IS  FROM  A NEW  INPUT  FILE.') 

CALL  BEGIN 

FROMWH=l 

RETURN 

CLOSE(UNIT=4) 

CL0SE(UNIT=6) 

CL0SE(UNIT=7) 

RETURN 

END 

SUBROUTINE  CSTMOD (CUNTT/LIFE,DISCRT, SPARES, OPER, HRSMNT^ WRATE 
1 ,i>ASSYR,NV,MCBF,MTTR,COST) 

DIMENSION  C0ST(6) 

INTEGER  PASSYR 
REAL  MCBF,MTTR 
RATE=DISCRT/100. 

IF(LIFE,EQ.O)GO  TO  4 
SPRPCT=SPARES/100. 

CAPINV=NM*CUNIT 

SPRINV=CAPINV*SPRPCT 

PV=I. 

IF(RATE.LT..01)GO  TO  2 
PV=RATE/(1.-(1./(1.+RATE)**LIF£)) 

COSTd  )=CAPIMV*PV 

COST(2)=SPRINV*PV 

COST(3)=OPER*NM 

C0ST(4)=KP3MNT*WRATE*NM 

CQST(5)=WRATE*MTTR/3600.*(PASSYR/MCBF) 

COST(6)=0. 

DO  3 J=l,5 

C0ST(6)=C0ST(6)-»-C0ST(J) 

GO  TO  1 
WRITE(4, 101)1 

FORMAT!'  NO  COST  DATA  FOR  OR  NON-EXISTENT  UNIT',13) 

CONTINUE 

RETURN 

END 


SUBROUTINE  MINV 
PURPOSE 

INVERT  A MATRIX 
USAGE 

CALL  MINV{ A,N,D,L,M) 

DESCRIPTION  OF  PARAMETERS 

A - INPUT  MATRIX,  DESTROYED  IN  COMPUTATION  AND  REPLACED  BY 
RESULTANT  INVERSE. 

N - ORDER  OF  MATRIX  A 
D - RESULTANT  DETERMINANT 
L - WORK  VECTOR  OF  LENGTH  K 


MINV  10 
MINV  20 
MINV  30 
MINV  40 
MINV  50 
MINV  60 
MINV  70 
MINV  80 
MINV  90 
MINV  100 
MINV  110 
MINV  120 
MINV  130 
MINV  140 
MINV  150 
MINV  160 
MINV  170 
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M 


WORK  VECTOR  OF  LENGTH  N 


10 

15 


20 


25 


REMARKS 

MATRIX  A MUST  BE  A GENERAL  MATRIX 


SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 


METHOD 

THE  STANDARD  GAUSS-JQRDAN  METHOD  IS  USED.  THE  DETERMINANT 
IS  ALSO  CALCULATED.  A DETERMINANT  OF  ZERO  INDICATES  THAT 
THE  MATRIX  IS  SINGULAR. 


SUBROUTINE  MI NV ( A, N,D, L, M ) 
DIMENSION  A(1  ),L(1),M(1) 


IF  A DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED/  THE 
C IN  COLUMN  1 SHOULD  BE  REMOVED  FROM  THE  DOUBLE  PRECISION 
STATEMENT  WHICH  FOLLOWS. 

DOUBLE  PRECISION  A/D/BIGA/HOLD 

THE  C MUST  ALSO  BE  REMOVED  FROM  DOUBLE  PRECISION  STATEMENTS 
APPEARING  IN  OTHER  ROUTINES  USED  IN  CONJUNCTION  WITH  THIS 
ROUTINE. 

THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 
CONTAIN  DOUBLE  PRECISION  FORTRAN  FUNCTIONS.  ABS  IN  STATEMENT 
10  must  be  CHANGED  TO  DABS. 


SEARCH  FDR  LARGEST  ELEMENT 

D = 1.C 
NK=-N 

DC  90  K=l/N 

NK=NK+N 

L( K)=K 

M(K)=K 

KK=NK+K 

bIGA=A(KK) 

DO  20  J=K/N 
1Z=N*(J-1) 

DO  20  I=K/N 
IJ=IZ+I 

IF(  ARS(EIGA)-  ABS(A(IJ)))  15,20,20 
BIGA=A(IJ) 

L(K)=I 

H(K)=J 

CONTINUE 

INTERCHANGE  ROWS 
J=L(K) 

IF(J-K)  35,35,25 
KI=K-N 


MINV  180 
MINV  190 
MINV  200 
MINV  210 
MINV  220 
MINV  230 
MINV  240 
MINV  250 
MINV  260 
MINV  270 
MINV  280 
MINV  290 
MINV  300 
MINV  310 
MINV  320 
MINV  330 
MINV  340 
MINV  350 
MINV  360 
MINV  370 
MINV  380 
MINV  390 
MINV  400 
MINV  410 
MINV  420 
MINV  430 
MINV  440 
MINV  450 
MINV  460 
MINV  470 
MINV  480 
MINV  490 
MINV  500 
MINV  510 
MINV  520 
MINV  530 
MINV  540 
MINV  550 
MINV  560 
MINV  570 
MINV  580 
MINV  590 
MINV  600 
MINV  610 
MINV  620 
MINV  630 
MINV  640 
MINV  650 
MINV  660 
MINV  670 
MINV  680 
MINV  690 
MINV  700 
MINV  710 
MINV  720 
MINV  730 
MINV  740 
MINV  750 
MINV  760 
MINV  770 
MINV  780 


A-A8 


oooooo  ooo  ooo  nooo  ono 


DO  30  1=1, N 

MINV  790 

KI=K1+N 

MINV  800 

H0LD=-A(KI) 

MINV  810 

J1=KI-K+J 

MINV  820 

A(KI)=A(J1) 

MINV  830 

30 

A(JI)  =H0LD 

MINV  840 
MINV  850 

INTERCHANGE  COLUMNS 

MINV  860 
MINV  870 

35 

I=H(K) 

MINV  880 

IF(I-K)  45,45,38 

MINV  890 

38 

JP=N*U-1) 

MINV  900 

DO  40  J=l,iN 

MINV  910 

JK=NK+J 

MINV  920 

JI=JP+J 

MINV  930 

H0LD=-A(JK) 

MINV  940 

A(JK)=A(JI) 

MINV  950 

40 

A(JI)  =H0LD 

MINV  960 
MINV  970 

DIVIDE  COLUMN  BY  MINUS  PIVOT  (VALUE  OF  PIVOT  ELEMENT  IS 

MINV  980 

CONTAINED  IN  BIGA) 

MINV  990 
MINVIOOO 

45 

IF(BIGA)  48,46,48 

MINVIOIO 

46 

D=0.0 

MINV1020 

RETURN 

MINV1030 

48 

DO  55  1=1, N 

MINV1040 

IF(I-K)  50,55,50 

MINV1050 

50 

1K=NK+I 

MINV1060 

A(IK)=A(IK)/(-bIGA) 

MINV1070 

55 

CONTINUE 

MINV1080 

MINV1090 

REDUCE  MATRIX 

MINVllOO 

MINVlllO 

DO  65  I=1,N 

MINV1120 

IK=NK+I 

MINV1130 

HOLD=A(IK) 

MINV1140 

1J=I“N 

MINV1150 

DO  65  J=1,N 

MINV1160 

lJ=IJ+;\ 

MINV1170 

IF(I-K)  60,65,60 

MINV1180 

60 

IF(J-K)  62,65,62 

MINV1190 

62 

KJ=IJ“I+K 

MINV1200 

A(1J)=H0LD*A(KJ)+A(IJ) 

MINV1210 

65 

CONTINUE 

MINV1220 

MINV1230 

DIVIDE  ROW  BY  PIVOT 

MINV1240 

MINV1250 

KJ= K-N 

MINV1260 

DO  75  J=1,N 

MINV1270 

KJ=KJ+N 

MINV1280 

IF(J-K)  70,75,70 

MINV1290 

70 

A(KJ)=A(KJ)/BIGA 

MINV1300 

75 

CONTINUE 

MINV1310 

MINV1320 

PRODUCT  OF  PIVOTS 

MINV1330 

MINV1340 

D=D*niGA 

MINV1350 

MINV1360 

REPLACE  PIVOT  BY  RECIPROCAL 

MINV1370 

MINV1380 

A(KK)=1.0/BIGA 

MINV1390 

A-49 


non 


flO  CONTINUE 

FINAL  ROW  AND  COLUMN  INTERCHANGE 

K = N 

100  K=(K-1) 

IF(K)  150,150,105 
105  I=L(K) 

IF(I-K)  120,120,108 
10b  JQ=N*(K-1) 

JR=N*(I-1) 

DO  110  J=1,N 

JK=JQ+J 

HOLD=A(JK) 

JI=JR+J 

A(JK)=-A(JI) 

110  A(JI)  =HOLD 
120  J=H(K) 

IF(J-K)  100,100,125 
125  KI=K-N 

DO  130  I=1,N 

KI=KI+N 

HOLD=A(KI) 

JI=KI-K+J 

A(KI)=-A(JI) 

130  A(JI)  =HOLD 
GO-  TO  100 
150  RETURN 
END 


SUBROUTINE  GELG 


PURPOSE 

TO  SOLVE 


A GENERAL  SYSTEM  OF  SIMULTANEOUS  LINEAR  EQUATIONS. 


GELG 
GELG 
GELG 
GELG 
GELG 
GELG 
GELG 
GELG 
GELG 
GELG 

THE  M BY  N MATRIX  OF  RIGHT  HAND  SIDES.  (DESTROYED)GELG 
ON  RETURN  R CONTAINS  THE  SOLUTION  OF  THE  E QUATIONS. GELG 


USAGE 

CALL  GELG{k,A,M,f.,EPS,IER) 

DESCRIPTION  OF  PARAMETERS 
R 


A 

M 

N 

EPS 

lEP 


THE  M BY  M COEFFICIENT  MATRIX.  (DESTROYED) 
THE  NUMBER  OF  EQUATIONS  IN  THE  SYSTEM. 

THE  NUMBER  OF  RIGHT  HAND  SIDE  VECTORS. 

AN  INPUT  CONSTANT  WHICH  IS  USED  AS  RELATIVE 
TOLERANCE  FOR  TEST  ON  LOSS  OF  SIGNIFICANCE. 
RESULTING  ERROR  PARAMETER  CODED  AS  FOLLOWS 


IER=0 

IER=-1 


IER=K  - 


REMARKS 

INPUT  MATRICES  R 


NO  ERROR, 

NO  RESULT  BECAUSE  OF  M LESS  THAN  1 OH 
PIVOT  ELEMENT  AT  ANY  ELIMINATION  STEP 
EQUAL  TO  0, 

WARNING  DUE  TO  POSSIBLE  LOSS  OF  SIGNIFI- 
CANCE INDICATED  AT  ELIMINATION  STEP  K+1/ 

WHERE  PIVOT  ELEMENT  WAS  LESS  THAN  OR 
EQUAL  TO  THE  INTERNAL  TOLERANCE  EPS  TIMES  GELG 
ABSOLUTELY  GREATEST  ELEMENT  OF  MATRIX  A.  GELG 

GELG 

GELG 

AND  A ARE  ASSUMED  TO  BE  STORED  COLUMNWISE  GELG 


MINV1400 
MINV1410 
MINV1420 
MINV1430 
MINV1440 
MINV1450 
MINV1460 
MINV1470 
MINV1480 
MINV1490 
MINV1500 
MINV1510 
MINV1520 
MINV1530 
MINV1540 
MINV1550 
MINV1560 
MINV1570 
MINV1580 
MINV1590 
MINV1600 
MINV1610 
MINV1620 
MINV1630 
MINV1640 
MINV1650 
MINV1660 
MINV1670 
MINV1680 
GELG  10 
,GELG  20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 


GELG 

GELG 

GELG 

GELG 

GELG 

GELG 

GELG 

GELG 

GELG 

GELG 

GELG 

GELG 

GELG 


A-50 


ooo  ooonoo  oooo 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


IN  RESP,  M*M  SUCCESSIVE  STORAGE  LOCATIONS.  ON  RETURN 

SOLUTION  MATRIX  R IS  STORED  COLUMNWISE  TOO. 

THE  PROCEDURE  GIVES  RESULTS  IF  THE  NUMBER  OF  EQUATIONS  M IS 
GREATER  THAN  0 AND  PIVOT  ELEMENTS  AT  ALL  ELIMINATION  STEPS 
ARE  DIFFERENT  FROM  0.  HOWEVER  WARNING  IER=K  - IF  GIVEN  - 
INDICATES  POSSIBLE  LOSS  OF  SIGNIFICANCE.  IN  CASE  OF  A WELL 
SCALED  MATRIX  A AND  APPROPRIATE  TOLERANCE  EPS,  IER=K  MAY  BE 
INTERPRETED  THAT  MATRIX  A HAS  THE  RANK  K.  NO  WARNING  IS 
GIVEN  IN  CASE  M=l, 


SUBROUTINES 

NONE 


AND  FUNCTION  SUBPROGRAMS  REQUIRED 


METHOD 

SOLUTION  IS  DONE  BY  MEANS  OF  GAUSS-ELIMINATION  WITH 
COMPLETE  PIVOTING. 


SUBROUTINE  G£LG( R, A,M, N, EPS, lER) 


DIMENSION  A(1),P(1) 
IF(M)23,23,1 


SEA 

lER 

PIV 

MM= 

NM= 

DO 

T3= 

IF( 

PIV 

1=L 

CON 

TQL 

A(I 


RCH  FOR  GREATEST  ELEMENT  IN  MATRIX  A 
= 0 
= 0. 

M*M 

N*M 

3 L=1,MM 
Af'S(A(L)  ) 

Tb-PIV)3,3,2 
= TB 

TIRUE 
= £PS*P1  V' 

) IS  PIVOT  ELEMENT.  PIV  CONTAINS  THE  ABSOLUTE  VALUE  OF  A(I). 


START  ELIMINATION  LOOP 
l.ST  = l 

DO  17  K=1,M 

TEST  ON  SINGULARITY 
IF(PIV)23,23,4 
4 IF(I£R)7,5,7 
t)  IF(PIV-TOL)  6,6,7 

6 IER=K-1 

7 PIVI=1./A(I) 

J=(T-1)/M 

I=I-J*M-K 

J=J+1-K 

I+K  IS  ROW-INDEX,  J+K  COLUMN-INDEX  OF  PIVOT  ELEMENT 

PIVOT  ROW  REDUCTIQ?v  AND  ROW  INTERCHANGE  IN  RIGHT  HAND  SIDE  R 

DO  8 L=K,NM,M 

LL=L-^I 

k(LL)=R(L) 

8 R(L)=T8 


GELG  330 
GELG  340 
GELG  350 
GELG  360 
GELG  370 
GELG  380 
GELG  390 
GELG  400 
GELG  410 
GELG  420 
GELG  430 
GELG  440 
GELG  450 
GELG  460 
GELG  470 
GELG  480 
GELG  490 
.GELG  500 
GELG  510 
GELG  520 
GELG  530 
GELG  540 
GELG  550 
GELG  560 
GELG  570 
GELG  580 
GELG  590 
GELG  600 
GELG  610 
GELG  620 
GELG  63C 
GELG  640 
GELG  650 
GELG  6o0 
GELG  670 
GELG  680 
GELG  690 
GELG  700 
GELG  710 
GELG  720 
GELG  730 
GELG  740 
GELG  750 
GELG  760 
GELG  770 
GELG  780 
GELG  790 
GELG  800 
GELG  810 
GELG  820 
GELG  830 
GELG  840 
GELG  850 
GELG  860 
GELG  870 
GELG  880 
GELG  890 
GELG  900 
GELG  910 
GELG  920 
GELG  930 


A-51 


oooo  noorj  no  nono 


IS  ELIMINATION  TERMINATED 
IF(K-M)9,10,18 

COLUMN  INTERCHANGE  IN  MATRIX  A 
9 LEND=LST-^V-K 
IE(J)12,12, 10 

10  II=J*M 

DO  11  L=LST,LEND 
TB=A(L) 

LL=L+II 

A(L)=ACLL) 

11  A(LL)=TB 

RON  INTERCHANGE  AND  PIVOT  ROW  REDUCTION  IN  MATRIX  A 

12  DO  13  L=L3T,MM,M 
!.f  =I  +I 

TB=PIV1*A(LL) 

A(LL J=A(L) 

13  A(L)=Tb 

SAVE  COLUMN  INTERCHANGE  INFORMATION 
A(LST)=J 

element  kEDUCTIQN  and  NEXT  PIVOT  SEARCH 
PIV=0. 

LST=LST+1 
J = 0 

DO  16  II=LST,LEND 
PIVI=-A(II) 

I3T=11+M 

J=J+1 

DO  15  L=IST,MM,M 
LL=L-J 

A(L)=A(L)+PIVI*A(LL) 

TB=ABS(A(L) ) 

IF(Tb-FIV )15,15,14 

14  PIV=TB 
1=L 

15  CONTINUE 

DO  16  L=K,NM,M 
LL=L+J 

16  K(LL)=R(LL)+PIVI*R(L) 

17  LST=LST+M 

END  OF  ELIMINATION  LOOP 


BACK  SUBSTITUTION  AND  BACK  INTERCHANGE 
lb  1F(M-1)23, 22,19 
19  IST  = MM-»-M 
LST=M-«-l 
DU  21  1=2, M 
II=LST-I 
IST=IST-LST 
L=IST-M 
L=A(L)-*-.5 
DO  21  J=II,NM,M 
TB=R(J) 

LL=J 

DO  20  K=IST,MM,M 
LL=LL+1 


GELG  940 
GELG  950 
GELG  960 
GELG  970 
GELG  980 
GELG  990 
GELGIOOO 
GELGlOlO 
GELG1020 
GELG1030 
GELG1040 
GELG1050 
GELG1060 
GELG1070 
GEXGIOBO 
GELG1090 
GELGllOO 
GELGlllO 
GELG1120 
GELG1130 
GELG1140 
GELG1150 
GELG1160 
GELG1170 
GELG1180 
GELG1190 
GELG1200 
GELG1210 
GELG1220 
GELG1230 
GELG1240 
GELG1250 
GELG1260 
GELG1270 
GELG1280 
GELG1290 
GELG1300 
GELGi310 
GELG1320 
GELG1330 
GELG1340 
GELG1350 
GELG1360 
GELG1370 
GELG1380 
GELG1390 
GELG1400 
GELG1410 
GELG1420 
GELG1430 
GELG1440 
GELG1450 
GELG1460 
GELG1470 
GELG1480 
GELG1490 
GELG1500 
GELG1510 
GELG1520 
GELG1530 
GELG1540 
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20 

TB=TB-A(K)*R(LL) 

GELG1550 

K=J+L 

GELG1560 

R(J)=R(K) 

GELG1570 

21 

R(K)=TB 

GELG1580 

22 

RETURN 

GELG1590 

GELG1600 

GELG1610 

ERROR  RETURN 

GELG1620 

23 

IEP=-1 

GELG1630 

RETURN 

GELG1640 

END 

GELG1650 

A-53/A-54 
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